Precision Kaon and Hadron Physics with KLOE 

F. Bossi, E. De Lucia, J. Lee-Franzini, S. Miscetti, M. Palutan 
and KLOE Collaboration^) 

Laboratory Nazionali di Frascati dell'INFN, Via E. Fermi 40, 1-00044 Frascati, Italia 



Summary. — We describe the KLOE detector at DA"1>NE, the Frascati 0-factory, 
and its physics program. We begin with a brief description of the detector design 
and operation. Kaon physics is a major topic of investigation with KLOE thanks in 
part to the unique availability of pure Ks,Kl, K beams at a 0-factory. We have 
measured all significant branching ratios of all kaon species, the Kl and K ± lifetimes 
and the K — > -k form factor's t dependence. From the measurements we verify the 
validity of Cabibbo unitarity and lepton universality. We have studied properties 
of light scalar and pseudoscalar mesons with unprecedented accuracy. We have 
measured the e + e~ — >7r + 7r _ cross section necessary for computing the major part of 
the hadronic contribution to the muon anomaly. The methods employed in all the 
above measurements as well as the (j> leptonic width, precision mass measurements 
and searches for forbidden or extremely rare decays of kaons and ^-mesons are 
described. The impact of our results on flavor and hadron physics to date, as well 
as an outlook for further improvement in the near future, are discussed. 
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1. Introduction 

Late in 1989 the Istituto Italiano di Fisica Nucleare, INFN, decided to construct an 
e+e~ collider meant to operate around 1020 McV, the mass of the ^-meson. The <£-meson 
decays mostly to kaons, neutral and charged, in pairs. Its production cross section peaks 
at about 3 /ibarns. Even with a modest luminosity, £=100 pb /s, hundreds of kaon 
pairs are produced per second. This collider, called a ^-factory and christened DA<I>NE, 
is located at the Laboratori Nazionali di Frascati, LNF, INFN's high energy physics 
laboratory near Rome [1]. 

11. DA$NE. - The DA$NE layout is shown in fig. 1. The heart of the collider are 
two storage "rings" in which 120 bunches of electrons and positrons are stored. Each 
bunch collides with its counterpart once per turn, minimizing the perturbation of each 
beam on the other. Electrons and positrons are injected in the rings at final energy, ~510 
MeV. Electrons are accelerated to final energy in the Linac, accumulated and cooled in 
the accumulator and transferred to a single bunch in the ring. Positrons require first 
accelerating electrons to about 250 MeV to an intermediate station in the Linac, where 
positrons are created. The positrons then follow the same processing as electrons. 

Since the stored beam intensities decay rapidly in time, the cycle is repeated sev- 
eral times per hour. A collider is characterized by its luminosity L defined by: event 
rate=£x cross section, beam lifetime and beam induced radiation background. The latter 
is large for short lifetimes, which also result in the average luminosity being smaller than 
the peak value, in a word, lower event yield and large, variable background. DA$NE is 
housed in the ADONE building at LNF, see fig. 2. 

DA<i>NE runs mostly at center of mass energy W = m<p w 1019.45 MeV. The <j> meson 
decays dominantly to charged kaon pairs (49%), neutral kaon pairs (34%), pn (15%), and 
i]-f (1.3%). A <p factory is thus a copious source of tagged and monochromatic kaons, 
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both neutral and charged. In most of the following we will assume that all analysis 
are carried out in the <fi- meson center of mass. While this is in fact what we do, it is 
appropriate to recall that the electron and positron collide at an angle of 7T-0.025 radians. 
The 0-mesons produced in the collisions therefore move in the laboratory system toward 
the center of the storage rings with a momentum of about 13 MeV corresponding to 
/30^O.O15, 70~1.OOO1. K mesons from 0-decay are therefore not monochromatic in the 
laboratory, fig. 3. 

All discussion to follow refers to a system of coordinates with the x-axis in the hor- 
izontal plane, toward the center of DA<J?NE, the y axis vertical, pointing upwards and 
the z-axis bisecting the angle of the two beam lines. The momentum of neutral kaons 
varies between 104 and 116 MeV and is a single valued function of the angle between 
the kaon momentum in the laboratory and the 0-momentum, i.e. the x-axis. Knowledge 
of the kaon direction to a few degrees allows to return to the 0-meson center of mass, 
fig. 3. The mean charged kaon momentum is 127 MeV. 




Fig. 2. - The Adone building. 
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Fig. 3. - Laboratory momentum vs angle to the z-axis for K -mesons. 



1'2. KLOE. - A state-of-the-art detector is needed to collect data at DA$NE, from 
which physics results can be extracted. In summer 1991, the KLOE Collaboration was 
initiated proposing a detector and a program. The KLOE detector was not formally ap- 
proved nor funded for another couple of years, because the authorities had not envisioned 
fully the magnitude of KLOE. By mid Summer 1998 however, KLOE was designed, con- 
structed, completely tested and, complete with all of its electronics for signal processing, 
event gathering and transmission, was practicing with Cosmic Rays. Between Christmas 
and 1999 New Year, KLOE was moved from its own assembly hall onto the DA$NE's 
South Interaction Region. 

The detailed KLOE saga, with its concomitant joys and travails, could be an enter- 
taining tale that should be written up somewhere. Our purpose in this paper, is instead, 
chapter 2, to extract and describe those properties of the detector, chapter 3, its per- 
formance, and chapter 4, the special kaon "beams" available at DA$NE, that made the 
physics measurements possible. 

In the subsequent seven years KLOE, after beginning with a trickle of beam, by 2006 
had produced data of such precision that the kaon, the first fundamental particle to 
introduce the concept of "flavor" into our current way of thinking, got its definitive 21st 
century portrait re-mapped with high precision, 60 years after its discovery. These are 
the subjects covered in the next chapters: chapter 5 details the new measurements of 
kaon physical parameters and decays; chapter 6, discusses the KLOE determination of 
V us of the CKM mixing matrix , unitarity and lepton universality tests, as well as search 
for new physics; chapter 7 is on quantum interferometry , and chapter 8 is on tests of 
CPT. Concurrently, capitalizing on the profuse production of light mesons resulting from 
electromagnetic transitions from the </>-mesons KLOE has improved scores of light meson 
properties often by orders of magnitude and these are discussed in chapter 9. Besides 
neutral kaon decays, symmetry tests have also been performed using light meson and 
charged kaon decays, see sections 9T.2 9'4.2, and 5'3. Finally, utilizing initial state 
radiation KLOE has measured the dipion cross section, ow, over the whole range from 
the 7T7T threshold to the cj> mass energy, see chapter 10. 

Most KLOE data were taken during two running periods. The 2001-2002 run yielded 
450 pb _1 of data, corresponding to approximately 140 million tagged K$ decays, 230 
million tagged Kl decays, and 340 million tagged K ± decays. Analysis of this data 
set is essentially complete and most of the results are described in this report. KLOE 
also took data in 2004-2005, after completion of a number of DA$NE modifications 
intended to increase the luminosity. About 2000 pb _1 were collected, giving a fivefold 
increase in statistics. In total, KLOE has collected during the years 1999 - 2005 ^2.5 
fb _1 luminosity. Analysis of the 2004-2005 data is very advanced and prognosis of the 
improvements in accuracy in the measurements that are statistically limited are included 
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either in the specific sections or in the final chapter 11. 



2. Detector 



The KLOE detector has a giant inner core, approximately 125 cubic meters, see fig. 4, 
and was built in Frascati in the KLOE assembly building, see fig. 5. It is on the same 




scale and complexity as the "general purpose detectors" of its time that were operating 
at the world's foremost laboratories such as LEP at CERN. If we consider that LEP 
operated at about one hundred times DA<I>NE's energy, producing much more energetic 
particles and with correspondingly greater multiplicities, KLOE's large size would seem 
unwarranted. We shall see however that while KLOE's complexity was necessitated by 
its intent of being a definitive high precision experiment, its size, instead, was dictated by 
the quirk of one of the two neutral kaons, called the Kl-, of living an extraordinarily long 
time (for a short lived particles) , about 51 nanoseconds. Kaons from 0-mesons decaying 
at rest, travel at approximately one fifth of the speed of light. This low speed is a great 
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bonus for many reasons. However the mean path travelled by a K^-meson, Xl — ^Pct is 
3.4 m. Thus if we want a detector which can catch 1 — l/e=63% of all decaying neutral 
long-lived kaons, it must have a radius of some three and a half meters. We compromise 
on two meters and catch some 40% of those decays. Less than that in fact, since some 
volume is inevitably lost around the beam interaction point, IP, and some more is lost 
at the outer edge, in order to be able to recognize the nature of the decays close to it. 

The major components of KLOE that enclose the decay volume around the IP arc, 
going outwards radially, a large tracking device to reconstruct the trajectory of charged 
particles, the Drift Chamber DC, surrounded by a hermetic calorimeter to measure the 
energy and the entry point of photons, the Barrel EMC and two End-cap EMC. The 
tracking chamber and the calorimeter are immersed in a mostly axial magnetic field of 
0.52 T, provided by a superconducting coil enclosed in its iron yoke. 




Fig. 5. - KLOE Assembly Hall. The KLOE control room is on the top floor. The Adone cupola 
is just visible behind. 



2T. Drift Chamber. - In the design of the KLOE tracking chamber one is constrained 
by the desire to catch all the charged secondary products from a decay and measure 
their properties with great precision, without jeopardizing the need to achieve the same 
for the neutral secondary products. Recalling that when a charged particle travels in 
a medium, a gas, it undergoes repeated collisions with electrons, releasing part of their 
energy, resulting in a string of ion pairs created marking its passage. Thus when the 
particle enters a drift chamber, DC (a gas filled chamber traversed by a large number of 
wires, some kept at +2000 volt, anode wires, some at ground), electrons of the ion pairs 
created along its trajectory drift to the positive voltage wires and from an avalanche 
multiplication mechanism a detectable signal appears at the wire's end. KLOE custom 
designed sophisticated electronics that are mounted at the wire's end which senses the 
tiny signals and provide us with the drift charge and time measurement. Furthermore, 
the total integrated charge collected at the wire's end gives information on the energy 
released from the initial particle. 
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Furthermore, particles of different mass travelling in the same medium and with 
the same momenta have different energy releases: for the same p they have different /3 
and approximately the specific ionization dE/dx oc 1//3 2 . Therefore this information 
is suitable to perform particle identification (PID). Finally, measuring the drift time 
locates the distance of the track from the anode wires which is then used to reconstruct 
the trajectory of the particles. 

A caveat, the actual path of the particle is altered by multiple scattering, the multiple 
scattering angle being oc l/\/X , where X is the radiation length which is roughly 
proportional to l/Z 2 , in gr/cm 2 (and one must not forget the density). Finally, while 
moving towards the anode, the drifting ionization undergoes diffusion, the amount of 
which is a function of the density of the gas mixture of the DC. 

After judicious balancing, in order to avoid the Kl^K$ regenerations (Kg generated 
from the strong interaction of a Kl with the traversed medium that might simulate CP- 
violating decays), a helium based gas mixture was chosen to minimize multiple scattering 
and density, despite the fact diffusion of drifting ionization is larger in He than certain 
other gases. The average value of the radiation length in the DC volume is Xq ~ 900 m, 
including the wire's contribution. 

Again, for the mechanical structure low-Z and low density materials were used to 
minimize Kl regeneration, tracks' multiple coulomb scattering and the absorption of 
photons before they can reach the calorimeter. To achieve this the chamber is constructed 
out of a carbon fiber composite. The chamber volume is delimited by an outer cylinder 
of 2 m radius, an inner cylinder of 25 cm radius and closed by two 9 mm thick annular 
cndplates 9.76 m radius of curvature and a 25 cm radius hole in the center. Due to their 
large size, the endplates were fabricated in two halves then glued together. About 52,000 
holes were drilled on each cndplate for the insertion of the feedthroughs holding the DC 
wires. The total axial load on the end plates, about 3500 kg, would deform the end 
plates by several mm at the end of stringing. The deformation problem was solved by 
using rings of rectangular cross-section which surround the L-shaped plate rims, made 
of the same material. The ring is used to apply a tangential pull to the spherical plate 
via 48 tensioning screws made of titanium alloy pulling on tapped aluminum inserts in 
the L-shaped flange. The net result of the pull is to approximate a complete sphere 
beyond the outer edge of the plate which greatly reduces the deformation. Each screw 
was instrumented with two strain gauges, calibrated in order to compensate for bending 
effects on the screws, so that the pulling force was always known. At the end of stringing 
the deviation from the spherical shape of the endplates was measured to be less than 1 
mm and the hole positions to be accurate at the level of 60 /xm. The KLOE DC is indeed 
a tour de force achievement thereafter copied in other experiments (albeit never to the 
same large dimensions). 

The KLOE DC [2] has ~52,000 wires, about f2,500 of them are sense wires, the 
remaining ones shape the electric field in 12,500 almost square cells located on cylindrical 
layers around the z-axis (defined by the bisector of the two crossing beams). Particles 
from the <f> decays are produced with small momenta and therefore track density is much 
higher at small radii. This motivated the choice of having the 12 innermost layers made 
up of 2 x 2 cm 2 cells and the remaining 48 layers with 3x3 cm 2 cells. For a particle 
crossing the entire chamber we get dozens of space points from which we reconstruct a 
track: the track of a pion from Kg — > ir + n~ decay has ~60 space points in average. 
An all stereo geometry allows the measurement of the z coordinate. Fig. 6 shows the 
wire geometry during the DC construction as illuminated by light. In this particular 
arrangement, shown in fig. 7 left, the wires form an angle with the z-axis, the stereo 
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angle e. Consecutive layers have alternating signs of this angle and its absolute value 
increases from 60 to 150 mrad with the layer radius while keeping fixed for all the layers 
the difference between R p and Rq distances (1.5 cm), defined in fig. 7 left. 




Fig. 6. - DC all stereo geometry. 



This design results in a uniform filling of the sensitive volume with almost-square drift 
cells, with shape slowly changing along z, which increases the isotropy of the tracking 
efficiency. 

To extract the space position from the measured drift time of the incident particle, we 
need to know the space-time relation describing the cell response. The whole DC can be 
described using 232 space-time relations only, parametrized in terms of the angles (3 and cb 
defined in fig. 7 right. The (3 angle characterizes the geometry of the cell, directly related 
to the electric field responsible of the avalanche multiplication mechanism. The <j> angle 
gives the orientation of the particle trajectory in the cell's reference frame, defined in the 
transverse plane and with origin in the sense wire of the cell. The presence of a signal on 
a wire and its drift time and distance is called hit. Once we have the hits, we can proceed 
with the event reconstruction using the wire geometry, space-time relations and the map 
of the magnetic field. This procedure follows three steps: (i) pattern recognition, (ii) 
track fit, and (iii) vertex fit. The pattern recognition associates hits close in space to 
form track candidates and gives a first estimate of the track parameters. Then track 
fit provides the final values of these parameters minimizing a \ 2 function based on the 
difference between the fit and the expected drift distances, residuals, as evaluated from a 
realistic knowledge of the cells' response (measured drift times and space-time relation). 
This is an iterative procedure because the cells' response depend on the track parameters. 
Finally the vertex fit searches for possible primary and secondary vertices, on the basis 
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Fig. 7. - Left: Wire geometry with the definition of stereo angle e between the wire of length L 
and the 2-axis. Right: Definition of f3 and <j> characterizing the shape of the cell and the angle 
of the incident track. 



of the distance of closest approach between tracks. 

Besides the drift time and position information, hits in the DC are also associated to 
a charge measurement directly related to the energy released in the cell from the incident 
particle. The amount of collected charge depends on the particle's track length in the 
crossed cell. The energy deposited in a sample of finite thickness fluctuates and exhibits 
a Landau distribution, with a long tail at high energies. Therefore to improve the energy 
resolution, the energy deposited is sampled many times for each track and a truncated 
mean technique is used. To use the charge information for particle identification purposes 
it is necessary to normalize the charge to the track length and define the specific ionization 
dE/dx. If the energy is measured in counts, we have on average 15 counts/cm for pions 
and muons from K^ 2 and K^ 2 decays, with momentum within 180 and 300 MeV. Instead 
K ± tracks exhibit a larger number, ^120 counts/cm. 

The drift chamber provides tracking in three dimensions with resolution in the bending 
plane better than 200 /xm, resolution on the z measurement of ~2 mm and of ~1 mm 
on the decay vertex position. The particle's momentum is determined from the curvature 
of its trajectory in the magnetic field with a fractional accuracy better than ^0.4% 
for polar angles larger than 45°. The calibration of the absolute momentum scale has 
a fractional accuracy of few 10~ 4 using several two- and three-body decays (e + e~ — > 
e + e - ,e + e~ — > Kl — ► ttIv 1 Kl — > 7r + 7r~7r°, K ± — > fj^i/,K^ — > tt^tt ) covering a 

wide momentum range. Deviations from nominal values of invariant mass, missing mass 
or momentum of the charged decay particle in the rest frame of the decaying particle 
are used as benchmarks for the calibration procedure. The resolution on the specific 
ionization measurement obtained using two- and three-body decays and Bhabha 
scattering events, is <~ 6% for pions and ~ 7% for electrons [3], using 60 truncated 
samples. 

The DC system's complexity, as well as the demands placed upon it requires the 
response of each cell to an incoming particle to be known accurately. The cell response 
is a function of its geometry and of the voltage applied to the wires. Thus the position 
of the DC, mechanical structure as well as wires, had to be surveyed, in situ, together 
with the other KLOE components. 
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The Slow Control monitoring system (sect. 2'6) fully embedded in the Data Acqui- 
sition process (sect. 2'5), sets and registers the voltages applied to the wires, sends the 
information to the read-out electronics and controls its proper functioning [4]. More- 
over the gas mixture parameters, namely composition, temperature and pressure, enter 
directly in the space-time relations used in reconstructing the particle trajectory, to ex- 
tract the space position from the measured drift time of the incident particle. The Slow 
Control System also monitors the gas parameters. 

To ensure the stability in time of the DC performance, the system has to be calibrated 
periodically by acquiring samples of cosmic-ray muons suitable for the measurement of 
the ^200 different space-time relations. The calibration program, incorporated into the 
KLOE online system, automatically starts at the beginning of each run and selects about 
80,000 cosmic-ray events. These events are tracked using the existing space-time relations 
and the average value of the residuals for hits in the central part of the cells is monitored. 
If the residuals exceed 40 /mi then about 300,000 cosmic-ray events are collected, with 
about 30 Hz rate, and a new set of calibration constants is obtained. At the beginning 
of each data-taking period a complete calibration of the space-time relations is needed 
together with the measurement of the DC time offsets, the latter from some 10 7 cosmic- 
ray events. Finally, during data taking the DC performances are monitored using selected 
samples of Bhabha scattering events. 

2'2. Electromagnetic Calorimeter. - The neutral decay products of the K^s, its short- 
lived partner K$ (lifetime <~ 89 picoseconds) and of those of their charged brethren, the 
K^'s, which we need to detect, are photons, either directly produced or from the decay 
of neutral pions into two photons. It is important to recall that photons for energies 
larger than a few MeV, as in most of our cases, interact with matter by electron positron 
pair production that in turn radiate photons when traversing matter, a process quite 
similar to pair-production. Thus, a photon, or an electron for that matter, travelling 
through dense, high Z matter, repeatedly undergoes radiation and pair production in 
cascade, until all its energy is expended in a so called shower of e + e~ and photons. 
From detection of these particles and their energy one measures the original photon's (or 
electron's) energy. This is done with a "calorimeter", which also usually pinpoints the 
position of the electromagnetic cascade (or shower) as well. 

The KLOE calorimeter, EMC, has been designed and built to satisfy a set of stringent 
requirements such as providing a hermetic detection of low energy photons (from 20 to 
500 MeV) with high efficiency, a reasonable energy resolution and an excellent time 
resolution to reconstruct the vertex for the Kl neutral decays. The calorimeter response 
has also to be fast since its signals are used to provide the main trigger of the events. 

The chosen solution for the EMC is that of a sampling calorimeter, composed of 
lead passive layers, accelerating the showering process, and of scintillating fiber sensing 
layers. It is made of cladded 1 mm scintillating fibers sandwiched between 0.5-mm- 
thick lead foils. The foils are imprinted with grooves wide enough to accommodate the 
fibers and some epoxy, without compressing the fibers. This precaution prevents damage 
to the fiber-cladding interface. The epoxy around the fibers also provides structural 
strength and removes light travelling in the cladding. In more detail, the basic structure 
of 0.5 mm, grooved Pb layers between which are embedded 1 mm diameter scintillating 
fibers are positioned at the corners of almost equilateral triangles. ^200 such layers are 
stacked, glued, and pressed, resulting in a bulk material. The resulting composite has 
a fibers:lead:glue volume ratio of 48:42:10 which corresponds to a density of 5 g/cm 3 , 
an equivalent radiation length X of 1.5 cm and an electromagnetic sampling fraction of 
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<~13%. The large ratio of active material to radiator, the frequent signal sampling and 
the special fiber arrangement result in a factor \/2 improvement in energy resolution with 
respect to calorimeter with slabs of equivalent scintillatordead ratio. The chosen fibers, 
Kuraray SCSF-81 and Pol. Hi. Tech. 00046, are blue emitting scintillating fibers with a 
fast emission time whose larger component is described by an exponential distribution 
with a r of 2.5 ns. In addition, the special care in design and assembly of the Pb-fiber 
composite ensures that the light propagates along the fiber in a single mode, resulting in 
a greatly reduced spread of the light arrival time at the fiber ends. Propagation speed 
in the fiber is <~ 17 cm/ns. These fibers produce more than 5xl0 3 photons for one MeV 
of deposited energy by ionizing particles. 3% of the light reaches the fiber ends. Fibers 
transmit 50% of the light over 2m(/i~ 4m) . 

This material is shaped into modules 23 cm thick (~15Ao). 24 modules of trapezoidal 
cross section are arranged in azimuth to form the calorimeter barrel, aligned with the 
beams and surrounding the DC. An additional 32 modules, square or rectangular in 
cross section, are wrapped around each of the pole pieces of the magnet yoke to form 
the endcaps, which hermetically close the calorimeter up to —98% of iir, see fig. 8. 
The unobstructed solid-angle coverage of the calorimeter as viewed from the origin is 
-94% [5]. 




Fig. 8. - EMC viewed before DC insertion. Note end-caps attached to pole pieces seen in profile. 



The fibers run parallel to the axis of the detector in the barrel, vertically in the end- 
caps, and are read out at both ends viewed by lucite light guides of area of 4.4x4.4 cm 2 . 
The other end of the guides are coupled to the photocathodes of fine mesh photomulti- 
plicrs of 1.5 inches diameter with quantum efficiency of —20%. The fiber length of 4 m 
corresponds to an average photoelectron yield of 35 p. e. /layer for a minimum ionizing 
particle crossing the calorimeter at the fiber center. The photomultiplicr tubes at each 
fiber's end transform the light into electric pulses. Their amplitudes, Ai, is proportional 
to the amount of energy deposited while the recorded times, U, are related to the time 
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of flight of the particle. For each cell the position of the readout elements and the dif- 
ference of the arrival times at the two fiber ends determine the shower position with a 
0(1) cm accuracy. The full stack of 23 cm (15 Xq) fully contains a shower of 500 MeV. 
4880 photomultiplier tubes view the 24+32+32 modules into which the calorimeter is 
subdivided for a grand total of 15,000 km used optical fibers. 

As a first step, the reconstruction program makes the average of time and energy 
of the recorded U and Ai for the two sides of each cell and compute the hit position. 
Corrections for attenuation length, energy scale, time offsets and light propagation speed 
are taken into account at this stage. A clustering procedure then groups together nearby 
clumps of energy deposition and calculates the average quantities over all the partici- 
pating cells. The calibration constants to transform ti and Ai from raw quantities to 
time in nanoseconds and energy in MeV are evaluated with dedicated on-line and off-line 
algorithms. 

The energy calibration starts by a first equalization in cell response to minimum 
ionizing particles (m.i.p) at calorimeter center and by determining the attenuation length 
of each single cell with a dedicated cosmic ray trigger. This is done before the start of 
each long data taking period. The determination of the absolute energy scale in MeV 
relies instead on a monochromatic source of 510 MeV photons: the e + e~ — > 77 sample. 
This last calibration is routinely carried out each 200- 400 nb _1 of collected luminosity. 

For the timing, the relative time offsets of each channel, T?, related to cable lengths 
and electronic delays and the light velocity in the fibers are evaluated every few days 
with high momentum cosmic rays selected with DC information. An iterative procedure 
uses the extrapolation of the tracks in the calorimeter to minimize the residuals between 
the expected and measured times in each cell. A precision of few tens of picoseconds is 
obtained for these offsets. 

To run at the design luminosity, DA<i>NE operated with a bunch-crossing period equal 
to the machine radio frequency (RF) period, Trf = 2.715 ns. Due to the spread of the 
particle's arrival times, the trigger is not able to identify the bunch crossing related 
to each event, which has to be determined offline. The start signal for the electronics 
devoted to time measurement (TDC) is obtained by phase-locking the trigger level- 1 
signal (sect. 2'4) to an RF replica, with a clock that has a period of 4 x Trf- The 
calorimeter times are then related to the time of flight of the particle from IP to the 
calorimeter, T to f, by T c i — T to f + 6c — N^Trf, where 6c is a single number accounting 
for the overall electronic offset and cable delay, and Nt, c is the number of bunch-crossings 
needed to generate the TDC start. The values of 6c and Trf are determined for each 
data taking run with e + e~ — > 77 events by looking at the distribution of T c ; — R c i jc where 
Rd/c is the expected time of flight: well separated peaks correspond to different values 
of Nb c . We arbitrarily define 6c as the position of the peak with the largest statistics, 
and determine Trf from the peaks' distances. Both quantities are evaluated with a 
precision better than 4 ps for 200 nb _1 of integrated luminosity. This measurement of 
Trf allows us to set the absolute calorimeter time scale to better than 0.1 %. During 
offline processing, to allow the cluster times to be related to the particle time of flight, we 
subtract 5c and determine, on an event by event basis, the "global" event start time T , 
i.e. the quantity Ni, c Trf- A starting value for all analysis is evaluated by assuming that 
the earliest cluster in the event is due to a photon coming from the interaction point. 
Further corrections are analysis dependent. 

The high photon yield and the frequent sampling enable cluster energies to be mea- 
sured with a resolution of (Je/E = 5.7% /y/E (GeV), as determined with the help of the 
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DC using radiative Bhabha events. The absolute time resolution a t = 57 ps/^E (GcV) 
is dominated by photoelectron statistics, which is well parametrized by the energy scal- 
ing law. A constant term of 140 ps has to be added in quadrature, as determined from 
e+e~ — > 77, radiative <j> decays and cj) — > -k + tt~it data control samples. This constant 
term is shared between a channel by channel uncorrelated term and a common term to 
all channels. The uncorrelated term is mostly due to the calorimeter calibration while 
the common term is related to the uncertainty of the event To, arising from the DA$NE 
bunch-length and from the jitter in the trigger phase-locking to the machine RF. By mea- 
suring the average and the difference of T c i — R c i/c for the two photons in <f> — > 7r + 7r~7r° 
events, we estimate a similar contribution of <~ 100 ps for the two terms. Cluster posi- 
tions are measured with resolutions of 1.3 cm in the coordinate transverse to the fibers, 
and, by timing, of 1.2 cm/ y/E (GeV) in the longitudinal coordinate. These performances 
enable the 27 vertex in Kl — ► -k + tt~tt decays to be localized with a w 2 cm along the 
Kl line of flight, as reconstructed from the tagging Kg decay. Incidentally, the thin lead 
layers used and the high photon yield allows high reconstruction efficiency for low energy 
photons. 

2'3. Beam Pipe and Quadrupole Calorimeter. - At DA$NE the mean Kg, Kl and K ± 
decay path lengths are As = 0.59 cm, \l = 3.4 m and A± = 95 cm. To observe rare Kg 
decays and K^Kg interference with no background from Kg — * Kl regeneration, a decay 
volume about the interaction point with r > 15 As must remain in vacuum. The beam 
pipe, see fig. 9, surrounds the interaction point, with a sphere of 20 cm inner diameter, 
with walls 0.5 mm thick made of a Be-Al sintered compound. This sphere provides a 
vacuum path ~ 10 x the K° amplitude decay length, effectively avoiding all Kg — > Kl 
regeneration. Permanent quadrupoles for beam focusing are inside the apparatus at 




Fig. 9. - Beam pipe at the KLOE interaction region. Just visible on either side are the final 
focus quadrupole ends. 



a distance of 46 cm from IP, and surrounding the beam pipe. The quadrupoles of 
this low-/? insertion are equipped with two lead/scintillating-tile calorimeters, QCAL, 
of <~ 5X° thickness, with the purpose of detecting photons that would be otherwise 
absorbed on the quadrupoles [6]. Specifically, the QCAL main task is to identify and 
reject photons from Kl — ► 37r° decays when selecting the CP violating Kl — ► 2tt° events. 
Each calorimeter consists of a sampling structure of lead and scintillator tiles arranged 
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in 16 azimuthal sectors. The readout is performed by plastic optical fibers, Kuraray 
Yll, which shift light from blue to green (wavelength shifter, WLS), coupled to mesh 
photomultipliers. The special arrangement of WLS fibers allows also the measurement 
of the longitudinal z coordinate by time differences. Although the tiles are assembled 
in a way which maximizes efficiency for Kl photons, a high efficiency is in fact also 
obtained for photons coming from the IP. This allows us to extend the EMC coverage of 
the solid angle down to 008 0=0.99 for prompt decays. The time resolution is related to 
the light yield and emission time/photoelectron of the Kuraray Yll fibers. For a m.i.p. 
crossing the whole stack of 16 tiles, we had obtained during tests a light yield of ~ 50 
photoelectrons for entry points close to the PM's side. Considering that a m.i.p. deposits 
an energy equivalent to that of a 75 MeV photon, this corresponds to <~ 0.5 p.e./MeV 
and in turn implies a time resolution for electromagnetic showers of 205 ps/-\/E(GeV). 
However during running, the presence of the B-fi cld reduc es the effective light yield by a 
factor 1.4, thus worsening the timing to 240 ps/ ^/E(GcV), which however is high enough 
to yield an efficiency above 90% for photon energies greater than 20 MeV. 

2'4. Trigger. - Event rates at DA$NE arc high; at the maximum luminosity of 10 32 
cm _2 s _1 up to 3000 s~ x and 30,000 Bhabha s _1 are produced within the KLOE accep- 
tance window. The trigger design [7] has been optimized to retain all 4> decays. Moreover, 
all Bhabha and 77 events produced at large polar angles are accepted for detector mon- 
itoring and calibration, as well as a downscaled fraction of cosmic-ray particles, which 
cross the detector at a rate of <~ 3000 Hz. Finally, the trigger provides efficient rejec- 
tion on the two main sources of background: small angle Bhabha events, and particle 
lost from the DA$NE beams, resulting in very high photon and electron fluxes in the 
interaction region. 

The trigger electronic logic system rapidly recognizes topologies and energy releases 
of interest, processing the signals coming both from EMC and DC. Since </> decay events 
have a relatively high multiplicity, they can be efficiently selected by the EMC trigger 
by requiring two isolated energy deposits (trigger sectors) in the calorimeter above a 
threshold of 50 MeV in the barrel and 150 MeV in the endcaps. Events with only two 
fired sectors in the same endcap are rejected, because this topology is dominated by 
machine background. Moreover, events with charged particles in the final state give 
a larger number of hit wires in the DC than do background events. The DC trigger 
employs this information, requiring the presence of ^15 hits in the DC within a time 
window of 250 ns from beam crossing. An event satisfying at least one of the two above 
conditions generates a first level trigger, Tl, which is produced with minimal delay and 
is synchronized with the DA$NE master clock. The Tl signal initiates conversion in the 
front-end electronics modules, which are subsequently read out following a fixed time 
interval of ~ 2/^s, which is ultimately driven by the typical drift distances travelled by 
electrons in the DC cells. In case of DC trigger, a validation of the first level decision is 
required by asking for ^120 hits within a time window. 

The KLOE trigger also implements logic to flag cosmic-ray events, which are recog- 
nized by the presence of two energy deposits above 30 MeV in the outermost calorimeter 
layers. For most data taking, such events were rejected after partial reconstruction by 
an online software filter. 

2'5. Data acquisition and processing . - At a luminosity of ^100 /ib _1 /s, the trigger 
rate is ^2000 Hz. Of these ~500 Hz are from <fi decays or Bhabha scattering. All digitized 
signals from the calorimeter, drift chamber, calibration and monitoring systems are fed 
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via fiber optics from the electronics, sitting on platforms mounted on the KLOE iron 
yoke, to the computers in the next building. The Data Acquisition architecture has been 
designed to sustain a throughput of 50 Mbytes/s all along this path. In the near-by 
control room one can see real time displays of reconstructed charged tracks traversing 
the DC and photons depositing energy clusters in the EMC. The display of two events 
is shown in fig. 10. 

The events are reconstructed, meaning that charged tracks and energy deposits that 
occurred at the same instant in time on the various subcomponents of KLOE are asso- 
ciated together, properly labelled and packed into a unique file. To optimize the event 
reconstruction we use continuously updated calibration constants for the time and en- 
ergy scales. Furthermore, background hits are acquired using a random trigger, which are 
recorded in order to be used in the simulation of KLOE events (MC). The MC simulation 
program reproduces accurately the detector geometry, the response of the apparatus and 
the trigger requirements , but of course cannot foretell the moment to moment variation 
of the background. 

All raw data, background, reconstructed and MC events are recorded on tape. KLOE 
has a huge tape library of ~800 TB capacity [4]. At reconstruction time, events are 
classified into broad categories, or "streams". Data summary tapes are produced with a 
separate procedure. 




Fig. 10. - KsKl events in KLOE. Left: both Kl and Ks decay to charged pions tt tt~ . Color 
patches indicate energy in the calorimeter. Right: 5 neutral particles. 71 trough 74 are showers 
from 4 photons from Ks — * — >47. The iCE,-mesons reaches the calorimeter before decaying. 
It literally crashes in the EMC, producing a huge energy release. 



Any of the KLOE analyses, even after event reconstruction is quite complex. It 
includes making data summary tapes and extensive modelling of the detector response. 
This requires generation of large number of Monte Carlo (MC) events for each process 
under study. "Offline" processes run in a farm of computers of ~ 150 kSPECint2K CPU 
power, all residing in the dedicated computer center. Both data and MC summaries 
(~80 TB) are cached on disk for analysis [8]. 

2'6. Slow control system. - KLOE has to be operated 24 hours a day, 7 days a week. 
Therefore it has to be kept under continuous control, both to guarantee efficient data 
taking, and for safety reasons. Parameters such as high and low voltage settings, the 
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Table I. - Summary table of the EMC and DC performances. 



temperature of the electronics components, or the status of the drift chamber gas system 
are set and monitored by a dedicated system fully embedded in the Data Acquisition 
process: the Slow Control. 

This system also records some machine parameters, and provides to the DA$NE team 
information about KLOE including our on line measurements of the luminosity and of 
the background levels. 

3. Reconstruction Performances 

KLOE provides a continuous monitoring of the machine working point, providing 
feedback to DA$NE continuously as well. The most important parameters are the beam 
energies and crossing angle, which are obtained from the analysis of Bhabha scattering 
events with polar angles above 45 degrees. The average value of the center-of-mass 
energy is evaluated online during data taking with a precision of ~ 50 kcV for each 200 
nb _1 of integrated luminosity. This determination is further refined with offline analysis 
to achieve a precision of ~ 20 keV, as discussed later. The position of the e + e~ primary 
vertex, with coordinates Xpy, Ypy, and Zpy, is reconstructed run- by- run from the same 
sample of Bhabha events. Xpy and Ypy are determined with typical accuracy of about 
10 fim, and have widths L(X) and L(Y) which are about 1 mm and few tens of microns, 
respectively. Zpy is also reconstructed online with 100 — 200 /im accuracy, but it has a 
natural width of 12 — 14 mm, determined by the bunch length itself. 

3'1. Overview of EMC and DC performances. - We note that while KLOE is a general 
purpose apparatus with good performances of the two main detectors, the DC and the 
EMC, its real strength manifests when joining the reconstruction capability of both detec- 
tors. To give an overall view at a glance, we first summarize in table I the parametrized 
resolution for neutral particles (mainly photons) and for charged tracks (mainly elec- 
trons/pions) for prompt decays and K^,K^. Resolution on invariant masses are also re- 
ported in the same table as well as the reconstruction capability for neutral and charged 
vertices. 

In KLOE, the single particle reconstruction efficiency can be determined from data 
control samples, as well as from Monte Carlo simulation. As an example, the photon 
reconstruction efficiency as a function of the photon energy can be estimated from <f> — ► 
7T + 7r~7r° events, where a redundant determination of the decay photon kinematics is 
possible using the charged pion momenta as measured from the DC. For charged particles 
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produced at the interaction point, the track efficiency can be measured from Kg — > 
-k + tt~ events selected with the requirement of having at least one track fulfilling the two- 
body decay kinematics. A high reconstruction efficiency is observed for both detectors 
on data and Monte Carlo. For each specific analysis, the event selection efficiency is 
evaluated from Monte Carlo. To take into account data-MC difference in the cluster/track 
efficiencies, the calculation is usually performed by weighting each particle in the final 
state with the ratio of data/MC reconstruction efficiencies. 

For particle identification, PID, we start extrapolating tracks on the calorimeter sur- 
face and associating tracks to clusters. The identification of photons is based on clusters 
not associated to any track and on the difference between expected and reconstructed 
time of flight, ToF, for straight trajectories. The PID of charged tracks, which at the con- 
sidered energy are electrons, pions and muons, is based on a combined usage of EMC and 
DC, by looking at ToF, E/p and cluster shape. At low momenta, P < 200 MeV, the ToF 
is usually the winning handle while at higher momenta variables based on cluster shape 
are used. Often these information are combined together in a likelihood identification 
variable [9] or using a neural network [10]. To identify charged kaons we can use instead 
dE/dx, given the excellent separation provided by the specific ionization measurement 
between kaons and the other charged particles produced in the apparatus. 

3'2. Absolute energy scale. - An improved determination of the center-of-mass energy, 
W, is obtained for each run by fitting the e+e~ invariant-mass distribution for Bhabha 
events to a Monte Carlo generated function, which includes radiative effects. Initial state 
radiation (ISR), where one or both initial colliding particles radiate a photon before 
interacting, affects the e + e~ center-of-mass energy and therefore the final state invariant 
mass. The ISR is mostly collinear to the beam, and in general is not detected. MC 
Bhabha events were generated using the Babayaga [11] event generator, which accounts 
for both final and initial state radiation. 

The absolute energy scale is calibrated by measuring the visible cross section for the 
(j) — > KsKl process. The cross section peak is fitted to a theoretical function [12], 
which depends on the <j> parameters, takes into account the effect of ISR, and includes 
the interference with the p(770) and the w(782) mesons. The <j) mass, total width, and 
peak cross section are the only free parameters of the fit, the p(770) and the w(782) 
parameters being fixed. The results of the fit to the data are shown in fig. 11. The 
4> mass value obtained from the fit is — 1019.329 ± 0.011 MeV, which has to be 
compared with = 1019.483 ± 0.011 ± 0.025 MeV, measured by CMD-2 at VEPP-2M 
[13]. For CMD-2 the absolute energy scale is accurately determined from beam induced 
resonant depolarization, so that we use the ratio m^ MD /m^ LOE = 1.00015 to correct 
our determination of the center of mass energy. This corresponds to a shift in the value 
of W of -150 keV. 

3'3. Luminosity . - The tagging capability of the calorimeter trigger allowed us to 
perform fast luminosity measurements, providing a robust estimate of the DA<I>NE lumi- 
nosity with few per cent accuracy. This was achieved by counting Bhabha events in the 
acceptance of the calorimeter barrel. For this purpose, we selected events with two trig- 
ger sectors fired on the barrel, imposing a higher discrimination threshold than the one 
used in the standard acquisition chain. The tight angular selection, which reduces the 
available statistics by a factor of — 6, was motivated by the need to keep under control 
the overwhelming machine background, which is concentrated at small angle. Additional 
cuts on the position of the fired sectors and on the time difference between the two sig- 
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Fig. 11. - Top: Cross section for e + e — > <j> — * KsKl as a function of the center-of-mass energy 
(line, fit). Bottom: Fit residuals. 



nals were added to improve rejection on the residual background, from cosmic-ray events. 
When operating at £=100 ^ib _1 /s, the luminosity measurement was updated every 15 s, 
with a statistical error of ~ 3%. Fig. 12 top shows the monthly integrated luminosity by 
KLOE since January 2001, while in the bottom the year integrated luminosity is shown. 
These data have been collected almost exclusively at the <f) peak. In early 2006, KLOE 
also collected ~200 pb _1 of data at yfs = 1000 MeV, at which a(e+e- -> cj> -> tt+tt^tt ) 
drops to 5% of its peak value. Finally, it is important to note that already in year 
2000 KLOE collected 20 pb _1 , producing quite a few publications which improved the 
knowledge in both hadronic and kaon physics [14, 15, 16, 17, 18]. 

A more accurate measurement of the integrated luminosity is performed offline [19], 
using Bhabha events in the polar angle range 55° < 6 < 125°, the so-called VLAB 
events. The effective cross section for these events, ~ 430 nb, is large enough to reduce 
the statistical error at a negligible level. The luminosity is obtained by counting the 
number of VLAB candidates, AVlab, and normalizing it to the effective Bhabha cross 
section, <Jy^ AB , obtained from MC simulation, after subtraction of the background, 5bk B : 

(1) [Cdt = ^^(l-5 hks ). 

J °Vlab 

The precision of the measurement depends on the correct inclusion of higher-order 
terms in computing the Bhabha cross section. For this purpose we use Babayaga, 
which include the QED radiative corrections in the framework of the parton-shower 
method. The quoted precision is 0.5%, although an updated version of this generator, 
BabayagaONLO [20] , has been released recently. The new predicted cross section decreased 
by 0.7% and the theoretical systematic uncertainty improved from 0.5% to 0.1%. 

The VLAB events are selected with requirements on variables that are well repro- 
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duced by the KLOE MC. The acceptance cut on the electron and positron polar angle, 
55° < # + ,_ < 125° is based on the calorimeter clusters, while the momentum selection, 
p+ : - > 400 MeV, is based on the DC determination. The background from /i+/i~(7), 
7r + 7r _ (7), and tt + tt~tt° events is well below 1% and is subtracted. All selection effi- 
ciencies (trigger, cluster and track reconstruction) are greater than 99%, as obtained 
from MC and confirmed with data control samples. Finally, corrections are applied on 
a run-by-run basis to follow the small variations in the center-of-mass energy and in the 
detector calibration. In the end, we quote a 0.3% uncertainty in the determination of 
the acceptance due to experimental effects. 

4. Kaon beams 

The neutral kaon pair from <fi — > K°K° is in a pure J PC — \~~ state. Therefore the 
initial two-kaon state can be written, in the (/>-rest frame, as 

\KK,t = 0) = {\K°,p)\K°,-p) - \K°,p)\K°,-p)) /V2 
= {\K s ,p)\K L ,-p) - \K L ,p)\K s ,-p))/V2, 

where the identity holds even without assuming CPT invariance. Detection of a Kg 
thus signals the presence of, "tags" , a Kl and vice versa. Thus at DA<I>NE we have pure 
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K$ and Kl beams of precisely known momenta (event by event) and flux, which can be 
used to measure absolute Ks and Kl branching ratios. In particular DA<f>NE produces 
the only true pure K$ beam and the only Kl beam of known momentum. A K$ beam 
permits studies of suppressed K$ decays without overwhelming background from the Kl 
component. A Kl beam allows lifetime measurements. Similar arguments hold for K + 
and K~ as well, although it is not hard to produce pure, monochromatic charged kaon 
beams. In the following sections we briefly discuss the criteria adopted to tag the K$, 
Kl, and charged kaon beams. 

4'1. Ks beam. - Because of the exceptional timing capabilities of the EMC and the 
slowness of the kaons from 0-decays, we use a ToF technique to tag ifg-mesons in a 
unique way. As already mentioned, neutral kaons have a velocity (3 — 0.2162 in the 4> 
frame and /3=0.193 to 0.239 in the lab, and 60% of them reach the calorimeter in a time 
>31 ns. T^i-mesons interact in the calorimeter with an energy release up to their mass, 
497 MeV. For an energy release of 100 MeV we have a time resolution of ~0.3 ns, i.e. 
about a ~1% accuracy on Kl velocity, shown in fig. 13. 



Fig. 13. - Kl velocity in <^-frame from time of flight. 



In KLOE, we call a Kl interaction in the EMC a "K-crash", and we use it as a tag 
for the Ks- Note that the above error on (3{Kl) corresponds to an error on the kaon 
energy of ~0.25 MeV or ^1 MeV for its momentum, with just one event. Adding the 
information about the position of the energy release, the direction of the Kl line of flight 
is determined with ~ 1° angular accuracy. Using the previous value of Kl momentum, 
and the run-average value of the <fi momentum as evaluated from Bhabha events, we 
determine the Ks momentum from ps = p$ — pl, with the same accuracy as for Kl ■ 

4'2. Kl beam. - All of the produced Ks mesons decay within few centimeters from the 
interaction point. Therefore, in principle, all of them can be detected in the apparatus. 
On the contrary, Kl's can travel several meters before decaying, thus escaping detection. 
From the above, it follows that a pure sample of nearly monochromatic -ftT^'s can be 
selected by identification of Ks decays. The best performances in terms of efficiency 
and momentum resolution are obtained by selecting Ks — > tt + tt~ decays. These are 
identified as two-track vertices close to the e + e~ interaction point, with invariant mass 
in a 5 MeV acceptance window around the nominal kaon mass. The Ks — » 7r + 7r~ decay 
provides also a good measurement of Kl momentum, pl — V4> ~ Vs, where p^ is again 
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the <f> momentum as determined from Bhabha scattering events. The resolution is ^0.8 
MeV on the momentum, while the Kl flight direction is determined within ~ 2°. 

More interestingly, events with both a Kg — > ir + n~ decay and a K-crash identified 
allow us to measure the DA$NE beam energy spread. This is extracted by comparing 
for each event the two independent determinations of center-of-mass energy that are 
evaluated using the Ks and Kl momenta, respectively. Using this technique we find a 
beam energy spread of ~220 keV, in agreement with the machine model expectation. 

For many KLOE measurements, it is of fundamental importance to accurately de- 
termine the exact amount of lost/retained decays, or, in other words, determining the 
Kl decay fiducial volume. For decays involving two charged particles this is obtained 
relatively easily by the observation of the vertex of the two related tracks in the drift 
chamber, measured with an accuracy of ~ 3 mm. The task is more difficult for totally 
neutral decay channels, but can still be achieved using the excellent timing accuracy of 
the EMC. As an example, we can locate the Kl^tt°tt decay point, the "neutral vertex", 
with accuracies of 0(2) cm, as illustrated in fig. 14, left. Each photon from the Kl decay 




L K (cm) 

Fig. 14. - Left: D is the Kl—*tt tt — >4y decay point. The Kl path is obtained from the time 
for each photon to arrive at the calorimeter, i(I — »A), and the Ks direction. Right: Resolution 
<j(Lk) on the determination of the Kl decay length using photon vertices, as a function of Lk, 
in Kl — > 7r + 7r~7r° events. The contributions from the uncertainties on the point of photon 
incidence on the calorimeter, the cluster time, and the Kl flight direction are shown as the 
dot-dashed, dashed, and dotted lines, respectively. 



point (D) defines a time-of-flight triangle: the first side is the segment from the interac- 
tion point to the Kl decay vertex (ID); the second is the segment from the Kl decay 
vertex to the centroid of the calorimeter cluster (DA); and the third is the segment from 
the interaction point to the cluster centroid (IA). The Kl direction is initially known, 
because the Kl decay is tagged. The photon vertex position is specified by the distance 
ID, which is determined from 

I A 2 + ID 2 - 2IA ■ ID = DA 2 

(3) ID/f3 K + DA = ct 7 , 

where i 7 is the cluster time and (3k is the Kl velocity. The position of the interaction 
point is obtained by backward extrapolation along the Kg flight path. The Kl decay 
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vertex position is evaluated for each neutral cluster separately, and the energy-weighted 
average of the values for each cluster is then taken as its final determination. The 
accuracy in the location of the photon vertex has been studied using Kl — > 7r + 7r~7r° 
decays, in which the decay position can be independently determined using clusters and 
tracks, with much greater precision in the latter case. The dependence of the position 
resolution on decay distance (Lk) is illustrated in fig. 14, right. 

4'3. K ± beam. - The </>-meson decays 50% of the times into K + K~ pairs; these 
are quasi anti-collinear in the laboratory, due to the small crossing angle of the e + e~ 
beams, with momentum ranging between 120 and 133 MeV/c and decay length A± <~ 95 
cm. While moving towards the DC volume, the average K ± momentum decreases to 
~100 MeV/c due to the energy loss in the beam pipe and DC inner wall materials. As 
for neutral kaons, the identification of a K T decay tags a K beam with known flux, 
thus giving the possibility to measure kaon absolute branching ratios. Charged kaons 
are tagged using their two-body decays, K ± — > and K ± — > 7r ± 7r°, accounting 

for ^85% of the total decay channels. These decays are observed as "kinks" in a track 
originating at the IP. The charged kaon track must satisfy 70 < pk± < 130 MeV/c. The 
momentum of the decay particle in the kaon rest frame evaluated using the pion mass, p*, 
is 205 and 236 MeV, for K ± decays to n ± ir° and ^v, respectively. Fig. 15 shows the p* 
distribution obtained assuming the decay particle is a pion. Cutting on the spectrum of 
the decay particle, we select about 1.5 x 10 e K + K~ events/pb -1 of integrated luminosity. 



o.o3 r 
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180 200 220 240 260 



Fig. 15. - Momentum spectrum in the kaon rest frame of the negative charged decay particle, 
assuming the particle has the pion mass for data (dots) and MC (lines). The distribution is 
normalized to unity. The two peaks correspond to pions and muons from K~ — > 7T _ 7r° (205 
MeV/c) and K~ — ► fj,~ > M (236 MeV/c). The muon peak is broadened by the use of the incorrect 
mass. 



5. Kaon physical parameters and decays 



Since the kaon's discovery about sixty years ago, knowledge of its properties, from 
masses, lifetimes, and how often each species decays into to which particles (the so 
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called branching ratios, BR), has been accumulated piecemeal literally over hundreds 
of individual experiments, performed at various laboratories all over the world. Each 
experiment has its own peculiarity, contingent upon the particle beam available and its 
own equipment's acceptance range, which can result in highly accurate measurements 
from a statistical point of view (namely based on thousands of data events) but with 
hidden corrections peculiar to it perhaps unknown even to its authors. The particle data 
group [21] in charge of compiling the results, often adopt a "democratic" way of weighing 
the various results according to their nominal errors and literally applying a scale factor 
to these errors when faced with two inconsistent measurements. This procedure skirts 
the problems of judging the validity of results obtained over time and sometimes ends 
up with a hodge-podge of inconsistent results forced into a mold of dubious scientific 
validity. 

The mission of KLOE is to measure most, if not all, of the properties of the kaon 
system to high accuracy, with a single detector simultaneously. One problem that con- 
sistently plagues the interpretation of older measurements is the lack of clarity about 
accounting for radiative contributions. All of our measurements of kaon decays with 
charged particles in the final state are fully inclusive of radiation. Radiation is auto- 
matically accounted for in the acceptance correction. All our MC generators incorporate 
radiation as described in [22]. 

Because of the availability of tagged and pure kaon beams of known momenta, KLOE 
is the only experiment that can at once measure the complete set of experimental inputs, 
branching ratios, lifetimes and form factor parameters for both charged kaons and long 
lived neutral kaons. In addition KLOE is the only experiment that can measure Kg 
branching ratios at the sub-percent level. 

5T. K a Mass. - Kaon masses are about 500 MeV, very close to one half the </>- 
meson mass, known to an accuracy of 0.019/1019.460 ^2 x 10~ 5 by the use of the g-2 
depolarizing resonances, [21]. The events <fi — > K$Kl offer a unique possibility to obtain 
a precise value of the neutral kaon mass. Indeed we observe that if the <f>- meson is at 
rest the kaon mass can be extracted from the kaon momentum using the relation: 

,, Am K p\ Ap K 2 Ap K 

(4) m K = \ -^-p 2 K ; ~ (3 Z 

m K m z K p K p K 

Since pk — 110 MeV, measuring it at 1% level, well within the KLOE capability, results 
in a measurement of the K° mass better than 0.1%. 50,000 events are enough to reach 
a statistical accuracy of about 1 keV. 

<f> mesons are produced with a momentum along the x axis, p^ = 12.5 MeV at DA<i>NE. 
From the measured momenta of the two pions from K$ — > n + ir~, we measure the K$ 
momentum. The momentum is given by pk l = P4, — Pk s j where p^ is the average 
(j> momentum measured with Bhabha events collected in the same runs. The center of 
mass energy of the KsKl pair (Wkk) is related to the kaon mass rriK, according to: 

W KK {m K ) = yj 2M| + 2E Ks E Kl - 2p Ks ■ p KL 

with 




E Ks = \/p 2 K s + m K 



E K L = \jv 2 K L +m 2 K . 
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On the other hand the collision center of mass energy W is computed from Bhabha 
events, as described in sect. 3. 

Corrections due to ISR have to be taken into account when relating W to Wkk ■ The 
correction function fx(W) has been evaluated using a full detector simulation where the 
radiation from both beams has been implemented, and Wkk is reconstructed as in the 
data. The expression of the radiator function has been taken from [23], including 0{a 2 ) 
corrections. The correction |1 — fx(W)\ is very small below the resonance, corresponding 
to a shift in Wkk of 40 keV. For W above the <f> mass, Wkk increases up to 100 keV. 
The neutral kaon mass is then obtained solving the equation: 

(5) W = f K (W) x W KK (m K ) 

The single event mass resolution is about 430 keV. Contributions to the mass resolution 
are: about 370 keV from experimental resolution, about 220 keV from beam energy 
spread, as measured by KLOE in agreement with machine theory, and about 100 keV 
from ISR. 

The systematic error due to the momentum miscalibration has been evaluated by 
changing the momentum scale in computing the pion momenta. A momentum miscal- 
ibration, 5p/p , translates to a miscalibration on the mass Smx/mx = 0.06<5p/p, in 
agreement with the above qualitative calculation. The momentum scale is obtained 
by using several processes covering a wide momentum range, from 50 to 500 MeV 
(Kl — > n + ir~ir°, Kl — > n£v, (f> — ► tt+tt^tt ). We obtain a fractional accuracy be- 
low 2 x 10~ 4 , in agreement with the estimate obtained using Bhabha's [8], which results 
in a systematic error 5mx of 6 keV. 

The systematic error coming from theoretical uncertainty on the radiator function 
has been evaluated considering the contribution from higher order terms in a. The 
correction function fx(W) has been evaluated by excluding the constant term in the 
0(a 2 ). The corresponding change in fx (W) is 1.3 x 10 -5 , corresponding to a variation 
on mx of 7 keV. Further checks have been made by using the function given in [24]: no 
significant differences were observed. Additional systematics come from the dependence 
of the measured mass from the value of W. They have been evaluated by comparing the 
mass values obtained with data collected at W < 1020 MeV and W > 1021 MeV, where 
the value of fx(W) is more than a factor two larger. The difference between the two 
mass values is mx(W < 1020) — mx{W > 1021) = 9 ± 10 keV, consistent with zero. 

Other sources of systematics are due to the uncertainties on the W calibration, i.e., 
the statistic and systematic error on and on obtained from our fit, discussed 

in sect. 3'2. The total contribution from these sources amounts to a mass uncertainty of 
15 keV. Systematic uncertainties are treated as uncorrelated. The result is [25]: 

(6) m K = 497.583 ± 0.005 stat ± 0.020 syst MeV. 

5'2. Kl lifetime. - Precision measurements of the lifetime of neutral kaons, especially 
in presence of many multibody decay modes, is particularly difficult since it is in general 
not possible to prepare monochromatic beams of neutral particles nor to stop them. 
KLOE enjoys the availability of such monochromatic beams. While the K$ lifetime is 
very well known, (0.8958 ± 0.0005) x lO" 10 s [21], this is not the case for the K L lifetime, 
tl- 

A total of approximately 13 million tagged Kl decays are used for the measurement 
of the four major Kl BRs, as discussed in sect. 5'5.1. Since the geometrical efficiency 
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for the detection of Kl decays depends on tl, so do the values of the four BRs: 
(7) BR(K L -> f)/ER<-°\K L /) = 1 + 0.0128 ns" 1 (r L - if >) , 

where BR^ ' is the value of the branching ratio evaluated for 77° , an assumed value for 
Tl- The four relations defined by eq. 7, together with the condition that the sum of 
all Kl BRs must equal unity, allow determination of the Kl lifetime and the four BR 
values independent of t l . This approach, followed in [26], gives tl = 50.72 ± 0.11 sta t ± 
0.13 sys t-stat i 0.33 sys t ns. 

In addition, the lifetime can be measured by observing the time dependence of the 
decay frequency to a single mode. Because of tagging, the lifetime can also be obtained 
by simply counting the total number of decays in a given time interval from a known 
i^Lsample. This second, independent input comes from our separate analysis of the 
proper decay-time distribution for Kl — > 37r° events, for which the reconstruction effi- 
ciency is high and uniform over a fiducial volume of ~0.4Al. The position of the Kl 
decay position is determined from the photon arrival times, with the technique described 
in sect. 4'2. The Kl proper decay-time, t* , is obtained from the decay length, Ik, divided 
by the Lorentz factor (3jc = PK L /fnK- The Kl momentum is evaluated event-by-event 
using the information from the Kg — ► tt + tt~ tagging decay. About 8.5 million decays are 
observed within the proper-time interval from 6.0 to 24.8 ns. We fit the distribution in 
this interval (fig. 16), and obtain tl — 50.92 ± 0.17 s t a t ± 0.25 sys t ns [27]. The system- 
atic error of 0.49% is presently dominated by the uncertainty on the dependence of the 
tagging efficiency on Ik, and by background subtraction. 

This latter measurement is included together with the results for the four main Kl 
BRs and the Kl lifetime in a fit to determine the Kl BRs and lifetime. We also use 
our measurements of BK(Kl -> 7r+7r-)/BR(i<: L /i3) and BR(K L -> jj)/BR(K l -> 3tt°), 
reported in sections 5'5.4 and 5'5.5, requiring that the seven largest Kl BRs add to 
unity. The only non-KLOE input to the fit is the 2006 PDG ETAFIT result BR(K L -> 
Tr°Tr°)/BR(KL — > 7r + 7r~) = 0.4391 ± 0.0013, based on relative amplitude measurements 
for K — > 7T7T. The results of the fit are presented in table II; the fit gives % 2 /ndf = 0.19/1 
(CL = 66%). The BRs for the KLei and Kl^ decays are determined to within 0.4% 
and 0.5%, respectively. 

5'3. if ± lifetime. - The measurements of the fsT ± lifetime listed in the PDG compi- 
lation, [21], exhibit poor consistency. The PDG fit has a confidence level of 1.5 x 10~ 3 
and the error on the recommended value is enlarged by a scale factor of 2.1. 
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Parameter Value Correlation coefficients 



BR(K e3 ) 


0.4008(15) 












BR(tf„ 3 ) 


0.2699(14) 


-0.31 










BR(3tt ) 


0.1996(20) 


-0.55 


-0.41 








BR(7T + 7T- 


7T°) 0.1261(11) 


-0.01 


-0.14 


-0.47 






BR(7r+7r- 


) 1.964(21) x 10~ 3 


-0.15 


+0.50 


-0.21 


-0.07 




BR(7r°7r ) 


8.49(9) x 10~ 4 


-0.15 


+0.48 


-0.20 


-0.07 


+0.97 


BR( 77 ) 


5.57(8) x 10~ 4 


-0.37 


-0.28 


+0.68 


-0.32 


-0.14 -0.13 


tl 


50.84(23) ns 


+0.16 


+0.22 


-0.14 


-0.26 


+0.11 +0.11 -0.09 


Table II. 


- Final KLOE measurements of 


main Kl 


BRs and 


TL- 





KLOE has measured the decay time for charged kaons in two ways both using events 
tagged by K^ 2 decays [28]. The first method was to obtain the proper time from the 
kaon path length in the DC. The second method was based on the precise measurement 
of the kaon decay time from the arrival times of the photons from kaon decays with a 7r° 
in the final state. These two methods reach comparable accuracy and allow us to cross- 
check systematics. The method relying on the measurement of the charged kaon decay 
length requires the reconstruction of the kaon decay vertex using DC information only. 
The signal is given by a K , moving outwards from the IP in the DC with momentum 
70 < p K ± < 130 MeV/c. The kaon decay vertex (V) has to be in the DC volume defined 
by 40 < yXy+yy < 150 cm, \zy \ < 150 cm. Once the decay vertex has been identified, 
the kaon track is extrapolated backwards to the IP into 2 mm steps, taking into account 
the ionization energy loss dE/dx to evaluate its velocity (3c. Then the kaon proper time 

/l_g2 

(i*) is obtained from: t* = Y^i^i = J2i p- i ^h- The reconstruction and selection 
efficiency of the decay vertex has been evaluated directly from data. The control sample 
is selected using calorimetric information only: neutral vertices are identified looking for 
two clusters in time fired by the photons coming from the 7r° — > 77 decay. The proper 
time distribution is fit between 16 and 30 ns correcting for the measured efficiency 
Resolution effects are taken into account in the convolution with the exponential decay 
function used to fit the t* distribution. The result we have obtained, which is the weighted 
mean between the K + and the K~ lifetimes, is r ± = (12.364+ 0.031 s t a t ± 0.031 sys t) ns. 
The fit window of this method covers 1.1 r . 

The second method relies on the measurement of the kaon decay time and uses charged 
kaon decays with a it in the final state, — > Xir°. In these decays the kaon time of 
flight is obtained from the time of the EMC clusters of the photons from the 7r° decay. We 
require the backward extrapolation of the tagging kaon track to the IP. Then exploiting 
the kinematic closure of the <p — > K + K~ decay, knowing the momentum of 4> and of the 
K + {K~ ) we can build the expected hclicoidal trajectory of the K~{K + ) on the signal side 
(virtual helix) . Stepping along this helix we look for the 7r° — > 77 decay. For each photon 
it is then possible to measure the kaon proper time t* using: t* = (t 7 — — — To) • \J\ — A K 
with f 7 the time of the cluster associated to the photon, r 1 the distance between Vq and 
the position in EMC of the photon cluster, T the production time of the <f> and (3\ the 
beta of the signal kaon. The efficiency has been evaluated directly from data. The control 
sample has been selected using DC information only, selecting the kaon decay vertex in 
the volume defined above. The proper time is fit between 13 and 42 ns, about 2.3 
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Fig. 17. - Charged kaon proper time distribution obtained with the two methods: the kaon 
decay length (Left) and the kaon decay time (Right). The yellow line shows the fit performed 
with the convolution of an exponential decay function and the resolution function. 



correcting for the efficiency and using the convolution of an exponential decay function 
and the resolution function. The weighted mean between the K + and the K~ lifetimes 
gives the result t ± = (12.337 ± 0.030 s t a t ± 0.020 sys t) ns. Given a statistical correlation 
of p = 30.7%, the weighted mean between the two charges and methods is: 



(8) = (12.347 ±0.030) ns 

Both fits are shown in fig. 17. The comparison of K + and K~ lifetimes is a test of 
CPT invariance which guarantees the equality of the decay lifetimes for particle and 
antiparticle. The average of the two methods are r + = (12.325 ± 0.038) ns and r~ = 
(12.374 ± 0.040) ns. From these measurement we obtain: t~ /t+ = 1.004 ±0.004. The 
result agrees well with CPT invariance at the 4 per mil level. 

5'4. Ks decays. - As already mentioned, central to Kg studies is the use of K-crash 
tag (sect. 4'1), thanks to which it has been possible to measure the K$ branching ratios, 
whose values span six orders of magnitude. 

5'4.1. Ks^ir + ir~ (7), ir°ir°. The Ks decays overwhelmingly (99.9%) into two pions: 
7T°7T° and 7r + 7r~. The ratio of the charged decay mode to that of the neutral, 1Z S = 
T(Ks — > it + 7t~ ("f))/T(Ks — * 7r°7r°), is a fundamental parameter of the Ks meson. 
First of all, it provides the BRs for the Ks — ► 7r°7r° and Ks — > 7r + 7r~ with only small 
corrections. The latter BR is a convenient normalization for the BRs of all other Ks 
decays to charged particles. From Rg one can also derive phenomenological parameters 
of the kaon system, such as the differences in magnitude and phase of the isospin 7 = 0,2 
7T7r scattering amplitudes. Finally, this ratio together with the corresponding one from 
the Kl determine the amount of direct CP violation in K — > tttt transitions. Prior 
to KLOE, the fractional error on 1Z S was 1.2%, as a result of an average on various 
measurements, each of ~5% accuracy. This averaging was somehow questionable, since 
the various experiments did not clearly describe their procedure for handling radiative 
events. During the last six years we performed two precise measurements of 1Z S with 
increasing systematical accuracy, and with proper treatment of radiative corrections. 
The final combined value has a precision of three parts per mil [29] . 
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Given the tag, the Ks — ► 7r + 7r~ (7) events are selected by requiring the presence of two 
tracks of opposite charge with their point of closest approach to the origin inside a small 
fiducial volume around the IP, and with momentum in the range 120 < p < 300 MeV. 
Kg — ► 7r°7r° events are identified by the prompt photon clusters from ir° decays. A 
calorimeter cluster is considered as a prompt photon if its velocity, evaluated by ToF, 
is compatible with (3=1 within the expected resolution. Moreover, it must not be 
associated to any track. To accept a Ks — > 7r°7r° event, three or more prompt photons 
are required, with a minimum energy of 20 MeV, to reduce contamination from machine 
background. 

To reach a high precision, all of the possible systematic effects have been studied in 
great detail, making use of several independent data control samples. Systematic errors 
can arise from imperfections in the detector simulation, limitations in the methods used 
to evaluate the detector efficiencies, and uncertainties in the cross sections and branching 
ratios used to estimate the fraction of background events. Moreover, a small difference 
in the tagging efficiency between charged and neutral decays, due to the different time 
response of the two categories of events, has been taken into account. Actually, the 
e + e _ interaction time, the so called To, is obtained from the fastest particle reaching 
the calorimeter, assuming a velocity (3 — 1 and a straight flight path starting from 
the interaction point (sect. 2'2). This assumption is correct for photons coming from 
Ks— >7r 7r°events, while most of the pions from their charged counterpart arrive ~3 ns 
later: therefore their T is delayed by one RF period, ^2.7 ns. As a consequence of this, 
the velocity of the Kl crash clusters is overestimated by ~10% for Ks^>n + ir~ events. 

Another subtle effect that has to be taken into consideration, is the different energy 
response of the EmC to pions and photons. This, in turn, causes a small difference in the 
trigger efficiencies for the two categories of events (^98.5% for 7r+7r~ , ~100% for 7r°7r°), 
which has to be properly taken into account. 

We obtain TZ^ — 2.2549 ± 0.0054. Using this result, and the measurement of T(Ks — > 
■nev)/T{Ks — ► n + n~(j)) discussed in 5.4.2, we extract the dominant Ks branching ratios 
(see table III). To this end, we exploit unitarity: the sum of the BR's for the irir and ttIv 
modes has been assumed to be equal to one, the remaining decays accounting for less 
than 10~ 4 . 

To extract the value of the phase shift difference Sq — 82 between 1 = and 1 = 2 
amplitudes we use the method described in reference [30]. 
In the isospin limit, we find: 

(9) S - 6 2 = (44.5 ±1.0)°. 

Effects of isospin breaking increase this value by <~ 10°, although the correction is deter- 
mined with a large error (see Ref. [30]). 

5'4.2. Ks — > irev(j). The Ks decays semileptonically less than one per cent of the 
time. To pick out such decays where the event contains an unseen neutrino is non- 
trivial. Yet, KLOE has isolated a very pure sample of ^13,000 semileptonic Ks decays 
and accurately measured the BRs for Ks — ► n + e~D(-f) and Ks — ► -K~e + v(~f) [31]. 

The basic steps in the analysis are: tag Ks decays by the Kl crash, apply a cut 
on the -kit invariant mass (this removes 95% of the Ks — > ir + n~ decays and reduces 
the background-to-signal ratio to about 80:1), impose several geometrical cuts to further 
improve the purity of the sample and, in particular, remove contamination by events with 
early n — > [iv decays. Finally, stringent requirements are imposed on the particle's ToF, 
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which very effectively separates electrons from pions and muons and allows the charge 
of the final state to be assigned. 

Fig. 18 shows the signal peak and the residual background in the distribution of 
Ae p = E m - lss — |pmiss| for the it~ e + v channel, where E m i ss and p mlss are respectively the 
missing energy and momentum at the vertex, evaluated in the signal hypothesis. For 
signal events, the missing particle is a neutrino and Ae p = 0. 

•i K -* it e + v 

t 

: * * DATA 
I ■ I 
i i * Signal 

* Background 

-20 20 40 60 

A £p (MeV) 

Fig. 18. - Ae p distribution for K$ — ► Tr~e + v candidates, showing signal and background. 
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The numbers of Trev decays for each charge state are normalized to the number of 
observed tt + it~ events, resulting in the ratios in the first column of table III. Using the 
result for TZ^ of the previous section (also in table III) , the absolute BRs for K3 — > tttt 
and Ks — > 7re^ reported in the second column of the table are obtained. These ratios 



Mode 


BR(mode)/BR(7r + 7r~) 


BR(mode) 


7T + 7T~ 




(69.196 ±0.024 ±0.045)% 


„0 
7T 7T 


1/(2.2549 ± 0.0054) 


(30.687 ± 0.024 ± 0.045)% 


■K~e + u 


(5.099 ± 0.082 ± 0.039) x 10~ 4 


(3.528 ± 0.057 ± 0.027) x 10" 4 


ir + e~i> 


(5.083 ± 0.073 ± 0.042) x 10" 4 


(3.517 ± 0.050 ± 0.029) x 10~ 4 




(10.19 ± 0.11 ± 0.07) x 10" 4 


(7.046 ± 0.076 ± 0.051) x 10^ 4 



Table III. - KLOE measurements of Ks branching ratios. 



give also the first measurement of the semileptonic charge asymmetry for the K$: 

(10) A S = nKs^«-e+v e )-r(K s ^«+e-V e ) = 1Q _ 3 
y ' SUM v ; 

The comparison of As with the asymmetry At, for Kr, decays allows test of the CP 
and CPT symmetries to be performed. Assuming CPT invariance, As =At, ^2Ree~ 
3x 1CU 3 , where e gives the CP impurity of the Ks, Kl mass eigenstates due to CP 
violation in AS=2 transitions (for an exact definition of e see sect. 8T). To evaluate As 
and At, without making any assumptions on CPT symmetry, we consider the standard 
decomposition of the semileptonic amplitudes [32] : 

A(K°^7r-l+is) =A (l-y) , 
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A(K° -> tt+1- v ) = A* {1 + y*)(x+ - x_y , 
A(K°^w + l-u) =A* {l + y*) , 

(11) A{K {) = A {1 - y){x+ + x-) , 

where x + describes the violation of the AS = AQ rule in CPT conserving (violat- 
ing) decay amplitudes, and y parametrizes CPT violation for AS" = AQ transitions. 
The difference between the charge asymmetries, 

(12) A s - A L = 4 (Rc6 + Rc:r_) 

signals CPT violation either in the mass matrix (6 term, see sect. 8T) or in the decay 
amplitudes with AS ^ AQ (Rea;_ term). The sum of the asymmetries, 

(13) A S + A L =4 (Rce-Rey), 

is related to CP violation in the mass matrix (e term) and to CPT violation in the 
A5* = AQ decay amplitude (y term). 

Using Al = (3.34 ± 0.07) x 10~ 3 [21] and our As measurement (eq. 10), we obtain 
from eq. 12 

(14) Rex_ +Re6= (-0.5 ± 2.5) x 10~ 3 . 

Current knowledge of these two parameters is dominated by results from CPLEAR [33] : 
the error on Re<5 is 3 x 10 -4 and that on Rex_ is 10~ 2 . Using Re<5 = (3.0 ± 3.3 sta t ± 
0.6 syst ) x 10~ 4 from CPLEAR, we obtain: 

(15) Rc.t_ = (-0.8 ± 2.5) x 10~ 3 , 

thus improving on the error of Rex- by a factor of four. 
From eq. 13 we obtain 

(16) Rce-Rcy = (1.2 ±2.5) x 10" 3 . 

We calculate Ree using values from [21] that are obtained without assuming CPT con- 
servation: Roe = |e| x cos0 e = (1.62 ± 0.04) x 10~ 3 . Subtracting this value from eq. 16, 
we find 

(17) Rey= (0.4 ±2.5) x 10~ 3 , 

which has precision comparable to that (3x 10~ 3 ) obtained from the unitarity relation 
by CPLEAR [34]. 

5'4.3. Ks~^jj- A precise measurement of the Ks — > 77 decay rate is an important test 
of Chiral Perturbation Theory (ChPT) predictions. The decay amplitude of Ks — > 77 
has been evaluated at leading order of ChPT [35], 0(p 4 ), providing a precise estimate of 
BR(K S -» 77) = 2.1 x 10~ 6 , with 3% uncertainty This estimate is - 30% lower with 
respect to the latest determination from NA48 [36] , thus suggesting relevant contributions 
from higher order corrections. 
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We measured the K$ — * 77 rate using 1.9 fb _1 of integrated luminosity. A MC 
background sample of comparable statistics and a large MC signal sample (equivalent 
to ~ 50 fb _1 of collisions) have been used in the analysis. After K-crash tag, we select 
~ 700 million events, out of which ~ 1900 are K$ — > 77 decays. 

The signal is selected by requiring exactly two prompt photons, with an efficiency of 
~ 83%. After photon counting, the background composition is dominated by Ks — * 27r°, 
with two photons undetected by the EMC. This background is strongly reduced, with 
negligible signal loss, by vetoing the events with photons reaching the QCAL in a 5 ns 
coincidence window with respect to the event-To. At the end of this selection, we are left 
with 157 x 10 3 events and a signal over background ratio S/B=l/80. Further background 
reduction is obtained by performing a kinematic fit to the event in the signal hypothesis, 
and using as constraints the total 4-momentum conservation, the kaon mass and the 
photon ToF. We retain events with x 2 < 20, reaching S/B^ 1/3 with an efficiency on 
signal of ~ 63%. After this cut we count the signal events by fitting the 2-D distribution 
of the photon invariant mass, M 77 , and of the photon opening angle in the K$ center 
of mass, #* 7 , and using MC signal and background shapes. We count AT (77) = 711 ± 35 
signal events out of 2740 events, with x 2 /dof = 854/826, corresponding to a probability 
P(x 2 ) ~ 24%. In fig. 19, the observed distributions of cos#* 7 and M 77 are compared 
with the simulated signal and background shapes as weighted by the fit procedure. To 
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Fig. 19. - Distributions of cos#* 7 (top) and A/ 77 (bottom) at the end of the analysis chain; black 
points are data, solid line is the sum of MC signal and background shapes as weighted by the 
fit. 



get the BR(i^5 — > 77), signal events are normalized to the number of Ks — > 27r° decays 
observed in the same data sample. 

These are selected, after K-crash tag, by requiring three to five prompt photons. 
An extensive study on systematics effects, which affect both the efficiency correction 
and the signal counting, has been carried out. As an example, the MC efficiency of 
the prompt photon selection has been corrected for possible differences in the photon 
detection efficiency between data and MC, and a systematic error has been evaluated 
which accounts for residual imperfections. Moreover, a control sample of Kl — > 77 events 
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decaying close to the IP has been selected, which has been used to check systematics 
from the kinematic fit selection and from data-MC differences in the EMC energy scale. 
Systematic errors on signal counting have been evaluated by: (i) performing a set of fits 
to the M 77 — #* 7 distribution in regions with an enhanced signal content, to evaluate 
the error induced by the uncertainty in the background shape, (ii) repeating the fit with 
different bin-size and applying an energy scale correction a factor of two greater with 
respect to that measured with the Kl — > 77 control sample. 
Finally, we obtain [37] 



which differs by 3ct from the previous best determination. Our result is also consistent 
with 0(p 4 ) ChPT prediction. 

5'4.4. Ks^>3ir°. The decay Kg — > 3tt° violates CP invariance. The parameter 77000, 
defined as the ratio of Kg to Kl decay amplitudes, can be written as 77000 = A(Kg — > 
37T°) / A(K l — > 37T°) = e + €qqq, where e quantifies the Kg CP impurity and e' 000 is due 
to a direct CP-violating term. Since we expect e' 000 << e [38], it follows that 77000 ~ e - 
In the Standard Model (SM), therefore, BR(K S -> 3tt°) ~ 1.9 x 1CT 9 to an accuracy of 
a few %, making the direct observation of this decay quite a challenge. 

We performed a search for the Kg — > 37r° decay using 450 pb _1 of integrated lumi- 
nosity. A MC background sample ~ 2.5 times larger than data and a high statistics MC 
signal sample have been used in this analysis as well. After K-crash tag, we select the 
signal by requiring six prompt photon clusters, and no tracks from the IP, which is useful 
to reject background from Kg — > n + ir~ with an early Kl — > 37r° decay. At this stage 
of the analysis, the residual background is dominated by Kg 7r°7r° events with two 
spurious clusters from shower fragments {splitting) or accidental coincidence with clusters 
produced by machine activity. Further background rejection is obtained by performing 
a kinematic fit to the event in the signal hypothesis, and cutting at x 2 /d.o.f < 3. 

After the previous rejection cuts, the signal search is performed by defining two x 2 - 
like discriminating variables, (3 — (2, optimized for Kg — > 37r° and Kg — > 2ir° selection, 
respectively. (2 is defined by selecting the four out of six photons which provide the best 
agreement in terms of kinematic variables with the Kg — > 2ir° hypothesis. (3 is defined 
as 



where Ami = mi — m^a is the difference between the nominal ir° mass and the invariant 
masses of the zth photon pair, chosen among the six clusters of the event. £ 3 is close 
to zero for a Kg — > 3n° event, and is expected to be large for a six-photon background 
event. A signal box is defined in the (3 — C2 plane, while side-bands are used to check 
the MC prediction for the background. At the end of the analysis chain we observe two 
candidates in the signal box, with an expected background of 3.1 ± 0.8 ± 0.4. The signal 
selection efficiency, after tagging, is ~ 24%. To obtain an upper limit for BK(Kg — ► 37r°), 
we normalize the number of signal candidates to Kg — ► 7r°7r° decays observed in the same 
data sample. For this purpose, we select events with three to five prompt photon clusters. 
Finally, we obtain [39] 



(18) 



BR{K S - 77) = (2.26 ± 0.12 stat ± 0.06 syst ) x 10 
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(19) 



BR(K S -» 3tt°) < 1.2 x 10 



-7 



at 90% CL 
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which is a factor of <~ 6 improvement with respect to the best previous limit [40]. This 
limit on the BR can be directly translated into a limit on 1 77000 1 : 



This result can be visualized in the complex plane as a circle of radius 0.018 centered at 
zero in the Re(77ooo), Ini(r/ooo) plane, and is a 2.5 times better improvement with respect 
to the result of [40] . 

5'4.5. Ks^>e + er . This decay is a flavor-changing neutral current process, suppressed 
in the SM, with an amplitude dominated by the two photon intermediate state. Using 
ChPT to C(p 4 ), one obtains the prediction BR(Ks^>e + e~) ~2xl0~ 14 [41]. A value 
significantly higher could indicate new physics. Prior to KLOE the best experimental 
limit on this decay was set by CPLEAR, BR(K s ^>e+e-) < 1.4 x 10~ 7 at 90% CL [42]. 

We performed a search for the Ks — > e + e~ decay using 1.9 fb _1 of integrated lumi- 
nosity. A MC background sample of comparable statistics and a large MC signal sample 
have been also used in the analysis. After K-crash tag, we search for the signal by requir- 
ing two tracks of opposite charge originating near the IP. The two tracks are required to 
have an invariant mass M ee , evaluated in the electron hypothesis, in a <~ 20 MeV window 
around the nominal Ks mass^ 1 ) This cut is particularly effective on Ks — > tt + tt~ events, 
which peak at M ee <~ 409 MeV. After it the background is dominated by Ks — ► n + ir~ 
with at least one pion wrongly reconstructed, and by <f> — > 7r + 7r~7r° events. The Ks events 
are strongly reduced by cutting on the track momentum in the Ks rest frame, which is 
expected to be ~ 206 MeV for Ks — ► 7r + 7r~ decays. To reject <j) — * 7r + 7r~7r° events we use 
instead |p m i SS | = \p# — pk s ~ Pk l \, where Ks and Kl momenta are evaluated from the 
charged tracks and from the K-crash tag, respectively. The value of \p m iss I peaks at zero 
for the signal, within the experimental resolution of few MeV. It spreads towards higher 
values for </> — > tt + tt^tt° events. The residual background, both from Ks — ► ir + n~ and 
4> — > 7T + 7r _ 7r decays, is rejected identifying the two electrons by ToF, and by using the 
properties of the associated calorimetric cluster. The reliability of the MC background 
simulation is checked after each step of the selection on the invariant mass sidebands. At 
the end of the analysis chain, we don't count any event in the signal box; the background 
estimate is also compatible with zero. The upper limit on the number of signal events is 
therefore N ee = 2.3 at 90% CL The signal selection efficiency, after tagging, is ~ 47%. 
Such performances in terms of exceptional background rejection (> 10 8 ) with an accept- 
able signal efficiency, have been achieved largely thanks to the very good momentum 
resolution of our DC. To obtain an upper limit for the BR(Ks — ► e + e~), we normalize 
N ee to the Ks — » 7r + 7r~ decays observed in the same data sample. For this purpose, we 
used the same selection criteria as for the measurement of BR(Ks — » 7r + 7r~), described 
in sect. 5'4.1. Finally, we obtain as a preliminary result [43] 



which represents a factor of ~ 15 improvement with respect to the best previous limit. 



( x ) With the invariant mass cut given above, we actually measure Ks — » e + e (7), with E* < 
20 MeV. 



(20) 




(21) 



BR(K S -> e+e-) < 9.3 x 10 



-9 



at 90% CL 
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5'5. Kl Decays. - A measurement of the absolute Kl BRs is a unique possibility of 
the ^-factory, where a pure sample of nearly monochromatic K^s can be selected by 
identification of a simultaneous K$ — > 7r + 7r~ decay. The absolute BRs can be deter- 
mined by counting the fraction of Kl's that decay into each channel, and correcting for 
acceptances, reconstruction efficiencies, and background. More specifically, we evaluate 
Kl BRs as: 



where Nf is the number of decays to / identified in the FV after background subtraction, 
A/tag is the number of tagged K^s, epy is the fraction of decays in the FV, e rcc is the 
reconstruction efficiency for channel /, and £tag/(etag) is the tag bias, i.e. the ratio of 
the Kg — > 7r + 7r~ reconstruction efficiencies for Kl —> f and independent of the Kl fate, 
respectively. We use a FV well within the DC, which gives an efficiency epv ~ 26%. 
Losses of Kl's from interactions in the beam pipe and chamber walls are taken into 
account in the evaluation of epv, as well as its dependence on the Kl lifetime (see 
sect. 5'2). 

5'5.1. Kl semileptonic and 3 pion decays. The four major decay modes, ir^e^v (K e s), 
7r ± /x =F ^ (K^), 7T + 7r~7r°, and 3tt° account for more than 99.5% of all decays, and are 
measured simultaneously from a sample of 13 millions tagged Kl decays. For the first 
three modes, two tracks are observed in the DC, whereas for the 37r° mode, only photons 
appear in the final state. The analysis of two-track and all-neutral-particle events is 
therefore different. 

Two-track events are assigned to the three channels of interest by use of a single 
variable: the smaller absolute value of the two possible values of A M7r = |p m i ss | — E misS} 
where p m i ss and E m i ss are the missing momentum and energy in the Kl decay, respec- 
tively, and are evaluated assuming the decay particles are a pion and a muon. Fig. 20, 
left, shows an example of a A Mir distribution. We obtain the numbers of K e3 , K^, and 
Tv + TT~ir° decays by fitting the A M7r distribution with the corresponding MC-predicted 
shapes. The signal extraction procedure is tested using PID variables from the calorime- 
ter. This is illustrated in fig. 20, right, where the A e7r spectrum is shown for events 
with identified electrons, together with the results of a fit with MC shapes. The ^3(7) 
radiative tail is clearly evident. The inclusion of radiative processes in the simulation is 
necessary to obtain a good fit, as well as to properly estimate the fully inclusive radia- 
tive rates. The reconstruction efficiency is ~ 54% for K e3 , <~ 52% for K^ 3 , and ~ 38% 
for 7T + 7r~7r°, as evaluated from MC. These efficiencies are then corrected to account for 
MC imperfections in reproducing the tracking efficiency The correction factors range 
between 0.99 and 1.03, depending on the decay channel and on the data taking period. 

The decay Kl — > 3ir° is easier to identify. Detection of > 3 photons originating at 
the same point along the Kl flight path is accomplished with the technique described 
in sect. 4'2. This selection guarantees very high reconstruction efficiency, ~ 99%, and 
small background contamination, ~ 1%. 

The tag bias is mostly due to a dependence of the calorimeter trigger efficiency on 
the Kl final state. Additional effects are introduced by the possible overlap between 
Ks and Kl charged products, which undermines Ks reconstruction performances for 
Kl decaying close to the IP. Both effects are reduced by additional requests to the tag 
selection, which give an overall tagging efficiency of <~ 9%. 
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Fig. 20. - Left: distribution of for a run set. Right: distribution of A e ^ for events with an 
identified electron, for the entire data set. 
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Our measured BR values for the four modes [26], corresponding to a reference value 
51.54 ns, are listed in table IV. The errors given for the absolute BRs do 



Parameter 


Value 




Correlation coefficients 


BR(°)(^ e3 ) 


0.4049(21) 


1 






0.2726(16) 


0.09 


1 


BR(°)(37r°) 


0.2018(24) 


0.07 


-0.03 1 


BR(°)(7r+7r-7r ) 


0.1276(15) 


0.49 


0.27 0.07 1 



Table IV. - KLOE measurements of principal Kl BRs assuming tl = 51.54 ns. 

not include the contribution from the Kr, lifetime value, which enters into the calculation 
of the FV efficiency. This dependence has been used in sect. 5 '2 to perform a combined 
fit to BRs and lifetime, including also the tj, measurement from the proper decay time 
distribution of Kt, — > 37r° decays. The results of this fit have been given in table II. 

5'5.2. Kl — > irlv^) form factors. The semileptonic kaon decay, see fig. 21, still pro- 

K x 




Fig. 21. - Amplitude for kaon semileptonic decay 



vides the best means for the measurement of \V US \, because only the vector part of the 
weak current contributes to the matrix element (n \ J a \ K }. In general, 



( 7T | J a \ K ) = /+ (t) (P + P ) a + f- (t) (P - P ) a , 
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where P and p are the kaon and pion four-momenta, respectively, and t = (P — p) 2 . The 
form factors (FF) / + and /_ appear because pions and kaons are not point-like particles, 
and also reflect both SU(2) and SU(3) breaking. Introducing the scalar FF fo(t), the 
matrix element above is written as 

(n(p)\u 7a s\K(P)) = f(0) ((P+ P ) a f + (t) + (P-p) a (/ (t)*p-/ + (i)^l 
(23) 

with Ak-k = fn 2 K ~ The / + and fo FFs must have the same value at t = 0. We 
have therefore factored out a term /(0). The functions f+(t) and fo(t) are both unity 
at t = 0. For vector transitions, the Ademollo-Gatto theorem [44] ensures that SU(3) 
breaking appears only to second order in m s — m Ut d- In fact, /+(0) differs from unity 
by only ~4%. The behavior of the reduced FFs f+(t) and fo(t) as a function of t can 
be measured from the decay spectra. At the level of accuracy reached by the present 
experiments, the parameterization used to extract the FF slopes is a relevant issue. 
If the FFs are expanded in powers of t up to t 2 as 

(24) f +At) = 1 + ^ J. + 1 . Kfi ^ ? 

four parameters (A' + , A" , A , and A ') need to be determined from the decay spectrum 
in order to be able to compute the phase space integral. However, this parametrization 
of the FFs is problematic, because the values for the As obtained from fits to the exper- 
imental decay spectrum are strongly correlated, as discussed in [45]. In particular, the 
correlation between A and A ' is —99.96%; that between and A" is —97.6%. It is 
therefore impossible to obtain meaningful results using this parameterization. 
Form factors can also by described by a pole form: 



(25) /+,„(*) = 



M 2 -t> 



which expands to 1 + t/M v s + (t/M v s ) 2 , neglecting powers of t greater than 2. It is 
not clear however what vector (V) and scalar (S) states should be used. 

Recent K e3 measurements [46, 47, 48] show that the vector FF is dominated by the 
closest vector (qq) state with one strange and one light quark (or Kir resonance, in an 
older language) . The pole-fit results are also consistent with predictions from a dispersive 
approach [49, 50]. We will therefore use a parametrization for the vector FF based on a 
dispersion relation twice subtracted at t = [49]: 



(26) /+(*) = exp 



Jr(A+ + #(*)) 



where H(t) is obtained using K — it scattering data. An approximation to eq. 26 is 



with p 2 and p% as given in table V. 
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Vector (eq. 27) 


Scalar (eq. 31) 


P2 X 10 4 


5.84 ± 0.93 


4.16 ±0.50 


P3 X 10 4 


0.30 ±0.02 


0.27 ±0.01 



Table V. - Constants appearing in dispersive representations of vector and scalar form factors. 



The pion spectrum in decay has also been measured recently [46, 51, 52]. As 
discussed in [45], there is no sensitivity to Ag. All authors have fitted their data using a 
linear scalar FF: 

(28) / ( t ) = i + Ao-L 

Because of the strong correlation between X' and Ag, use of the linear rather than the 
quadratic parameterization gives a value for Ao that is greater than X' by an amount 
equal to about 3.5 times the value of Ag. To clarify this situation, it is necessary to 
obtain a form for f (t) with at least t and t 2 terms but with only one parameter. 

The Callan-Trciman relation [53] fixes the value of scalar FF at t — Ak-k (the so- 
called Callan-Treiman point) to the ratio of the pseudoscalar decay constants fx/ fir- 
This relation is slightly modified by SU(2)-breaking corrections [54]: 

(29) /o(A ^ ) = !_l_ + AcT , 

where Act is of the order of 1CP 3 . A recent parametrization for the scalar FF [49] 
allows the constraint given by the Callan-Treiman relation to be exploited. It is a twice- 
subtracted representation of the FF at t = Ak-k and t = 0: 

(30) f (t) = cxp (-L- log(C - G(t)) 

such that C = /o(Aktt) and /o(0) = 1. G(t) is derived from Kit scattering data. As 
suggested in [49] , a good approximation to eq. 30 is 

t , A 2 ±p 2 / t \ 2 , Ao 3 ±3p 2 Ao±P3 / * 



(31) f ( t ) = l + \ — , . , 

with p 2 and p$ as given in table V. The Taylor expansion gives logC = XqAk-k /m 2 ± 
(0.0398 ±0.0041). Eq. 31 is quite similar to the result in [55]. 

To measure K e ^ form-factor slopes, we start from the same sample of Kt, decays 
to charged particles used to measure the main Kl BRs. We impose additional, loose 
kinematic cuts and make use of ToF information from the calorimeter clusters associated 
to the daughter tracks to obtain better particle identification. The result is a high- 
purity sample of 2 million Kl — > Kev decays. Within this sample, the identification 
of the electron and pion tracks is certain, so that the momentum transfer t can be 
safely evaluated from the momenta of the Kl and the daughter tracks. We obtain the 
vector form-factor slopes from binned log-likelihood fits to the t distribution. Using the 
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quadratic parameterization of eq. 24, we obtain 



and A" 



(25.5±1.5 stat ±1.0 syst )xl0- 



(1.4 ± 0.7 sta t ± 0.4 syst ) x 10 , where the total errors are correlated with p = 
—0.95. Using the pole parameterization of eq. 25, we obtain My = 870±6 s t a t±7 syst MeV. 
Evaluation of the phase space integral for Kl — » irev decays gives 0.15470±0. 00042 using 
the values of A' + and A" from the first fit and 0.15486 ± 0.00033 using the value of My 
from the second; these results differ by ~0.1%, while both fits give y 2 probabilities of 
<92%. The results we obtain using quadratic and pole fits are manifestly consistent. 

The measurement of the vector and scalar form-factor slopes using Kl — > n^iv de- 
cays is more complicated. First, there are two form factors to consider, and since all 
information about the structure of these form factors is contained in the distribution of 
pion energy (or equivalently, i), the correlations between form- factor slope parameters 
are very large. In particular, it is not possible to measure Aq for any conceivable level of 
experimental statistics [45]. Second, at KLOE energies, clean and efficient n/fi separa- 
tion is much more difficult to obtain than good ir/e separation. However, the form- factor 
slopes may be obtained from fits to the distribution of the neutrino energy E v , rather 
than to the distribution in t. E v is simply the missing momentum in the Kl — > -k^lv 
decay evaluated in the Kl rest frame, and requires no tt/ n assignment to calculate. A 
price is paid in statistical sensitivity: the E v distribution is related to the t distribution 
via an integration over the pion energy, and so the statistical errors on the form-factor 
slope parameters will be 2-3 times larger when the E v distribution, rather than the t 
distribution, is fit (assuming that the fit parameters are A' + , A" , and Ao). 
About 1.8 million decays were accepted. We first fit the data using equations 24 and 28 
for the vector and scalar FFs respectively. The result of this fit is shown in fig. 22. 
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Fig. 22. - Residuals of the fit (top plot) and E v distribution for data events superimposed on 
the fit result (bottom plot). 



We obtain 

A; = (22.3 ± 9.8 stat ± 3.7 syst ) x 10" 3 
(32) A'j = (4.8 ± 4.9 stat ± 1.6 syst ) x 10" 3 
A = (9.1±5.9 stat ±2.6 syst ) x 10~ 3 

with x 2 /dof = 19/29, and correlation coefficients as given in the matrix. Improved 
accuracy is obtained by combining the above results with those from our Kl € 3 analysis. 
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Wc then find: 

A' + = (25.6 ± 1.5 sta t ± 0.9 syst ) x 10~ 3 / 1 -0.95 0.29 \ 

(33) A" = (1.5 ± 0.7 sta t ± 0.4 syst ) x 10~ 3 1 -0.38 
A = (15.4 ± 1.8 stat ± 1.3 syst ) x 10" 3 V 1 / 

with x 2 /dof = 2.3/2 and the correlations given in the matrix on the right. Finally, the 
same combination of if e3 and results has been performed to take advantage of the 
recent paramcterizations of the FFs based on dispersive representations (eqs. 26 and 
30). We perform a fit to the values obtained for A+, A" , and A that makes use of the 
total error matrix as described above, and the constraints implied by equations 27 and 
31. Thus, the vector and scalar FFs are each described by a single parameter. Dropping 
the " ' " notations, we find: 

A+ = (25.7 ± 0.4 stat ± 0.4 syst ± 0.2 param ) x 10~ 3 

(34) , 
A = (14.0 ± 1.6 stat ± 1.3 syst ± 0.2 param ) x 10~ 3 

with % 2 /dof = 2.6/3 and a total correlation coefficient of —0.26. The uncertainties arising 
from the choice of parameterization for the vector and scalar FFs are given explicitly. 

The values of the phase space integrals for K(s decays are listed in table VI, for both 
determinations of the FF slopes (equations 33 and 34). We note that the use of the 



Parametrization 


I(K° e3 ) 








eq. 33 


0.15483(40) 


0.10271(52) 


0.15919(41) 


0.10568(54) 


eq. 34 


0.15477(35) 


0.10262(47) 


0.15913(36) 


0.10559(48) 



Table VI. - Phase space integrals for decays, obtained using FF slopes as determined with 
quadratic/linear (first row) and dispersive (second row) parameterizations of the FFs. 

dispersive parameterization changes the value of the phase space integrals by at most 
~ 0.09% with respect to what obtained using quadratic and linear parametrizations for 
vector and scalar FFs, respectively. This reflects into a systematic uncertainty from FF 
modelling of less than one per mil on the determination of /+ (0) V us , which is evaluated 
using the phase space factors from the FF dispersive parameterization. 

Finally, from the Callan-Treiman relation, eq. 29, we compute /+(0) = 0.967 ± 0.025 
using fx/fn = 1-189 ± 0.007 from a recent lattice calculation [57], and Act = (—3.5 ± 
8.0) x 10~ 3 . Our value for /+(0), although with a rather large error, is in agreement with 
the present best lattice determination, /+(0) = 0.9644 ± 0.0049 [58]. 

5'5.3. Radiative decay. Two different processes contribute to photon emission 

in kaon decays: inner bremsstrahlung (IB) and direct emission (DE). DE is radiation 
from intermediate hadronic states, and is sensitive to hadron structure. The relevant 
kinematic variables for the study of radiation in K13 decays are E* the energy of the 
radiated photon, and 9*, its angle with respect to the lepton momentum in the kaon 
rest frame. The IB amplitude diverges for E* — > 0. For K e ^, for which m e w 0, the 
IB spectrum in 8* is peaked near zero as well. The IB and DE amplitudes interfere. 
The contribution to the width from IB-DE interference is 1% or less of the purely IB 
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contribution; the purely DE contribution is negligible. To disentangle the IB and DE 
components, we measure the double differential spectrum d 2 T/dE*d9* and compare the 
result with expectations from Monte Carlo generators. In the ChPT treatment of [59], 
the photon spectrum is approximated by 

(35) + 

The DE contributions are summarized in the function f(E*), which represents the de- 
viation from a pure IB spectrum. All information on the strength of the DE terms is 
contained in the parameter (X). Theoretical predictions on this parameter suffer by large 
uncertainties, due to the poor knowledge of ChPT low energy constants. In contrast, a 
fit to the E*-9* spectrum allows us to measure for the first time a value for (X). From 
the fit, we extract also 

T(K eZl ; E; > 30 MeV , 9; > 20°) 

(36) R = 



K. 



•3(7) 



The value of this ratio has been computed at 0(p 6 ) in ChPT, leading to the prediction 
R = (0.963 ± 0.006(X) ± 0.010) x 10~ 2 . 

The experimental distribution in (£"*, 9*) has been fit using the sum of four indepen- 
dently normalized MC distributions: the distributions for K e ^ events from IB satisfying 
(not satisfying) the kinematic cuts E* > 30 MeV and 9* > 20°, the distribution cor- 
responding to the function f(E*) 1 and the physical background from Kl — ► ir + ii~ and 
Kl — * n^v events. Using a sample of 9 000 events, and 3.5 million -^3(7) f° r 

normalization, we find [60]: 

R = (0.924 ± 0.023 stat ± 0.016 syst ) x 10~ 2 
(37) {X) = -2.3±1.3 stat ±1.4 syst . 

The dependence of R on (X), predicted by ChPT, can be used to further constrain the 
possible values of R and (A") from our measurement. The constraint is applied via a fit, 
which gives R = (0.944 ± 0.014) x 10~ 2 and (X) = -2.8 ± 1.8, with x 2 /ndf = 0.64/1 
(P = 42%). The resulting value of R may be compared with the value quoted in [59], 
calculated for (X) = -1.2 ± 0.4. 

5"5.4. Kl — > n + n~ . CP violation was discovered in 1964 through the observation of 
the decay Kl^tt + it^ [61]. The value of BR(Al — > n + n~) is known today with high 
accuracy from the results of many experiments [21]. In the SM, CP violation is naturally 
accommodated by a phase in the quark mixing matrix [62, 63]. BR(ii£^7r + 7r~), together 
with the well known values of BR(iis^7r + 7r~), rg, and tl, determines the modulus of 
the amplitude ratio 



(38) |r?+_| = ^T(K L -> 7r+7r-)/r(^ s -> tt+tt"), 

which is parameterized as 77+ _ = e + e'. Here e quantifies the Kl CP impurity and e' 
is due to a direct CP-violating term. Since we know that e' ~ 10 _3 e [21], it follows 

that a measurement of |^_| | can be directly compared with the SM prediction for e. We 

obtain a precise determination of BR(ifL — ► tt + tt~) by measuring the ratio BR(Kl — > 
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-k + it~)/BR(Kl — > n^fi^v). This approach is particularly convenient, since the values 
of the tagging efficiencies for the Kl^tt + 'k~ and Kl — » it\iv decays are very similar, 
and the related systematic uncertainties do cancel in the ratio. For this analysis, we 
used the same data sample as for Kl BRs determination (sect. 5'5.1), and our value of 
BR(Kl — > w ± fi^u) has been used to extract BR(Kl — > 7r + 7r~)( 2 ). 

Given a Kg — ► 7r + 7r~ tagging decay, we select events with two tracks with opposite 
charge belonging to the same vertex along the Kl line of flight. A fiducial volume 
in the DC is chosen to cou nt Kl decays. The best variable to select Kl — > n + Tr~ 
decays is y 1 E 2 niss + \p m iss\ 2 , where the missing 4-momentum is evaluated using the tag 
information for Kl, and DC for Kl decay particles, which are assumed to be pions. We 
count 7r + 7r~ events by fitting the \J E 2 miss + \p m i SS \ 2 distribution, and using MC for the 
signal and the background shapes. The result of this fit is shown for a fraction of the 
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Fig. 23. - Distribution of \/ Ef nias + \p m iss\ 2 , with fit to MC distributions for different decay 
channels. Left, full fit range; right, signal region. 



whole data sample in fig. 23; the signal region is expanded in the same figure, right. 
The systematic uncertainty in the result is dominated by distortions introduced by MC 
in the simulation of signal and background shapes, which are particularly sensitive to 
momentum resolution. We obtain [64]: 

(39) BR(K L -» tt+O = (1.964 ± 0.012 stat ± 0.017 syst ) x 1(T 3 

This branching ratio measurement is fully inclusive of final-state radiation, and includes 
both the inner bremsstrahlung and the CP-conserving direct emission (DE) components. 
Using our measurements of BK(Ks — > tt + tt^) and tl, the value of rg from the PDG[21], 
and subtracting the contribution of the CP-conserving DE process [65] from the inclusive 
measurement of BR(iCL -» Tr+Tr"), we obtain |r? + _| = (2.219 ± 0.013) x 10~ 3 . Finally, 
using the world average of Re(e'/e) = (1.67 ± 0.26) x 10~ 3 and assuming arge' = arge, 
we obtain |e| = (2.216 ± 0.013) x 10~ 3 . The value of |e| can also be predicted from the 
measurement of the CP-conserving observables AM(Bd), AM(B S ), V u b, and V c b- This 
is particularly interesting to test the mechanism of the CP violation in the SM. For 



( 2 ) Statistical correlations have been taken into account in the determination of the measure- 
ment error. 
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this purpose, we use the prediction |e| = (2.875 ± 0.455) x 10~ 3 obtained in [66]. No 
significant deviation from the SM prediction is observed. Notice however that, due to 
the large uncertainties on the computation of the hadronic matrix element corresponding 
to the K-K mixing, the theoretical error is much larger than the experimental one. 

5'5.5. Kl—* 77- The Kl — ► 77 decay provides interesting tests of ChPT. This is 
because in the SU(3) limit it takes contribution from 0(p 6 ) terms only [67]. Moreover, 
a precise measurement of Kl — * 77 is needed to compute the absorptive part of the 
Kl — ► decay rate. This could be used to constrain the dispersive part of the 

same decay, which in turn is related to the parameter Vtd of the CKM matrix. From 
<~ 300 pb _1 of integrated luminosity, we measure the ratio Y(Kl — * 1i)/^{Kl — * 37r°). 
Using Kg — * 7r + 7r _ as a tagging decay, we select the signal by requiring two energetic 
photons belonging to a single vertex along the Kl line of flight. A good rejection power 
over the Kl — ► 37r° background is achieved by profiting of the two-body decay kinematics. 
After the selection, we count the signal by fitting the two-photon invariant mass, M 77 , 
with a combination of MC distribution for signal and background (fig. 24). From the 




Fig. 24. - Distribution of M 77 after the selection cuts; dots are data, histograms are MC shapes 
for signal (M 77 ~ tuk) and background (M 77 < ran)- 



previous fit, we count 22 185 ± 170 signal events, with efficiency ~ 81%. Finally, we 
obtain [68] 

(40) t&IZw) = (2 ' 79 ± °- 02stat ± °- 02s y^ x 10_3 - 

Using our result for BR(K L -> 37r°) decay, a decay width T(K L -» 77) = (7.2 ± 0.1) x 
10~ 12 eV is evaluated. This value is in reasonable agreement with 0(p 6 ) ChPT predic- 
tions, provided the value of the pseudoscalar mixing angle is close to our measurement 
of dp = ip P - arctan V2 = (-13.3 ± 0.3 stat ± 0.7 sys ± 0.6 th )° (section 9 - 2.2). 

5'6. if decays. - The charged kaons, K^, decay mostly into [i^v (K^a) and ^ir 
(K n 2), about 60% and 20% of all decays respectively, and the semileptonic modes account 



44 



F. BOSSI, E. DE LUCIA, J. LEE-FRANZINI, S. MISCETTI, M. PALUTAN ETC. 



for about 8%. In the following the measurement of the main branching ratios of K ± done 
at KLOE will be reported. These are absolute branching ratios determined tagging with 
two-body decays, easily identified as peaks in the p* distribution (sect. 4'3). A residual 
dependency of the tagging criteria on the decay mode of the tagged kaon is still present 
and it is accounted for in the final branching ratio evaluation. The few per mil accuracy 
on BR(i^2) and BKiK^) has been obtained tagging with K~ and using K + for signal 
search, neglecting corrections to the BRs from nuclear interactions (NI) of the kaon 
(<tni{K + ) <~ <tni(K~)/10 2 ). All the measurements are inclusive of final-state radiation. 
The results on BR(iQ^) and K ± lifetime t± (eq. 8) have been used to evaluate V us 
(sect. 6'2.1) and, together with the results from neutral kaons, to test lepton universality 
(sect. 6'2.2). From the BR(K^ 2 ) and t± values we have measured V us /V u d (sect. 6'3) 
and combining this result with the values of V us and V u d we have tested CKM unitarity 
(sect. 6'5). Leptonic two-body decays K+ 2 and K+ 2 allow us to put bounds on New 
Physics expected in Super-Symmetric extensions of the SM (sees. 5'6.5 and 6'4). 

5'6.1. K + — ► fi + v(-f). The KLOE measurement of this absolute BR [69] is based on the 
use of K~ — > [T v decays for event tagging. This choice minimizes possible interference in 
track reconstruction and cluster association of the tagging K~ and the tagged K + decays. 
The large number of K^ 2 decays allows for a statistical precision of ~0.1%, while setting 
aside a generous sample for systematic studies. For all the tagged events, the search for a 
positive kaon moving outwards from the IP in the DC with momentum 70 < px± < 130 
Mev/c is performed. Then kaon decay vertices in the DC volume, 40 < Rv < 150 cm, 
are selected. The number of signal counts is extracted from the distribution of p* , the 
momentum of the charged decay particle evaluated in the kaon reference frame and in 
the pion mass hypothesis, between p* min — 225 MeV/c and p* Max — 400 MeV/c. This 
distribution is shown in fig. 25 left. The spectrum shows the contamination, for a total 
amount of 2%, from K^ 2 and K + — > ir a l + u events. All the background sources have 
one neutral pion in the final state and therefore can be identified by the observation of 
two photon's clusters of equal time and compatible with the ir° mass. This selection 
provides the sample to evaluate on data the p* distribution of the background. Also 
the distribution for the signal events can be obtained from data control sample selected 
using EMC information only. This distribution is used together with the shape of the 
background sources to fit the overall spectrum and to perform background subtraction, 
fig. 25, center. Fig. 25 right shows the spectrum after background subtraction. The 
branching ratio is then obtained dividing the signal counts by the number of tagged 
events and correcting for the reconstruction and selection efficiency. The efficiency has 
been determined directly on data using a control sample of events selected exploiting 
their signature in the EMC. The control sample is constituted by events with K~ 2 ,^ 

providing the tag and ^2(7) selected using EMC information only. This critcrium is 
largely independent from the selection procedure based on DC information that have 
been used for obtaining the signal sample. 

In a sample of four million tagged events, KLOE finds ~865,000 signal events giving: 

(41) BR(K + -» fi+iy (7)) = 0.6366 ± 0.0009 sta t ± 0.0015 syst . 

This measurement is fully inclusive of final-state radiation and has a 0.27% uncertainty. 

5'6.2. K + — > 7r + 7r n (7). The measurement of the absolute BR of this decay, inclusive 
of radiation, is performed tagging with both K~ 2 and K~ 2 decays, providing a pure K + 
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Fig. 25. - Left. Monte Carlo spectra of p* for K + decays, showing contributions from various 
channels. Center. Distribution of p* . Right. Distribution of p* after background subtraction. 
The shaded area is used to count K + — > fi + v events. 



beam for signal search. The selection of K^ 2 (l) decays uses DC information only, with 
K + and its decay vertex selected as done in the K + — ► analysis. Loose cuts 

on p* and on the momentum difference between the kaon and the charged secondary 
track are applied to reject K ^ 3ir decays and if split tracks. The signal count is 
then extracted from the fit of the p* distribution in the window starting from p* ut = 180 
MeV/c (see fig. 26). This spectrum exhibits two peaks, the first at about 236 MeV/c 
from K~X 2 decays and the second at about 205 MeV/c from K^ 2 decays. Lower p* values 
are due to three body decays. The momenta of the charged secondaries produced in the 
kaon decay have been evaluated in the kaon rest frame using the pion mass hypothesis. 
Therefore the K^y, peak appears to be symmetric while the one is asymmetric due 
to the incorrect mass hypothesis used. 

The fit to the p* distribution is done using the following three contributions: K^, 

350000 
*jj 300000 

c 

> 250000 
200000 
150000 
100000 
50000 


1: 

Fig. 26. - Fit of the p* distribution. Left: black dots are Data and grey histogram is the fit 
output. Right: the three contributions used to fit the Data are shown: K^, K^2 and three-body 
decays. 
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and three-body decays. For and K n2 components we use shapes obtained from 
data control samples selected using EMC information only. For the contribution from 
three-body decays we use the MC distribution. Fig. 26 left shows the result of the fit 
of the p* distribution performed on the if~ 2 -tagged data sample. Fig. 26 right shows 
the three contributions: K v2 , K^ 2 and three-body decays. Using a total number of 
12,113,686 if~ 2 -tagged events we obtain 818, 347± 1, 912 signal counts. From the sample 
of 9,352,915 if ~ 2 -tagged events we get 621, 612± 1, 678 signal counts. The reconstruction 
and selection efficiency has been evaluated on data from a control sample selected using 
EMC information only, to avoid correlation with the DC driven sample selection. Once 
a tagging K~ 2 decay has been identified, the control sample selection is given by K + 
decays with a 7r° in the final state identified via the reconstruction of n° — > 77 decays. 
Corrections to the efficiency accounting for possible distortions induced by the control 
sample selection have been evaluated using MC simulation. The final efficiency evaluation 
is strongly related to the charged kaon lifetime r± via the geometrical acceptance. The 
BR depends on t± as: BR/BR/ ' =1-0.0395 ns _1 (t±-t±^) with reference value t±^ = 
12.385 ± 0.024 [21] used in the final BR evaluation. 

The weighted average, accounting for correlations, of the absolute branching ratios 
obtained using events tagged by K~ 2 and K~ 2 decays is a measurement with 4.6 per mil 
accuracy [70]: 

(42) BR(K+ -> 7r+7r°(7)) = 0.2065 ± 0.0005 sta t ± 0.0008 syst . 

The value reported is shifted by -1.3% (^2cr) with respect to the PDG06 fit value [21] 
and has a 20% improvement in the fractional accuracy. 

5'6.3. K ± semileptonic decays. The values of BR(iT ± — > 7r°e ± ;/(7)) and BR(.ff ± — > 
7r°/x ± f (7)) are each determined from four independent measurements: K + and K~ de- 
cays tagged by K — + [iv and K — > tttt . In the analyzed data set about 60 million tag 
decays were identified. The signal selection asks for a decay vertex in the DC volume. 
After removal of K n2 decays from the signal sample exploiting kinematics, the ir° is re- 
constructed from the two 7s which provide the kaon decay time ft e cav . This sample is 
composed mainly of semileptonic decays, residual two-body and K — > tx^tt^tx^ decays. 
To reject K n2 events with an early tt^ 1 — > \x^v decay, we cut on the lepton momentum 
evaluated in the center of mass of the ^ using the muon mass hypothesis. Then to isolate 
K e s and decays, the lepton is identified using a time of flight technique. Requiring 
the charged decay track to point to an energy deposit in EMC, the kaon decay time is 
given by tf^ y = ti ept - \Jpj ep t + m ie P t- Measuring the time of the cluster associ- 
ated to the charged decay particle ti ept , its momentum pi ept and track length L; ep t , we 
can evaluate mf ept by imposing t^ ay = tf^ v . The mf ept distribution is shown in fig. 27: 
the K e z component shows as a narrow peak around zero while the component is the 
peak around the value. The signal count is extracted from a constrained likelihood 
fit to this distribution, using a linear combination of K e ^ and shapes and of back- 
ground sources all taken from MC simulation. The fit result superimposed to the data 
distribution is also shown in fig. 27. The reconstruction and selection efficiency has been 
measured using MC simulation and corrected for relevant data/MC differences. Control 
samples have been selected to measure tracking and calorimeter clustering efficiency as 
a function of a suitable set of variables to perform data/MC comparisons. 
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Fig. 27. - Distribution of mf ept , for events selected as decays. Also shown is a fit to MC 
distributions for signal and background. 



The BRs depend on the value used for r± as 

(43) BR(X ± -> /)/BR (0) (X ± /) = 1 - 0.0364 ns" 1 (t± - t^A . 

Using the current world average value of the lifetime t±^ and averaging over the 
available samples, we get[71]: 

BR (0) (Af±) = (0.04965 ± 0.00038 sta t ± 0.00037 syst ) 

(44) BR (0) [K±z) = (0.03233 ± 0.00029 sta t ± 0.00026 syst ). 

These results, completely inclusive of final-state radiation, have a fractional accuracy 
of 1.1% for the K e % and of 1.2% for the decays. The total errors are correlated 
with 0.627 and do not include any contribution from the uncertainty on t±. Inserting in 
eq. 43 the KLOE result for t± from eq. 8, we obtain the values of fC ± semileptonic BRs 
reported in table VII, which have been used for our evaluation of V us . 

5'6.4. K —f 7r ± 7r°7r°(r'). This measurement is based on a sample of ~ 3.3 x 10 7 
tagged K ± [72], using K^ 2 and K^ 2 decays. The signal selection requires a decay vertex 
in the DC volume and at least four energy deposits in EMC, on-time with respect to 
the decay vertex, and p* < 135 MeV. Finally we require the sum of the energies of the 
most energetic photons to be J^^i < 450 MeV. The background sources are K^ 2 ,Kf 3 



Parameter 


Value 




Correlation coefficients 




BR(A- e3 ) 


0.04972(53) 


1 






br(a- m3 ) 


0.03237(39) 


+0.63 


1 






12.347(30) ns 


-0.10 


-0.09 


1 



Table VII. - KLOE measurements of K semileptonic decays and lifetime. 
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and radiative decays like K ± — » 7r ± 7r°7, in which a spurious cluster in EMC is paired 
with the cluster of the radiated photon mimicking the second 7r° . The total background 
fraction has been estimated from MC to be less than 1%. The overall efficiency is given 
by the product of the efficiency to reconstruct the kaon track, the track of the charged 
secondary and the decay vertex, times the efficiency of having at least 4 clusters fulfilling 
the energy and timing requirements. Each component has been measured separately for 
each tag and charge. Data control samples from EMC selection have been used for the 
evaluation of efficiencies involving DC variables and viceversa. Corrections evaluated 
using MC simulation account for possible contamination of the control samples. 

We get 52253 and 30798 r' decays from 1.9925xl0 7 events tagged by K± 2 and 
1.2753xl0 7 tagged by K^ 2 decays respectively. The final branching is: 

(45) BR(if ± -» 7r ± 7r°7r°) = (1.763 ± 0.013 stat ± 0.022 syst )% 

5'6.5. Rk =r(_ft' ± — > e ± v e )/r(K ± — ► /i ± ^ M ). For this measurement we decided to 
perform a "direct search" for — ► e ± ^ e and — ► n^v^ decays without tagging, 
to keep the statistical uncertainty on the number of — > e ± i/ e counts below 1%. 
The presence of a one-prong decay vertex in the DC volume with a secondary charged 
decay track with relatively high momentum (180-^270 MeV) is required. To distinguish 
— ► e ± ^ e and K ± — ► /i ± ^ M decays from the background we use the lepton mass 
squared obtained from the momenta of the kaon and the charged decay particle, 

assuming zero neutrino mass. This provides a clean identification of the K — > \i 
sample, while to improve background rejection in selecting — > e ± z/ e decays we need 
to use the information from the EMC. A particle identification based on the asymmetry 
of the energy deposits between the EMC planes, on the spread -Er.m.s. of these energy 
deposits and on the position of the plane with the maximum energy release is used. 

The number of signal events is extracted from a likelihood fit to the two-dimensional 
-Er.m.s. vs ^iop distribution, using MC shapes for signal and background. The shape 
for K — > ez/(7) have been evaluated considering both the Inner Bremsstrahlung (IB) and 
Direct Emission (DE) contributions. For signal count only events with photon energy in 
the kaon rest frame E* < 20 MeV have been considered. With this choice the contribution 
from DE is negligible and we measure the IB contribution, the only one entering in the SM 
prediction of Rk- The number of signal events obtained from the fit is Nk e 2 = 8090 ± 156 
and the projection of the fit result on the axes is compared to data in fig. 28 left. 
The number of K ± — > /x ± i/ A , events in the same data set is extracted from a similar fit to 
the Mjgp distribution, where no PID cuts are applied. The fraction of background events 
under the muon peak is estimated from MC to be less than one per mil. The number 
of K ± — ► [i^v^ events in the same data set is 499 251 584±35403. Using the number of 
observed — > e ± f e and K ± — ► [i^v^ events, we get the preliminary result [73]: 

(46) R K = (2.55 ± 0.05 ± 0.05) x 10~ 5 . 

This value is compatible within errors with the very precise SM prediction Rk = (2.472± 
0.001) x 10~ 5 . Recently it has been pointed out that, in a SUSY framework, sizeable 
violations of lepton universality can be expected in K\ 2 decays [74] from the couplings 
to the charged Higgs boson H + . With this Yukawa coupling, the dominant non-SM 
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Fig. 28. - Left: Distribution of the squared lepton mass of the charged secondary track M lcp : 
filled dots are data, open dots the result from the maximum-likelihood fit. Right: Exclusion 
limits at 95% CL on tan/3 and the charged Higgs mass m H ± from Rk for different values of 
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contribution to Rk modifies this ratio to: 
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With the lepton flavor violating term A13 of the order of 1CP 4 — 10~ 3 , as expected from 
neutrino mixing, and moderately large values of tan (3 and m# ± , SUSY contributions may 
therefore enhance Rk up to a few percent. Fig. 28 right shows the strong constraints on 
tan/3 and given by the KLOE Rk result. For a moderate value of A13 w 5 x 10~ 4 , 

the region tan/3 > 50 is excluded for charged Higgs masses up to 1000 GeV/c 2 at 95% 
CL. 

5'6.6. Fit to K BRs.. We fit the six largest BRs and the lifetime r using our 
measurements of r± (eq. 8), BR(K+ 2 ) (eq. 42), BR(K+ 2 ) (eq. 41), BR(K±) (eq. 44) and 
BR(_ft' ± — > 7r ± 7r°7r°) (eq. 45), with their dependence on r±, together with BR(_ft' ± — > 
^-k+tt-) = 0.0550 ± 0.0010 from the PDG04 average( 3 ) [75], with the sum of the BRs 
constrained to unity. The fit results [70], with x 2 /ndf = 0.59/1 (CL=44%), are shown in 
Table VIII. We can also evaluate BR(i4T ± — > 7r ± 7r + 7r _ ) by using our measurements of the 
above listed BRs and imposing the constraint ^BR(X ± -> /) = 1. With BR(A'+ 2 ) = 
0.63660 ± 0.00175, BR( J fl' ± -> ^^n") = 0.01763 ± 0.00025 and BR{K+ 2 ) = 0.20681 ± 
0.00094, BR{Kf 3 ) = 0.04972 ± 0.00053 and BR(^ 3 ) = 0.03237 ± 0.00039, evaluated 
at t± equal to our measured value 12.347±0.030 ns, we get BR(if ± — > 7r ± 7r + 7r _ ) = 
0.0568 ± 0.0022. This result is in agreement with the PDG04 average. 



6. V us an d lepton universality 



( 3 ) PDG '06 gives the result of their constrained fit but not the average of the data 
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Parameter 


Value 






Correlation coefficients 




BR(K^) 


0.6376(12) 










BR(^ 2 ) 


0.2071(9) 


+0.48 








BRi-K+n+n-) 


0.0553(9) 


-0.48 


+0.21 






BR(A~ e3 ) 


0.0498(5) 


+0.37 


-0.13 


+0.16 




BR(tf„ 3 ) 


0.0324(4) 


+0.34 


-0.12 


+0.15 +0.58 




BR(7r ± 7r t V°) 


0.01765(25) 


-0.11 


+0.05 


-0.05 +0.04 +0.04 




t± (ns) 


12.344(29) 


-0.15 


-0.21 


-0.07 -0.06 -0.05 


-0.015 



Table VIII. - Results of the fit to K± BRs. 



61. Introduction. - While much emphasis is placed on the search for New Physics, 
we still lack precise information on the validity of certain aspects of the SM itself. One 
of the postulates of the SM is the form of the Lagrangian for the charged- current weak 
interactions, which can be written as 

(48) C cc = -°-W+(U L V CK m7 q D l + e Ll a v eL + HT^^l + f Ll a v TL ) + h. c. 

v2 

Two important properties are evident: there is only one coupling constant for leptons 
and quarks, and quarks are mixed by the Cabibbo-Kobayashi-Maskawa matrix [62, 63], 
Vqkm, which must be unitary. The 4-fermion Fermi coupling constant Gf is related to 
the gauge coupling g by Gf = .g 2 /(4\/2m^). 

Precise measurements of leptonic and semileptonic kaon decay rates can provide infor- 
mation about lepton universality. Combined with results from nuclear (3 decay and pion 
decays, such measurements also provide information about the unitarity of the mixing 
matrix. Ultimately they tell us whether quarks and leptons do indeed carry the same 
weak charge. The universality of electron and muon interactions can be tested by measur- 
ing the ratio Y(K — > -k^v)/T{K — ► irev). The partial rates V (K — ► irev) &ndT(K — ► tt/j,p) 
provide measurements of g* \V US \ , which, combined with g 4 \V u d\ 2 from nuclear (3 decay 
and the muon decay rate, test the unitarity condition iVud^ + lVus^ + lVubl 2 = 1. We recall 
that in 1983 it was already known that \V ub \ 2 < 4 x 1CT 5 [76] and today |V„ b | 2 <~1.5 x 10~ 5 
[21]. We will therefore ignore | V^t, | 2 in the following. The ratio T(K — > [iv)/T(-n — > iiv) 
provides an independent measurement of |V^ S | 2 / | | 2 . 

To perform these tests at a meaningful level of accuracy, radiative effects must be care- 
fully corrected for. Strong-interaction effects also introduce form factors, which must be 
calculated from first principles, or measured whenever possible. Usually, corrections for 
SU(2)- and SU(3)-breaking must also be applied. Recently, advances in lattice calcula- 
tions have begun to catch up with experimental progress. 

6'2. Via semileptonic kaon decays . - Assuming lepton universality the muon decay 
rate provides the value for the Fermi constant, Gf=1-16637(1) x 1CU 5 GeV~ 2 , [21]. The 
semileptonic decay rates, fully inclusive of radiation, are given by 

(49) Y(K mi) ) = Cl ^ K S ™ l^l 2 W°) 2 ^ (! + 5 k (2) + ^) 2 ' 
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In the above expression, the index K denotes K° — > ^ and K ± — » 7r° transitions, 
for which C|- = 1 and 1/2, respectively, is the appropriate kaon mass, and 

Sew is the universal short-distance electroweak correction [77] ( 4 ). Conventionally, 
/ + (0) = f+ v (0), and the mode dependence is encoded in the 5 terms: the long- 
distance electromagnetic (EM) corrections, which depend on the meson charges and 
lcpton masses, and the SU(2)-breaking corrections, which depend on the kaon charge 
[79] . Iki is the integral of the Dalitz plot density over the physical region and includes 
|/+,o(*)| 2 (sect. 5'5.2). Iki does not account for radiative effects, which are included in 

rEM 
°Kl ■ 

The experimental inputs to eq. 49 are the semileptonic decay rates (tables II, III, 
VII) and the phase space integrals (table VI), which have been discussed in previous 
sections. SU{2)— breaking and EM corrections which are used to extract /+(0) V us are 
summarized in table IX. The SU(2)— breaking correction is evaluated with ChPT to 
0(p 4 ), as described in [79]. The long distance EM corrections to the full inclusive decay 
rate are evaluated with ChPT to 0(e 2 p 2 ) [79], and using low-energy constants from [80]. 
The quoted results have been evaluated recently [81], and include for the first time the 
K^z channels for both neutral and charged kaons. 



Channel 




°Kt 


K° e3 





0.57(15)% 







0.80(15)% 


4 


2.36(22)% 


0.08(15)% 


K p3 


2.36(22)% 


0.05(15)% 



Table IX. - Summary of SU "(2) — breaking and EM corrections used to extract /+(0) K 



6'2.1. /+(0) Vus and V us . Using all of the experimental and theoretical inputs discussed 
above, the values of /+(0) V us have been obtained for Kl c3 , Kl u3 , Ks e 3, K^, and K^ 3 
decay modes, as shown in table X and in fig. 29. It is worth noting that the only external 
experimental input to this analysis is the K$ lifetime. All other experimental inputs are 
KLOE results. 

The five different determinations have been averaged, taking into account all known 
correlations. We find [82] 



(50) f+(0)V us = 0.2157 ±0.0006, 

with x 2 / ndf = 7.0/4 ( 13%). To evaluate the reliability of the S'L7(2)-breaking cor- 
rection, a comparison is made between separate averages of f+{0)V us for the neutral 
and the charged channels, which are 0.2159(6) and 0.2145(13). With correlations taken 
into account, these values agree within l.ler. Alternatively, an experimental estimate of 



( 4 ) Often in literature the numerical factor 192 is replaced by 768 for which the phase space is 
normalized to 1 for all final particles masses vanishing [78]. The phase space integrals listed in 
sect. 5'5.2 are evaluated accordingly to the chosen convention. 
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Channel 


f+(0)V ue 




Correlation coefficients 




K Le3 
K Se3 
K% 


0.2155(7) 
0.2167(9) 
0.2153(14) 
0.2152(13) 
0.2132(15) 


1 

0.28 
0.16 
0.07 
0.01 


1 

0.08 
0.01 
0.18 


1 

0.04 1 
0.01 0.67 


1 



Table X. - KLOE results for /+(0) V ua - 



S K is obtained by comparing the neutral result with the charged one evaluated with- 
out correcting for SU{2)— breaking. We get Sfxp 2 ^ = 1.67(62)%, which is in agreement 
with the value estimated from theory (table IX). 

The previous determination of f + (0)V us with 2.8 per mil fractional accuracy, can be 
used together with the /+(0) value from [83], /+(0) = 0.961(8), to extract V us = 
0.2245 ± 0.0020. However, lattice evaluations of /+(0) are rapidly improving in pre- 
cision. For example, the RBC and UKQCD Collaborations have recently obtained 
/+ (0) = 0.9644(49) from a lattice calculation with 2 + 1 flavors of dynamical domain-wall 
fcrmions [58]. Using their value of /+(0), our results give V us = 0.2237(13). 

6'2.2. Lepton universality. Comparison of the values of f+(0)V us for K e 3 and 
modes provides a test of lepton universality. Specifically, 



(51) 



(/+(0)K,)' 3 . 



CXp 



r M 3 43(1 + s Ke ) 



Me (f + (o)v us ) 2 



e3,exp 



r e3 1^(1 + 5 



where Ski is shorthand for 5 K U ^ + S K ^f. By comparison with eq. 49 and given Gp = 
g 2 1 rriyy) , r ^e is equal to the ratio g^/g 2 , with ge the coupling strength at the 
W — > Iv vertex. In the SM, r Me = 1. Averaging between charged and neutral modes, we 
find 



(52) 



1.000 ±0.008. 



K 



Le3 



L|i3 



K 



Se3 



K 



e3 



M3 



0.213 0.215 0.217 
Fig. 29. - KLOE results for /+(0) Vus. 
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This has to be compared with the sensitivity obtained in it — > iv decays, 
(vOtt = 1.0042(33) [84], and in r leptonic decays, (r^) T = 1.000(4) [21]. 

6'3. Via K — > ^iv decay. - High-precision lattice quantum chromodynamics (QCD) 
results have recently become available and are rapidly improving. The availability of 
precise values for the pion- and kaon-decay constants f n and fx allows use of a rela- 
tion between T(K^ 2 )/r(n^ 2 ) and |V^ S | 2 / |V^d| 2 , with the advantage that lattice-scale 
uncertainties and radiative corrections largely cancel out in the ratio [85]: 



'k) 



, N T(K u2M ) \V US \ 2 ft m K (l-ml/m 2 K , 

53 \ M (7) / = ttt^ 9 X 0.9930 ± 0.0035 , 

1 ' r(7r„ 2(7) ) |K d | 2 A 2 m K (l- ml/ml) 2 1 ' 

where the uncertainty in the numerical factor is dominantly from structure-dependent 
corrections and may be improved. Thus, it could very well be that the abundant decays 
of pions and kaons to [lv ultimately give the most accurate determination of the ratio of 
1 14s | to | V u d | ■ This ratio can be combined with direct measurements of \V u d\ to obtain 
\V US \. From our measurements of BR^if^) (eq. 41) and t± (cq. 8), and using r(7r /j2 ) 
from [21], we evaluate: 

(54) Vus/Vud x f K /U = 0.2766 ± 0.0009. 

Using the recent lattice determination of f K / ' f v from the HPQCD/UKQCD collabora- 
tion, f K /f n = 1.189(7) [57], we finally obtain V us /V ud = 0.2326 ± 0.0015. 

6'4. Bounds on New Physics from K — > jiv decay. - A particularly interesting observ- 
able is the ratio of the V us values obtained from helicity suppressed and helicity allowed 
modes: R m = \V ua (K e2 )/V us (K e3 ) x V ud (0 + -» + )/V ud (n f _ t2 )\- This ratio, which is 
equal to 1 in the SM, would be affected by the presence of scalar density or extra right- 
handed currents. A scalar current due to a charged Higgs H + exchange is expected to 
lower the value of Rg 2 3, which becomes [86]: 



(55) R m 



1-^11 



m 2 H+ \ m 2 K+ J 1 + 0.01 tan /3 



tan 2 /? 



with tan (3 the ratio of the two Higgs vacuum expectation values in the MSSM. In ad- 
dition, in this scenario both + — > + nuclear beta decays and are not affected, 
and the unitarity constraint for this modes can be applied. To evaluate Re 2 3, we fit 
our experimental data on and K? 3 decays, using as external inputs the most recent 
lattice determinations of /+(0) [58] and fx/ fit [57], the value of V ud from [87], and 
V 2 d + V MS (Ki3) 2 — 1 as a constraint. We obtain 

(56) R m = 1.008 ±0.008, 

which is one o above the SM prediction. This measurement can be used to set bounds 
on the charged Higgs mass and tan/3. Fig. 30 shows the region excluded at 95% CL in 
the mfj+ — tan j3 plane. 

The measurement of BR(B — > tv) [88] can also be used to set bounds on the m H + — 
tan P plane, which are shown in fig. 30. While the B — > tv can exclude quite an extensive 
region of this plane, there is an uncovered region corresponding to the change of sign of 
the correction. This region is fully covered by our result. 



54 



F. BOSSI, E. DE LUCIA, J. LEE-FRANZINI, S. MISCETTI, M. PALUTAN ETC. 




100 200 300 400 500 



Fig. 30. - Excluded region in the m H + — tan/3 plane by the measurement R123', the region 
excluded by B — > rv is also indicated. 
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Fig. 31. - KLOE results for |K, S | 2 and |V^»/V^d| 2 , together with |Vud| 2 from /3-decay measure- 
ments. The ellipse is the la contour from the fit. The unitarity constraint is illustrated by the 
dashed line. 



6'5. Test of CKM . - To test the unitarity of the quark mixing matrix, we combine all 
the information from our measurements on K^, K e ^, K^, together with superallowed 
+ — > + nuclear f3 decays (fig. 31). The best estimate of \V US \ 2 and \V u d\ 2 can be 
obtained from a fit to our results V us — 0.2237(13) and V us /V u d = 0.2326(15), together 
with V ud = 0.97418(26) [87]. The fit gives \V US \ 2 = 0.0506(4) and \V ud \ 2 = 0.9490(5), 
with x 2 /ndf = 2.34/1 (13%) and a correlation coefficient of 3%. The values obtained 
confirm the unitarity of the CKM quark mixing matrix as applied to the first row. We 
find 

(57) \V US \ 2 + \V ud \ 2 - 1 = -0.0004 ±0.0007 (~ 0.6a) 

In a more conventional form, the results of the fit are: 

, |Ks I = 0.2249 ±0.0010 

58 

v ' Kd= 0.97417 ±0.00026. 
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Imposing unitarity as a constraint, |K( S | 2 + IV^dj 2 = 1, does not improve the accuracy. 
One should also keep in mind that while lattice results for /+(0) and Jk/ fir appear to be 
converging and are quoted with small errors there is still a rather large spread between 
different calculations. If we were to use instead /+(0) = 0.961 ± 0.008 as computed in 
[83] and still preferred by many authors, we find \V US \ — 0.2258 ± 0.0012 which is less 
precise but satisfies more closely unitarity. 

7. Quantum interferometry 

7'0.1. Quantum coherence. As stated in sect. 4'1, the neutral kaon pair from <p — > 
K°K° are produced in a pure quantum state with J PC = 1 . We evolve the initial state 
in time, project to any two possible final states f± and /2, take the modulus squared, and 
integrate over all t\ and t<z for fixed At = ti—tz to obtain (for At > 0, with T = Tl+Ts)'- 



(59) 



If uh (&t) = 7^\(fi\Ks)(f2\Ks)\ 2 x 



\m\ 2 e~ rLAt + \m\ 2 e- rsAt 



- 2(1-C)k||%|e- 



-(r L +r s )Ai/2 



C0s(A77lAt + 02 — 4>l) 



The last term is due to interference between the decays to states /i and and £=0 in 
quantum mechanics. Fits to the At distribution provide measurements of the magnitudes 
and phases of the parameters r\i = (fa \Ki)/{fi \Ks), as well as of the K^-Ks mass 
difference Am and the decay widths Tl and Tg- 

Such fits also allow tests of fundamental properties of quantum mechanics. For exam- 
ple, the persistence of quantum-mechanical coherence can be tested by choosing f± = /a. 
In this case, because of the antisymmetry of the initial state and the symmetry of the 
final state, there should be no events with At = 0. Using ~ 300 pb" 1 of integrated lumi- 
nosity, we published an analysis of the At distribution for K$Kl — * 7r + 7r~7r + 7r~ events 
that establishes the feasibility of such tests. The At distribution is fit with a function of 
the form of eq. 59, including the experimental resolution and the peak from Kl — * Kg 
regeneration in the beam pipe. The results are shown in fig. 32. Observation of Q ^ 
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.300 
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regeneration 
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Fig. 32. - At distribution for 4> — > KsKl — * ir + %~ ir + tt~ events. 



would imply loss of quantum coherence. The value of £ depends on the basis, K°-K° or 



56 



F. BOSSI, E. DE LUCIA, J. LEE-FRANZINI, S. MISCETTI, M. PALUTAN ETC. 



Ks-Kl used in the analysis. Using the K°-K° basis we find [89] 

(60) C - (0.10 ± 0.21 stat ± 0.04 syst ) x 10- 5 

which means no violation of quantum mechanics. The statistical error has been already 
reduced by a factor of two analyzing <~ 1 fb _1 more of integrated luminosity. 

7'0.2. Quantum interference and quantum gravity. In the context of a hypothetical 
quantum gravity theory, CPT violation effects might occur in correlated neutral kaon 
states [90, 91], where the resulting loss of particle-antiparticle identity could induce a 
breakdown of the correlation of state imposed by Bose statistics. As a result, the initial 
state can be parametrized as: 

\i) = (\K°,p)\K°,-p) - \K°,p)\K°-P)) 

(61) +^(\K\p ) \K\-P) + \K\p)\K\-P)) , 

where w is a complex parameter describing a completely novel CPT violation phe- 
nomenon, not included in previous analyses. Its order of magnitude could be at most 

1 /9 

\uj\ <~ [(m| r /mpi anc k)/Ar] <~ 10~ 3 with Ar = T$ — Tl. Inserting the modified initial 
state in the expression of the KgKj^ — > 7r + 7r~7r + 7r~ decay intensity, and fitting the At 
distribution as in the case of the previous analysis, we obtained the first measurement of 
the complex parameter uo [89]: 

(62) Re(w) = (1.1±5.3 ± 0.9) x 10~ 4 Im(w) = (3.4±|;g ± 0.6) x 10~ 4 . 

This can be translated into an upper limit on the absolute value, u> < 2.1 x 10~ 3 at 
95% CL. Also in this case the statistical error has been reduced by a factor of two 
analyzing ~ 1 fb _1 more of integrated luminosity, thus reaching the interesting Planck's 
scale region. 

8. Test of CPT 

8T. CPT SM . - The three discrete symmetries of quantum mechanics, charge con- 
jugation (C), parity (P) and time reversal (T) are known to be violated in nature, both 
singly and in pairs. Only CPT appears to be an exact symmetry of nature. Exact 
CPT invariance holds in quantum field theory which assumes Lorentz invariance (flat 
space), locality and unitarity [92]. These assumptions could be violated at very large 
energy scales, where quantum gravity cannot be ignored [93]. Testing the validity of 
CPT invariance probes the most fundamental assumptions of our present understanding 
of particles and their interactions. The neutral kaon system offers unique possibilities for 
this study. 

Within the Wigner-Weisskopf approximation, the time evolution of the neutral kaon 
system is described by 



(63) 



i^(t) = mt) = (M- l -TMt), 
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where M and T are 2x2 time-independent Hermitian matrices and $f(t) is a two- 
component state vector in the K°-K° space. Denoting by rriij and the elements of 
M and T in the K°-K° basis, CPT invariance implies 



(64) toh=to 2 2 (orm Ko =m Io ) and 
The eigenstates of eq. 63 can be written as 
1 



rn = r 



22 



(° r r K„=% ) 



(65) K s ,l 



-ilm(mi 2 ) - |lm(ri 2 ) ± \ 



[(l + e s .L)K a ±(l-e s . L )K° 



m K a ~ m K 



( T Ko- T K a ) 



m L - m s + i(T s -T L )/2 



e±S, 



such that S — in the limit of exact CPT invariance. Unitarity allows us to express the 
four entries of T in terms of appropriate combination of kaon decay amplitudes Ai : 



(66) 



Y^A(f)Mf)\ 
f 



i,j = 1,2 = K°,K°, 



where the sum runs over all the accessible final states. Using this decomposition in eq. 65 
leads to the Bcll-Steinberger relation [94]: a link between Re(e), Im(<5) and the physical 
kaon decay amplitudes. In particular, without any expansion in the CPT-conserving 
parameters and neglecting only 0(e) corrections to the coefficient of the CPT- violating 
parameter (5, we find 



(67) 



Vr s -r L 



+ i tan 0sw 



Re(e) 



- ilm(5) 



A L (f)A* s (f), 



where As(f) and Al{J) are the physical K s and K L decay amplitudes, and (f> sw = 
arctan (2(m£ — ms)/(Ts — ^l))- The advantage of the neutral kaon system is that only 
the 7T7r(7), iririr and ir£v decay modes give significant contributions to the right-hand side 
of eq. 67. The solution to the unitarity relation in eq. 67 is: 



(68) 



Rc(e) 

i + kl 2 

Im(5) 



1 

N 



1 + k(1-26) (1 
(1 — k) tan^sw 



k) tan <psw 
-(l + «) 



S 4 Re(a,) 
S i Im(a i ) 



where the a parameters are related to the product of Ks and Kl decay amplitudes to 
ttit(j), tttttt and ttIi> final states, k = ts/tl, b = BR(Kl — » irlv), and 



(69) 



N = (1 + k) 2 + (1 - k) 2 tan 2 ^ S w - 26k(1 + k) . 



Recently we published improved results for the two parameters Re(e) and lm(5) us- 
ing eq. 67 and our measurements of neutral kaon decays [95]. Our analysis benefits in 
particular from three measurements: i) the branching ratio for Kl decays to 7r + 7r~[64] 
(sect. 5'5.4), which is relevant for the determination of Re(e); ii) the new upper limit on 
BR(ifs^7r 7r 7r°) [39] (sect. 5'4.4), which is necessary to improve the accuracy on Im(5); 
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and iii) the measurement of the K$ semileptonic charge asymmetry As [31] (sect. 5'4.2), 
which allows, for the first time, the complete determination of the direct contribution 
from semileptonic channels, without assuming unitarity. The experimental inputs to the 
determination of the a parameters of eq. 68 are all of the KLOE neutral kaon measure- 
ments (tables II and III), the values of ts, ttll —ms, 4>+-, 0oo from [21], and the tt + tt~j 
and 7r + 7r~7r° amplitudes from [96] and [97], respectively. Using all of these results we 
obtain the values reported in table XI; fig. 33 and fig. 34 show the 68% and 95% CL 
contours in the complex plane Re(a») — Im(a!j). 



Parameter Value 



a n+n - ((1.115 ± 0.015) + i(1.055 ± 0.015)) x 10~ 3 

0,0,1, ((0.489 ± 0.007) + i(0.468 ± 0.007)) x 10~ 3 
a w+n - lDE ((5±7) + t(6±7)) x 10~ 7 

ctV+^o ((0±2) +i(0±2)) x HT 6 

Ko^o] < 7 x 10" 6 at 95% CL 

a nlv ((0.3 ± 0.6) + i(-1.8 ± 1.8)) x 10~ 5 



TABLE XI. - Values of a parameters of eq. 68 for all relevant kaon decay channels. 
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Fig. 33. - Representation of a 7r + n - and a^o^o in the complex plane. In each case, the two 
ellipses represent the 68% and the 95% CL contours. 



Inserting all of the information in eq. 68 we finally obtain: 

(70) Re(e) = (159.6 ± 1.3) x 10~ 5 , Im(<5) = (0.4 ± 2.1) x 10~ 5 . 

The allowed region in the (Re(e), Im(<5)) plane at 68% CL and 95% CL is shown in fig. 35, 
left. Our results, eq. 70, improve by a factor of two the previous best determinations. 
The limits on Im(<5) and Re(e) can be used to constrain the mass and width difference 
between K° and K° via 

(71) 8 = i(mK ° ' m * 0) + 1 f K ° ' ^ cos^e^[l + 0(e)}. 

1 S — 1 L 

The allowed region in the AM = {m^o — rng-o), Ar = (Tk° — T^o) plane is shown in 
fig. 35, right. Since the total decay widths are dominated by long-distance dynamics, 
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Fig. 34. - Bounds of the a values for the three-body decays of Ks, Kl- Note that the same 
scale is used for all plots. 



in models where CPT invariance is a pure short-distance phenomenon, it is useful to 
consider the limit T K o = T^ - In this limit, neglecting CPT— violating effects in the 
decay amplitudes, we obtain the following bound on the neutral kaon mass difference: 

(72) -5.3 x 10~ 19 GeV < m K o - m R0 < 6.3 x 10~ 19 GeV at 95 % CL . 
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Fig. 35. - Left: allowed region at 68% and 95% CL in the Ree, Im<5 plane. Right: allowed region 
at 68% and 95% CL in the AM, AT plane. 



60 



F. BOSSI, E. DE LUCIA, J. LEE-FRANZINI, S. MISCETTI, M. PALUTAN ETC. 



8'2. CPT and Lorentz symmetry breaking. - CPT invariance holds for any realistic 
Lorentz-invariant quantum field theory. However a very general theoretical possibility for 
CPT violation is based on spontaneous breaking of Lorentz symmetry, as developed by 
Kostelecky [98] , which appears to be compatible with the basic tenets of quantum field 
theory and retains the property of gauge invariance and renormalizability (Standard 
Model Extensions, SME). For neutral kaons, CPT violation would manifest to lowest 
order only in the parameter 5, with dependence on the kaon 4-momentum: 

(73) 5 w i sin0svi/e^ sw 7K(Aao - (? K ■ Aa) /Am, 

where -fx and 0k are the kaon boost factor and velocity in the observer frame, and 
Aa^ are four CPT- and Lorentz-violating coefficients for the two valence quarks in the 
kaon. Following [98], the time dependence arising from the rotation of the earth can 
be explicitly displayed in eq. 73 by choosing a three-dimensional basis (X,Y,Z) in a 
non-rotating frame, with the Z axis along the earth's rotation axis, and a basis (x, y, z) 
for the rotating (laboratory) frame. The laboratory frame precesses around the earth's 
rotation axis Z at the sidereal frequency Q. Defining % the angle between Z and the 
positron beam direction, the CPT violating parameter 5 may then be expressed as: 

2tt 



s(\p\,e,t) = j-J^ * 5(p,t)dcj> 



= t Ik { Aa + (3 K Aa z cos 6 cos x 

Am 

(74) +/3if(Aay sin x cos # sin fit + Aax sinxcosflcosQt)} , 

The Aao parameter can be measured through the difference between K$ and Kl charge 
asymmetries, As — Al- Each asymmetry is integrated over the polar angle 9, thus 
averaging to zero any possible contribution from the terms proportional to cos# in eq. 74, 
giving the relation: 



(75) 



A s - A L ~ 



4Re {i sin (f>s\ye 



IK 



Am 



Aa ■ 



From our measured value of As [31] (sect. 5'4.2) and a preliminary evaluation of Al, we 
get As — Al = (—2 ± 10) x 10~ 3 , which gives as a preliminary evaluation of the Aa 
parameter [10]: 



(76) Aa = (0.4 ± 1.8) x 10~ 17 GcV . 

With the analysis of the full data sample, an accuracy a(Aa a ) ~ 7 x 10~ 18 GcV could 
be reached. 

The parameters Aax,Y,z can be measured using <fi — ► KsKl — * tt + -k~ ,n + n~ events, 
by fitting the decay intensity / (n + n~ (cos 6 > 0), n + n~ (cos 6 < 0); At), where the two 
identical final states are distinguished by their forward or backward emission. A prelim- 
inary analysis, based on ~ 1 fb _1 , yields [10]: 

Aa x = (-6.3 ±6.0) x 10~ 18 GeV 
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Aa Y = (2.8 ± 5.9) x 10" 18 GeV 
(77) Aa z = (2.4 ± 9.7) x 10~ 18 GeV , 

which represents the first determination of Aaz to date. 
9. Light meson spectroscopy 




BR 0.1% 

Fig. 36. - Graphical representation of <f> meson decaying to lower lying meson states. 

When running DA$NE at the <f> peak, the dominant decay channels of the <fi meson 
are K K pairs. However, as shown in fig. 36, the 4> meson decays with large rates also to 
other interesting hadronic final states. In this chapter, we mainly discuss the so called <f> 
radiative decays i.e. electromagnetic transitions (</> — > meson + 7) to other mesons, which 
result in a great profusion of lighter scalar and pseudoscalar mesons. The transition rates 
are strongly dependent on the wave function of the final-state meson. They also depend 
on its flavor content, because the is a nearly pure ss state and because there is no 
photon-gluon coupling. These radiative <j> decays are unique probes of meson properties 
and structure. This hadron source, coupled with KLOE 's great versatility, allows to 
perform many hadronic experiments, with great precision. Moreover, we can study with 
high statistics the properties of vector mesons. Due to space constraint, we will limit 
ourselves to mentioning only a handful of these measurements, just enough to give a 
taste of hadron physics with KLOE. 

91. rj meson. - The BR for the <p -> 777 decay is 1.3%. In 2.5 ib" 1 of KLOE data, 
there are ~ 100 million 77 mesons which, for most of the final states, are clearly identified 
by their recoil against a photon of E = 363 McV. The jj's decay predominantly into 
two photons or three pions. The large production rate and the clean 77 signature made 
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possible to search for rare or forbidden 77 decays even with reduced samples. A good 
example is provided by the decay 77 — > n + -K~ which, being the 77 an isoscalar, violates 
both P and CP. With the first 350 pb" 1 , KLOE had searched for evidence of this 
decay in the M M distribution of e + e~ — > 7r + 7r~7 events where the photon is emitted 
at large polar angles (9 > 45°). No peak is observed in the distribution of M n7T in the 
vicinity of m v . The corresponding limit is BK(rj — ► it + it~) < 1.3 x 10~ 5 at 90% CL 
[99], which improves the previous limit by a factor of 25. Similarly, the decay 77 — > 37 
violates C. KLOE has conducted a search for this decay using 410 pb _1 and has set the 
limit BR(rj — > 37) < 1.6 x 10 -5 at 90% CL, the most stringent result obtained to date 
[100]. In the following subsections, we will instead show precise measurement on the 77 
meson such as the mass and the description of the dynamics for the 77 — > 37r decays. We 
conclude by showing preliminary BR measurements for very rare but observed decays 
which arc still in progress. 

9T.1. 77 mass. KLOE has performed the most precise determination to date of the 
77 mass, on the basis of 17 millions of <fi — > 777, r\ — > 77 decays. The <f> — > 777 events 
are selected by requiring at least three energy deposits in the barrel with a polar angle 
6> 7 : 50° < # 7 < 130°, not associated to a charged track. A kinematic fit imposing 
energy-momentum conservation and time of flight of photons equal to the velocity of 
light is done for all 3 7's combination of N detected photons. The combination with 
the lowest \ 2 is chosen as a candidate event if x 2 < 35. The inputs of the fit are the 
energy, the position and the time of the calorimeter clusters, the mean position of the 
e + e~ interaction point, x^, and the total four-momentum, p^, of the colliding e + e~ pair. 
The mean values of and p^ are determined run by run using e + e~ — ► e + e~ events 
(almost 90000 events for each run, allowing a very precise determination of the relevant 
parameters). 

The events surviving the cuts are shown in the Dalitz plot, fig. 37-lcft, where three 
bands are clearly visible. The band at low m 2 7 is given by the <p ~> i"°7j 71-0 ~ > 77, 
while the other two bands are (p — ► 777, 77 — > 77 events. Using the cuts indicated by the 
black line shown in the Dalitz plot we select pure samples of 77, ir° — > 77 events. The 
resulting m 77 spectrum for the 77 (fig. 37-right) is well fitted with a single gaussian of 
a ~ 2.1 McV/c 2 , thus showing that the kinematic fit improves of a factor ~ 20 the mass 
resolution obtained by calorimetric reconstruction only. 

To determine the systematic error associated with the mass measurement, we have 
evaluated the uncertainties on all the quantities used in the kinematic fit and their effect 
on the fitted value. A sample of e+e~ — > 7r + 7r~7 events has been used to check the mean 
position of the interaction point and the alignment of the calorimeter with respect to 
the Drift Chamber. The absolute energy scale of the calorimeter and the linearity of the 
energy response was checked using both the e + e~ — > e + e~j and the 7r + 7r _ 7 events. A 
linearity of better than 2% is found and the absolute scale results set to better than 1%. 
Since the kinematic fit overconstrains the photon energy with the cluster positions, the 
systematic uncertainties related to the energy reconstruction are small and result in a 
77 mass error of 4 kcV for the scale and 4 keV for the linearity. For the same reason, a 
larger effect can be expected if the calorimeter does not have angular uniformity in </> or 9. 
Misalignment of single modules in the calorimeter has been checked resulting in a spread 
of about 10-15 keV on the 77 mass. While the chosen cut on \ 2 has negligible effect, the 
systematics due to the particular choice of the event selection cut on the Dalitz plot, 
shown in fig. 37-left, was determined to contribute 12 keV to the error. The measured 
value of the mass is, however, very sensitive to the center of mass energy of the 777 system. 
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Fig. 37. - Left: Dalitz plot distribution of the selected 37 events. The photons are sorted 
according their energies E\ < E2 < -E3. Right: Invariant mass distribution of the 7172 pair in 
the 77 region. The fitted function is the continuous line. 



Due to initial state radiation emission (ISR), the available center of mass energy is a bit 
lower than the one of the e + e~ beams, W, measured using e + e~ — * e + e~ events. A 
variation of 100 keV of the measured mass value is obtained from the MC simulation 
which includes ISR emission. Due to the large value of this correction, the rj mass has 
been determined as a function of W both in data and with the MC simulation+ISR. The 
resulting data-MC spread of 8 keV is assumed as systematic error. 

All the above mentioned studies on systematics have been also done for the tt° mass 
using the <fi — * 7r°7 events and for the ratio of the two masses r = m^jm^a. The mass 
scale is set by the knowledge of the absolute scale of W determined using the CMD2 
value as shown in the K mass measurement section. We obtain [25]: 

(78) m„o = (134.906 ± 0.012 stat ± 0.048 syst ) MeV 



(79) m„ = (547.874 ± 0.007 sta t ± 0.031 syst ) MeV 

the 7T° mass value is in agreement with the world average [21] within 1.4 a. 
As a check of this result, we use the measured ratio 

(80) r = = 4.0610 ± 0.0004 stat ± 0.0014 syst 

and the world average value of the 7T° mass, m r o = (134.9766 ± 0.0006) MeV to derive 
rn v = (548.14±0.05 s tat±0.19 sys t) MeV, consistent with the results quoted above although 
affected by a larger systematic error. This is due to the soft energy spectrum of the 
photons from ir° which results in a worse position resolution. 

While being the most accurate result today, our measurement of the rj mass is in good 
agreement with the recent determinations based on ij decays [101, 102]. Averaging those 
mass values and our result we obtain m, = 547.851 ± 0.025 MeV, a value different by ~ 
10 a from the average of the measurements done studying the production of the 77 meson 
at threshold in nuclear reactions [103]. 
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9'1.2. Dynamics of the decays rj — > 3tt. The decay 77 — > 37r violates isospin invariance. 
Electromagnetic contributions to the process are very small and the decay is induced 
dominantly by the strong interaction via the u, d mass difference. The 77 — > 3n decay 
is therefore an ideal laboratory for testing Chiral Perturbation Theory, ChPT. We have 
studied the dynamics of the decay rj — ► 7t + tt~?t using about 1.4 million of such events. 
From a fit to the Dalitz plot density distribution we made precise determinations of the 
parameters that characterize the decay amplitude. In addition, we have also obtained 
the most stringent tests on C-violation from left-right, quadrant and sextant integrated 
asymmetries. [104] 

A three body decay( 5 ) is fully described by two variables. We can choose two of the 
pion energies (E + ,E_,Eq) in the 77 rest frame, two of the three combinations of the 
two-pion masses squared (m 2 i __, m 2 _ 0l TOq + ) also called (s,t,u). Note that E + is linear 
in rn^_Q and so on, cyclically. The Dalitz variables, X, Y are linear combinations of the 
pion energies: 

E -mo 3 r 2 1 

where Q is the reaction "Q-value" . The decay amplitude is given in [105] as: 

(82) A ^ u) -^^- m ^^Mr 

where A 2 = (m 2 — fh 2 )/(m 2 — m 2 ) and m = l/2(m„ + ra d ) is the average u, d quark 
mass. — 92.4 MeV is the pion decay constant and Wl(s, t, u) the amplitude we would 
like to know. From eq. 82 it follows that the decay rate for rj — > ir + ir~ir° is proportional 
to A -4 . The transition rj — > 3tt would be therefore very sensitive to A if the amplitude 
9JI is known. At lowest order in ChPT: 

(83) m(s,t,u) = 3s ~ Am l . 

m 2 -ml 

and the decay width at leading order is T l0 (77 — > 7r + 7r~7r ) = 66 eV [105] to be compared 
with the measured width 295 ±16 eV[21]. A one-loop calculation within conventional 
ChPT improves considerably the prediction giving r nl ° (77 — ► tt + 7t~7t°) ~ 167 ± 50 cV 
but is still far from the experimental value. Higher order corrections help but do not 
yet bring agreement with measurements of both total rate and Dalitz plot slopes. Good 
agreement is found combining ChPT with a non perturbative coupled channels approach 
using the Bethe Salpeter equation [106]. Therefore a precision study of the 77 — > 3tt Dalitz 
plot, DP, is highly desirable. The amplitude squared is expanded around X = Y = in 
power of X and Y 

(84) \A(X,Y)\ 2 = l + aY + bY 2 + cX + dX 2 + eXY + fY 3 + ... 



( 5 ) Both rj and 7r are spinless, therefore there is no preferred direction. 
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The parameters (a, b, c, d, e, /, ... ) can be obtained from a fit to the observed DP density, 
see Fig. 38 and should be computed by the theory. Any odd power of X in A(X, Y) 
implies violation of charge conjugation. 

The event selection consists of a pre-selection, the requirement of one charged vertex 
inside a small region together with three photons and a total photon energy below 800 
MeV. A constrained kinematic fit imposing four-momentum conservation and correct 
time of arrival for the photons is then applied and only events with a % 2 probability 
greater than 1% are kept for further analysis. We finally required that the recoil photon 
energy lies between 320 and 400 MeV and we add two further constraints on the sum 
of charged pion energies and on the invariant mass of the two softest photons. All 
selection criteria and efficiencies were checked by data and Monte Carlo. The overall 
selection efficiency, taking into account all the Data-MC corrections is found to be e = 
(33.4 ± 0.2)%. The expected background, obtained from MC simulation is as low as 
0.3%. After the background subtraction the observed number of events in the Dalitz plot 
is N bs — 1-337 ± 0.001 million events. The results, from fitting the matrix element of 




Fig. 38. - Dalitz-plot distribution for the whole data sample. The plot contains 1.34 millions of 
events in 256 bins. 

Dalitz-plot as shown in eq. 84, including the statistical uncertainties coming from the fit 
and the estimate of systematics are: 

a = -1.090 ± 0.005(s*ai)tg-g?jj(s2/s*); 

b = 0.124 ± 0.006(stat) ± OMO(syst); 

d = 0.057 ± O.Q06(stat)t^ll(syst); 

f = 0.14 ± O.Ol(stat) ± 0.02(syst). 

Our fit raises several points unexpected from previous, low statistics analyses. In 
particular: (i) the fitted value for the quadratic slope in Y, b, is almost one half of 
simple Current Algebra predictions (b = a 2 /4), indicating that higher order corrections 
are probably necessary, (ii) the quadratic term in A, d, is unambiguously found different 
from zero, and (iii) the same applies for the large cubic term in Y. 

We have also fitted the Dalitz plot with a different parametrization of the event density 
which takes into account the final state ir-ir rescattering. Since strong interactions are 
expected to mix the two isospin 1=1 final states of the •;/ — > 37r decay, it is possible 
to introduce a unique rescattering matrix R which mixes the corresponding 1=1 decay 
amplitudes [107]. Thus, we have made a fit using the alternative parametrization of the 
Dalitz plot density distribution, from which it is possible to extract the Dalitz plot slope 
a parametrizing the 37r° decay amplitude. The result a = —0.038 ±0.002 is in reasonable 
agreement with the direct measurement by KLOE of a = -0.027 ± O^l^g [108]. 
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While the polynomial fit of the Dalitz plot density gives valuable information on the 
matrix element, some specific integrated asymmetries as defined in fig. 39 are very sen- 
sitive in assessing the possible contributions to C-violation in amplitudes with fixed AI. 
In particular left-right asymmetry tests C-violation with no specific AI constraint; quad- 
rant asymmetry tests C violation in AI = 2 and sextant asymmetries tests C violation 
in AI = 1. 



L/R Q S 




Fig. 39. - Definition of the r\ Dalitz plot asymmetry tests: left Alr, center Aq, right As 

In the following we present results on asymmetries which use four times the statistics 
entering the Particle Data Group's fits and are the most stringent tests to date. We 
first obtained from MC the efficiency for each region of the Dalitz plot as the number 
of events reconstructed in a particular region divided by the number of events generated 
in the same region, thus taking resolution effects into account. Since no asymmetries 
were introduced in the Monte Carlo, we have checked that the asymmetries estimated 
with our technique in Monte Carlo are compatible with zero. We then evaluated the 
asymmetry on data by counting events in the regions, subtracting the background and 
correcting for efficiency . The measured asymmetries are as follows: 

A LR = (0.09 ± 0.10(stat.)tg;$|(syst.)) x 10~ 2 

A Q = (-0.05±0.I0(stat.)t^(syst.)) x 10~ 2 

A s = (0.08±0.10(stat.)lo:?|(syst.)) x 10~ 2 

9'1.3. rj — > 7r°77. The decay r\ — > 7r°77 is particularly interesting in chiral perturbation 
theory. There is no 0(p 2 ) contribution, and the 0(p 4 ) contribution is small. At KLOE, 
the decay can be reconstructed with full kinematic closure and without complications 
from certain backgrounds present in fixed-target experiments, such as ir~p — > 7r°7r°n. 
Using 450 pb _1 of data, we obtain the preliminary result BR(ij — * 7r°77) = (8.4 ± 
2.7 ± 1.4) x 10~ 5 , which is in agreement with several 0(p 6 ) chiral perturbation theory 
calculations discussed in a recent 77 workshop [109]. The final measurement using the 
complete set of data is forthcoming in 2009. 

9*1.4. 77 — > 7r + 7r~ e + e~. This reaction is interesting because it can test a par- 
ticular formulation [110] of CP violation beyond the SM, through the measurement of 
the angular asymmetry between the pions and electrons planes. The analysis is based 
on the reconstruction of the invariant mass of four charged particles recoiling against 
a monochromatic photon in <p — > 777 events. The main background is represented by 
<f> — > 7r + 7r~7T° events, with either the ir° undergoing a Dalitz decay or one of the two 
photons from the ir° decay converting in the beam pipe. The same happens for <ft — > 777, 
with 77 — > 7r + 7T~7r° with the ir° Dalitz decay, and for <fi — > 777 with 77 — > 7r + 7r _ 7 events 
with photon conversion. These backgrounds are well reduced by kinematic cuts. Signal 
and background counts are extracted from the fit of the invariant mass of the four tracks 
events, using MC spectra for the signal and the various backgrounds, fig. 40. 
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Fig. 40. - Fit of the invariant mass of the four tracks selected. Dots, data; line, total MC. On 
the left signal MC is dark. On the right the hatch corresponds to — > 77-7 backgrounds; the 
continuous line is for all other sources of background. 



Using a data sample of ^622 pb" 1 , about 700 signal events have been observed while 
previous experiments collected less than 20 events. A preliminary estimate of the BR is 
(24 ± 2 stat ± 4 syst ) x 10 -5 . Since approximately 3 times more statistics is available, the 
CP violation test aforementioned will be accessible [111]. 

9'2. rj meson. - The magnitude of BR(0 — > 7/7) is a probe of the ss content of 
the 77'. Furthermore the ratio R^ = BR(0 — > ?/7)/BR((/> — * 777) can be related to the 
pseudoscalar mixing angle tpp in the basis {| uu + dd) /-\/2, | ss}}, offering an important 
point of comparison for the description of the 77-77' mixing in extended chiral perturbation 
theory [112]. 

9*2.1. 77' decays. At KLOE , the <p — ► 7/7 decay is identified either by the final state 
with 77 and charged pions, such as 77' — > r]i{ + Ti~ with 77 — ► 3ir°, or by the 77 and neutral 
pions final state, 77' — ► rjir°ir° with 77 — > 7r + 7r~7r°. The <f> — > 777 decay in the channel 
where 77 — ► 7r + 7r~7r which, as we have learnt from the previous section, is practically 
background free, is used as the normalization process. 

To select <f> — > rj'j events, we require, in addition to seven prompt photons, a vertex 
in a cylindrical region around the interaction point formed by two oppositely charged 
tracks. Then we perform a kinematic fit requiring energy-momentum conservation and 
times and path lengths to be consistent with the speed of light for photon candidates. 
From about 1.4 billion <f> collected, the final number of 77' events from the two selected 
decay chains is, after background subtraction, N v > 7 = 3405 ±61. The selection efficiency 
for detecting the events, evaluated from MC simulations and checked extensively with 
data sub samples, is e v > = (23.45 ± 0.16)%. 

We obtained = (4.77 ± 0.09 stat ± 0.l9 sys ) x 10~ 3 from which we derive BR(<p -> 
7/7) = (6.20 ± 0.11 s t a t ± 0.25 sys ) x 10 -5 . Scant knowledge of intermediate branching 
ratios prevents us to benefit from having twentyfold statistics increase over our previous 
BR determination [18]. 

9'2.2. Pseudoscalar mixing angle. The value of can be related to the pseudo-scalar 
mixing angle. The 77 — 77' system can be parametrized in terms of just an angle only in 
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the quark-flavor basis: 

\t] > = cosipp\qq > +sin(pp\ss > 
\rf > = — sinfp\qq > + cosipp\ss > 

where \qq >— (1/V2)\uu + dd >. Using the approach of ref.[113, 114], where the 
SU(3) breaking is taken into account via constituent quark mass ratio m s /fh, can be 
parametrized as 



(85 ) R , = BR ^ H til = cot2 vp U m > °™ tan ^ ( «V 



BR(</> — > 777) \ m C's sin2iy5p / \p n 

where tpy = 3.4° is the mixing angle for vector mesons; Pr,^') is the 77(77') momentum 
in the <fi center of mass; the two parameters Cms and C's take into account the OZI- 
rule effect on the vector and pseudoscalar wave-function overlap, they are 0.91±0.05 and 
0.89±0.07 respectively, for m a /m=1.24±0.07 [114]. Then we obtain: 



(86) 4> P = (41.4 ± 0.3 stat ± 0.7 sys ± 0.6 th )° 

The theoretical uncertainty on the mixing angle has been evaluated from the maximum 
variation induced from the spread of the ratio m s /fh, Cns and Cs values. In the tra- 
ditional approach to the mixing, the f] — r\' meson are parametrized in the octet-singlet 
basis; in this basis the value of the mixing angle becomes: dp = ipp — arctan\/2 = 
(-13.3±0.3 stat ±0.7 sys ±0.6 th )°. 



9'2.3. 7/ gluonium content. The QCD involves quanta, the gluons, which are expected 
to form bound states, and mix with neutral mesons. While the 77 meson is well understood 
as an SU(3)— flavor octet with a small quarkonium singlet admixture, the 77' meson is 
a good candidate to have a sizeable gluonium content. If we allow for a 77' gluonium 
content, we have the following parametrization: 

(87) |t/ >= X v i \qq > +Y v ,\ss > +Z V > \gluon > 

where the Z v i parameter takes into account a possible mixing with gluonium. The 
normalization implies X"i, + Y^, + = 1 with 

X v i = cos 4>g sin ipp 

(88) Y v > — cos (j)G cos fp 

Z v > = sin 4> G 

where 4>g is the mixing angle for the gluonium contribution. A possible gluonium content 
of the 77' meson corresponds to a non zero value for Zz, , that implies 



(89) 
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and eq. 85 becomes: 



(90) i? = cot 2 ip P cos 2 (j) G 1 



m s C N s tan (g v \ I p v 



m C's sin 2ipp J \ p v 

Combining our result with other experimental constraints, we estimate the glu- 
onium fractional content of rjf meson as Z 2 = 0.14 ± 0.04 and the mixing angle ipp — 
(39.7 ±0.7)° [115]. 

9'3. Scalar mesons, /o(980) and ao(980). - The composition of the scalar mesons with 
masses below 1 GeV, <r(600), /o(980) and do (980), is not yet well understood. KLOE is 
well suited for this task, since the radiative (j) decays into two pseudoscalar mesons are 
dominated by a scalar meson exchange, <j) — > Sj, with S = a, /o,ao, and we have now 
collected literally millions of such decays. A detailed proposal on how to elucidate the 
nature of these mesons with KLOE was one of the first papers JLF had written, after her 
arrival to LNF in 1991 [116]. Since then, theoretical models explaining how to describe 
these scalars and how to calculate their decay rate have proliferated greatly. For this 
review, we only pick two models which are illustrated below. The investigation on these 
mesons is experimentally carried out by measuring branching ratios and studying the 
mass-spectra or the event density in the Dalitz plot. 

Branching Ratios: The branching ratio for the decays 4> — > /07 — > 7T7T7 and <j) — > 
fl o7 ~^ nir 0r y are suppressed unless the /o and ao have significant ss content. The BRs 
for these decays are estimated to be of ~ 10~ 4 if the /o and ao are qqqq states, in which 
case they contain an ss pair, or of ~ 10~ 5 if the /o and ao are conventional qq states. If 
the /o and ao are KK molecules, the BR estimates are sensitive to assumptions about 
the spatial extension of these states. 

Invariant Mass distributions: Fits to the irir (rjir ) invariant-mass distributions can 
also shed light on the nature of the /o (ao) since they provide an estimate of the couplings 
of these mesons to the </> , g^s-f, and/or to the final-state particles. The rate expression 
for <fi — > 57 El-transitions contains a factor , where E 1 is the energy of the radiated 
photon, as required by considerations of phase space and gauge invariance. As a result, 
the invariant-mass distributions are cut off above and develop a long tail toward 
lower mass values. The fit results therefore strongly depend on the scalar meson masses 

Kaon loop 
Y 



f a -/^ 

\ w ^^""I'l V y 

K / , aT\ 71 No structure 




Fig. 41. - Feynman diagrams used to describe the tytv and nit mass spectra for the cf> — > Sj 
decay in the kaon-loop (left) and No structure (right) models. 



and widths and on the assumed specific model for the decay mechanism. Because of 
the proximity of the fo and ao masses to the KK threshold, and because these mesons 
are known to couple strongly to KK, the kaon-loop model, fig. 41 left, is often used 
[117]. In this model, the decay proceeds through a virtual K + K~ pair emitting the 
photon and subsequently annihilating into a scalar. This loop function damps the 
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behavior. The transition amplitude depends on the coupling of the scalar meson both to 
the pseudoscalar meson in the final state and to the K + K~ pair in the loop function. The 
propagator include finite width corrections which are relevant close to the KK threshold. 
A second model, called no-structure [118], sketched in fig. 41 right, describes the process 
as a point-like interaction where the dynamics of the scalar production is absorbed in 
the g^s-y coupling. The amplitude is described by a Breit-Wigner propagator, with a 
mass dependent width, which accounts for analytical continuation under irir and KK 
thresholds, and a complex polynomial describing a continuum background. 

9'3.1. (f> — > 7r°7r°7 and <fi — ► ?77r 7 decays.. The first KLOE studies of the <fi — ► 7r°7r°7 
[16] and <p ~> r/7r°7 [17] decays were performed with the 17 pb _1 of data collected in 
2000. All BRs quoted therein are in agreement with, but are much more precise, than 
earlier measurements. 

The amplitudes contributing to <p — * 7r°7r°7 include <j> — ► S-y, with S = fo or er, and 
4> — ► /0°7r° with p° — ^ 7r°7. The final state was selected requiring five prompt 7 in the 
event and pairing the photons to get the best value for the pion masses. The energy 
resolution was improved by performing a kinematic fit to the event. In this final state, 
more than 50% of the events were due to the not-resonant e + e~ — > un process, with 
uu — > 7r°7. This background was largely removed by directly cutting around m w on the 
reconstructed ny invariant mass. After this rejection, the M„ distribution was fit with 
a function including in the amplitude the terms describing the process <j) — > S7, <fi — ► pn, 
and their interference. We obtained a BR(> -> 7r°7r°7) of (1.09 ± 0.03 ± 0.05) x 10~ 4 . 

In analogy, for the ao case, the amplitudes contributing to the process are <j) — * a y 
and 4> — > p 7T°, with p° — > 777. The not-resonant process e + e~ — » W7r° with w — > 777, 
constitutes a practically negligible background. The final states corresponding to two 
different rj decay channels were studied: rj — > 77 (57 final state) and r\ — » 7r + 7r _ 7r° 
(7r + 7r _ 57 final state). The distributions of both samples were fit simultaneously with 
a function analogous to that discussed above. The values obtained for BR((/> — > rjir^y) 
were (8.51 ± 0.51 ± 0.57) x 10~ 5 and (7.96 ± 0.60 ± 0.40) x 10" 5 , for the 5 7 and tt+tt-^ 
samples respectively. 

The analyses of the <j> — > 7r°7r°7 and — » t]-k°7 decays based on a later data set, benefit 
of an increase in statistics by a factor of ^30 such that we could improve background 
removal and analysis cuts. 

For the 7r°7r°7 final state, the analysis is done without cutting on the ir°y invariant 
mass, which could bias the residual M W7r spectra, and plotting all events with same 
topology in the Dalitz plots shown in fig. 42 left. Two clear bands appear related to 
the not-resonant process which constitutes the larger background to the scalar term for 
M wvr < 700 MeV. Above this threshold, the /o(980) contribution is almost background 
free. We estimate with theory the event density on the Dalitz-plot by adding to the 
amplitudes 4> — > £7 and <f) — ► pir° the complete VDM description of the not-resonant 
process. The interferences between all diagrams are also taken into account [119]. We 
then fit the Dalitz plot density, by folding theory with reconstruction efficiencies and 
with a smearing matrix describing the probability for an event to migrate from a bin 
to another one. In particular, we fit data with (1), an improved k-loop model [120], 
and (2) with the no-structurc model [118]. In the improved k-loop model,the fo and a 
description are strongly coupled with the energy dependence of the rnr scattering phase. 
Ten different parametrizations are introduced. We find that, without adding the a in 
the model, it is impossible to adequately describe our data-set. In fig. 43 the function 
describing M T7I , as extracted by our best fit result for the k-loop model, is shown. The 
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Fig. 42. - left, M,^,,^ Dalitz plot in logarithmic scale after background subtraction. The two 
bands in the low region represent the e + e~ — » um° background, right, a e + e -_ fv o 7r o 1 near 
the /o mass region. 



contribution of each process is also represented. A large coupling of the /o(980) to kaons 
is obtained with large theory errors related to the parametrization choice (see table XII). 
Discussing with the authors of ref. [120], some typos in the parametrization were found. 
After correcting for it, a new preliminary set of values is extracted which much reduced 
theory errors, referred as 7r°7r°7 update in table XII. All of the above indicates a not- 
conventional structure for the /o meson. For the no-structure model, the situation is 
really different. The /o coupling to kaons get reduced and the introduction of the a does 
not substantially modify the results. However, the physical interpretation gets somehow 
hidden from the introduction of the continuum background in the parametrization. 



Parameter 


7T + 7T 


7 


7r°7r°7 


7r°7r°7 updates 


m fo (MeV) 


980^ 


987 


976.8 ±0.3 flt lag S ys + 10.1th 


984.7 ± 1.9th 


9f K+K- (GeV) 


5.0^ 


6.3 


3.76 ± 0.04fl t s y s -o'.Js th 


3.97 ± 0.43th 


g fon+ „- (GeV) 


3.0 -f 


4.2 


i Ar> _i_ n ni +0.01 +0.03 

— 1.4,3 ± U.Ulflt _ 06 sys _ .60 th 


-1.82 ± 0.19th 


J f K+K- 


2.2 -f- 


2.8 


R a 4- n 1 +°- 2 +° 3 

D.y =t U.iflt _Q.i sys _3.g th 




9f K+K- (GeV) 


1.6-=- 


2.3 


0.40 ± 0.04fl t toU sys 




g fo „ +1T - (GeV) 


0.9^ 


1.1 


1.31±0.01 fit t° f 3 sys 




"f K+K- 


2.6- 


4.4 


0.09±0.02 flt t° :tt sys 




50/07 (GeV -1 ) 


1.2 -+ 


2.0 


2.61 ± 0.02fl t lo'.08 sys 





Table XII. - Comparison of fo parameters estimates between k-loop and no-structure models. 
In the last column, we report the preliminary update of the KL fit to the 7r°7r°7 final state. 
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Fig. 43. - M W07TQ spectrum coming from the KL fit result. 

By integrating the scalar term contribution, we determine a value for the BR which 
slightly depends on the theory model used in the fit. We quote a BR to (f> — > 7t 7t°7 [119] 
of: 

BR(0 5 7 7T°7T° 7 ) = (l.07t°;°J(fit) tg$(syst) t™l(mod)) x 10" 4 . 

For the r]n°j final state, same analysis strategy of the previous work was followed with 
the advantage of using a more realistic detector response and simulation of the machine 
background in the official KLOE Monte Carlo. The physics background was estimated 
after event preselection by fitting MC shapes to data. At the end of the analysis, we have 
a signal efficiency of ~ 40% (20%) and a residual background of 55% (15%) for the sample 
with ?7 — > 7T (?7 — > 7t + 7t~7t°). The absence of a major source of interfering background 
allows to get the BR directly from event counting. We get a preliminary value for the 
BR(<t5 -> r/7T°7) of (6.98±0.10±0.23 S!/S t) x 10" 5 and (7.05±0.10±0.21) x 10~ 5 , for the 5 7 
and 7r + 7T~57 samples respectively, which corresponds to a decrease of the central value 
with respect to our previous measurement of ~ 12% with an improvement in precision 
better than a factor 3. In fig. 44, the rjir invariant mass distribution, M^o, is shown 
for both final states. A simultaneous fit to the spectra is also performed to extract the 
relevant ao coupling. Also for this meson, a large coupling to kaons is found. 

9'3.2. 0^7T + 7T _ 7. KLOE has also published a study of the decay <j> — > / o7 — » n + -K~-f 
[121]. The /o's charged mode was extremely difficult to identify because only a small 
fraction of the e + e~ — ► tt + 7T~7 events involve the /o; the principal contributions are 
from events in which the photon is from ISR or FSR. The analysis is performed on the 
M nv distribution for events with large photon polar angle where the ISR contribution 
is much reduced. The M W7r distribution was fit with a function composed of analytic 
expressions describing the ISR, FSR, pir contributions and two terms that describe the 
decay <f) - * 5*7 - * 7r7r 7 an d its interference with FSR which could be constructive or 
destructive. 

Fig. 45 shows the M M distribution, with the result of the kaon-loop fit superimposed. 
The overall appearance of the distribution is dominated by the radiative return to the p; 
the /o appears as the peak-like structure in the region 900-1000 MeV. Both the kaon-loop 
and no-structure fits strongly prefer destructive interference between the Sj and FSR 
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Fig. 44. - Distribution of M V7T o for cj> — » riiv 'y events. Right. ir + ir tt°tt° final state. Left. Five 
photon final state. Points, data; solid line, fit result. 



amplitudes. The kaon-loop fit gives coupling values in reasonable agreement with those 
obtained from the fit to the KLOE <fi — > 7r 7r°7 data discussed before. Integrating the 
appropriate terms of the fit functions gives values for BR(0 — * /07) in the neighborhood 
of 2 x 1CT 4 . 

Another important test in this final state is the pion forward backward asymmetry. 
The 7r + 7r~ pair has a different charge conjugation eigenvalue depending on whether it 
is produced from FSR and /o(980) (C = +1) or ISR (C = —1). An interference term 
between two amplitudes of opposite charge conjugation gives rise to C-odd terms that 
change sign by the interchange of the two pions and results in an asymmetry A c : 

_ A(6> 7r + > 90°) - N{8^ < 90°) 
c ~ N(9^+ > 90°) + N(6^ < 90°) ' 



Fig. 46 shows the comparison of the asymmetry data with the expectations based on 
the KL parametrization of the /o(980). [122, 123] The A c data distribution is reproduced 
only by calculations including the /q(980). 




600 800 1000 

Fig. 45. - Distribution of M„+ v - spectra for <f> -—> n + n~"/ events. 
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Fig. 46. - Pion forward backward asymmetry distribution. Data: black dots; Simulation 
FSR+ISR: open triangles; Simulation FSR+ISR+ scalar(KL): open squares. Right: zoom of 
the /o mass region. 



9'3.3. 4> — ► KKj. A fundamental consistency check of the study of <f> radiative decays, 
as well as being another important tool for the investigation of the s quark content of 
0o (980) and ao(980), is the measurement of final states involving kaons. The search of 
4> — > KK'-f, which is expected to proceed mainly through the [/o(980) + a (980)]7 in- 
termediate state, adds relevant information on the scalar meson structure. Theoretical 
predictions for the BR spread over several orders of magnitude, although the latest eval- 
uations essentially concentrate in the region of 10 -8 (see [124] and references therein). 
Experimentally, this decay has never been observed. Using 1.4 fb _1 collected at the cj) 
resonance and an equivalent MC statistics to model the background, the KLOE exper- 
iment searched for this decay using the <f> — > K$Ks"f — > ir + ir~ ti + it~ 7 decay channel 
[124]. The significant signal reduction (24% of the total rate) is compensated by a very 
clean event topology. The MC signal has been generated according to phase-space and 
radiative decay dynamics. The selection cuts applied to reject the main background 
component {<f> — > K$Kl with an ISR or FSR photon) have been optimized using MC 
only: no events survive the analysis cuts, while one event is found in data. The resulting 
upper limit on the branching ratio of <j> —> KsKsj < 1.8 x 10 -8 @ 90% CL rules out 
most of the theoretical predictions. 

9"4. Vector mesons. - 

9'4.1. <j) leptonic widths. The 0-leptonic widths provide information on the (^-structure 
and its production cross section in annihilations. They are necessary for decay branch- 
ing ratio measurements and estimates of the hadronic contribution to vacuum polar- 
ization. There is no direct measurement of the leptonic width. The cross sections for 
processes e + e~ — > e + e~, can be written as o~ ee — a ee .^+Gi nt and = 07^,0+0-;^, 

with ai n t, the interference term, given by: 

(91) (Tint = TjTTo 2 \2 , w2r2 / /«(6 l )dcoS0. 

m {W 2 - ml ) 2 + W 2 Tl J cos dmm 
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where W is the energy in the collision center of mass (CM), 9 min and 9 max define the 
acceptance in the polar angle 9 (see later) and where T« = T ee for e + e~ — ► e + e~ and 
Tu = \Zr ee r MAl /£ for e + e~ — » The £ term takes into account for the phase space 

correction. KLOE used data collected at CM energies W of 1017.2 and 1022.2 MeV, i.e. 
77i0 ± r^/2, and at the <j> peak, PV=1019.7 MeV. For e+e~ — > we measure the 

cross section. Since Bhabha scattering is dominated by the photon exchange amplitude, 
the interference term is best studied in the forward-backward asymmetry, Afb, defined 
as 



(92) A FB = ^-^. 

of + &b 

where up and ob are the cross sections for events with electrons in the forward and 
backward hemispheres. 

We find r ee =1.32±0.05±0.03 kcV and v /r ee r w =1.320±0.018±0.017 kcV. These re- 
sults, compatible with T ee =Tp^, provide a precise test of lepton universality. Combining 
the two results gives T u (4>) = 1.320 ± 0.023 keV[125]. 

9'4.2. </> — ► 7T + 7r _ 7r° . The decay of the <j) meson to 7r + 7r~7r , with a branching ratio 
(BR) of <~15.5%, is dominated by the pn intermediate states [126] p + n~ , p~7r + , and 
p°ir° with equal amplitudes. CPT invariance requires equality of the masses and widths 
of p + and p~ , while possible mass or width differences between p° and p* 1 are related to 
isospin-violating electromagnetic effects. Additional contributions to e+e~ — > 7r + 7r~7r° 
are the so called "direct term", <j> — » 7r + 7r~7r°, and e + e~ — > lutt , u> — > 7r + 7r~. We use 
some 2 million of e+e~ — ► 7r + 7r~7r° events to determine the masses M(p H °) and widths 
T(p H °) of the three charge states of the p-meson, from a fit to the Dalitz plot density 
distribution [127]. If we define the variables x = T+ - T~ and y = T°, where T + • 
are the kinetic energies of the three pions in the center of mass system (CM), the Dalitz 
plot density distribution D(x,y) is given by: 

(93) D(x,y) oc \p* x pl\ 2 \A p7I + A dir + A^\ 2 

where p± are the ^ momenta in the CM and A p7T , ^4dir and A uv are the three amplitudes 
described above and containing the dependence on the masses and widths of the p- 
mesons. 7r + 7r~7r° events are selected by asking for two non-collinear tracks with opposite 
sign of curvature and polar angle 6 > 40° which intersect the interaction region. The 
acollinearity cut (A6 < 175°) removes e + e~7 events without incurring an acceptance loss 
for the signal. The missing mass, M miss = (E^ - - E w - ) 2 - - p w + - p v - | 2 , 
where E and p are laboratory energies and momenta, is required to be within 20 MeV 
of the 7T° mass. This corresponds to an effective energy cut of < 20 MeV on the total 
energy radiated because of initial state radiation (ISR) . Two photons in the calorimeter 
are also required, with opening angle in the 7r° rest frame cos# 77 < —0.98. The Dalitz 
plot variables x and y are evaluated using the measured momenta of the charged pions, 
boosted to the center of mass system: x = — E*_ and y = E^ — E+ — El — m w o = T w o . 
Efj, and p^, are measured run by run using Bhabha scattering events. The resolution on 
x and y is about 1 MeV over the full kincmatical range. Three x 2 fits to the Dalitz plot 
density distribution have been performed: a) a fit assuming CPT and isospin invariance, 
i.e. m p a = m p + = 777 p - , r p o = T p + = T„- ; b) a fit assuming only CPT invariance i.e. 
777 p + = ?Tjp- , T p + = T p - ; and finally c) without assumptions on masses and widths. The 
results of the fits are shown in table XIII. 
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parameter 


fit (a) 


fit(b) 


fit(c) 


x 2 b(x 2 )] 


1939 [12%] 


1914 [21%] 


1902 [26%] 


m p o 




775.9 ±0.5 ±0.5 


775.9 ±0.6 ±0.5 


m p+ 


775.8 ±0.5 ±0.3 


775.5 ±0.5 ±0.4 


776.3 ± 0.6 ± 0.7 


m p - 






774.8 ± 0.6 ± 0.4 


r p0 




147.3 ±1.5 ±0.7 


147.4 ± 1.5 ± 0.7 




143.9 ±1.3 ±1.1 


143.7 ±1.3 ±1.2 


144.7 ± 1.4 ± 1.2 








142.9 ± 1.3 ± 1.4 



Table XIII. - Fit results. Masses and widths are in MeV. The amplitudes are scaled in such a 
way that a p —l. \ 2 values and probabilities of the fits are given in the first row. 



The convergency of fit (a) shows that the experimental distribution is consistent with 
CPT and isospin invariance. The p masses and the widths are close to recent results [21]. 
The results obtained for mass and width differences, including correlations between the 
parameters, are: m p a — m p ± = 0.4 ± 0.7 ± 0.6 MeV,m p + — m p - = 1.5 ± 0.8 ± 0.7 MeV, 
T> -r p ± = 3.6 ±1.8 ±1.7 MeV, T p + -T p - = 1.8 ±2.0 ±0.5 MeV. Note that the neutral 
and charged p masses are essentially equal within statistics. 

9'4.3. (j) ^ uk. Using <~ 600 pb _1 collected at center of mass energies between 1000 
and 1030 MeV, the KLOE experiment studied the production cross section of e+e~ — ► 
7r7r~7r 7r° and e + e~ — > 7r°7r°7 processes [128], which mainly proceed through the non- 
resonant p — ► W7T° intermediate state. The dependence of the cross section on the center 
of mass energy is parametrized in the form <r(W) = cto(W)|1 — Z m ^ \ , where Z is the 
complex interference parameter between the <j> decay amplitude and the non-resonant 
processes, whose bare cross section is ao(W). m^, and D$ represent the mass, the 
width and the inverse propagator of the (f> meson respectively. A model independent 
parametrization is used to describe the increase of the non-resonant cross section with 
W, which is linear in this energy range: <Jo(W) = <t + a'(W — m^). 

In fig. 47 data points with the superimposed fit function are shown for both channels. 

By fitting the interference pattern for both final states, a preliminary ratio T(u> — > 
ir°j)/T(uj — > ir + TT~iT ) = 0.0897 ± 0.0016 is obtained. Since these two final states rep- 
resent the 98% of the w decay channels, this ratio and the sum of the existing BR mea- 
surements on rarest decays are used to extract the main u> branching fractions imposing 
unitarity. BR(w -> tt+tt-tt ) = ( 90.24±0.18 )% and BR(w -» 7r° 7 ) = ( 8.09±0.14 )% are 
extracted, with a correlation of -71%. The parameters describing the e + e~ — ► 7r + 7r~7r°7r° 
reaction around allow to derive the branching fraction for the OZI and G-parity vi- 
olating decay 4> — > u>7r°, which is found to be: BR(</> — > lott ) = (4.4 ± 0.6 ) x 10~ 5 . 

10. The hadronic cross section and a p 

lO'l. The muon anomaly a p . - The observation of the anomalous magnetic moment 
of the electron helped drive the development of quantum electrodynamics (QED). The 
value of the muon anomaly, a pi is (m p /m e ) 2 40,000 times more sensitive than that 
of the electron to high-mass states in the polarization of the vacuum. The Muon g — 2 
Collaboration (E821) at Brookhaven has used stored muons to measure a p to 0.5 ppm 
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Fig. 47. - Cross section fit results for the e + e~ — ► 7r7r~7r 7r° (top-left) and e + e~ — > 7r°7r°7 
(bottom-left) channels. Black dots are data, solid line is the resulting fit function. Branching 
fraction for the two main to decay channels (right). The square is the KLOE fit result, while 
the dot is the PDG constrained fit result. The shaded regions are the 68% CL 



[129]: 

(94) a M = 9ti ~ 2 = 116,592,080 ±60 x HT 11 . 

The muon anomaly receives contributions from QED, weak, and hadronic loops in the 
photon propagator, as well as from light-by-light scattering. The lowest-order hadronic 
contribution is ajj ad = —7000 x 10~ n , with an uncertainty of —60 x 10~ n . To the extent 
that this uncertainty can be reduced, the E821 measurement offers a potential probe of 
new physics at TeV energy scales. 

The low-energy contribution to a{^ ad cannot be obtained from perturbative QCD. 
However it can be calculated from measurements of the cross section for e + e~ annihilation 
into hadrons via the dispersion integral 

(95) a^ d = -3 J2 / K{s) a(e+e- /) ds, 

47r f Jsth(f) 

where K(s), the QED kernel, is a monotonic function that varies approximately as 1/s, 
with s the squared center-of-mass collision energy, see Fig. 48, left. This amplifies the 
importance of the cross-section measurements at low energy. Approximately two-thirds 
of the integral is contributed by the process e+e~ — > n + -K~ for < 1 GeV, i.e., in the 
vicinity of m p . Before the arrival of the KLOE results we discuss below, calculations of 
fl had f rom wQj-Jd e + e _ data (dominated at low energies by the CMD-2 measurement of 
[13]) led to a value of a p —2.4 a lower [130] than the final value reported by E821. 

10'2. Measurement of cr x + T - and a p at KLOE. - DA$NE is highly optimized for 
running at W — and no large variation of W are accessible. However, initial state 
radiation, ISR, via the process e+e~ — > 7r + 7r _ 7 naturally provides access to hadronic 
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Fig. 48. - left: connection between a(mr'y), H, a(nn), K and a M(J . Right, H, radiator function, 
ow 7 with pion form factor set to one. 



states of lower mass. The cross section cr„ 



over the entire interval from threshold to 



rritf, is related to the distribution of s v = for tt + tt 7 events, 



(96) 



ISR 



where H, the "radiator function," describes the ISR spectrum. Note the subscript ISR 
on the differential 7r + 7r~7 cross section. This is very important because the contribution 
from final state photon emission, FSR is of the same order as that from the ISR process 
unless one imposes special fiducial volume cuts to suppress it. 




On the other hand, as is depicted in fig. 49 one must remember eventually to include 
FSR in the final evaluation of aj^ ad . To correctly calculate H and estimate the effects of 
FSR, an accurate simulation of the a- K7ri is critical. The radiator function used by KLOE 
to obtain o^n is provided by the Phokhara code [131] by setting the pion form factor 
F^(Mt T . k ) = 1, see Fig. 48, right. The Phokhara event generator has been continuously 
improved by inclusion of next-to- leading-order ISR (two initial-state photons), leading- 
order FSR and ISR-FSR [132]. 

Furthermore, has to be corrected for the running of the fine structure con- 
stant [133] (vacuum polarization) and for shifting from M W7r to the virtual photon mass, 
Af 7 » for those events with both an initial and a final state radiated photon. The Phokhara 
version of Ref. [134] contains also these contributions and it is used for the evaluation of 
the acceptance correction and for all reconstruction efficiencies. 
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10'2.1. Extraction of a™ and pion form factor. The differential cross section of 
e+e~ — > 7T7r7 as a function of M„ is obtained by subtracting, bin by bin, the back- 
ground events N^kg from the number of observed events N b s , correcting for the selection 
efficiency e se ;(M^ 7r ), and finally dividing by the total luminosity obtained as described 
in sect. 3'3: 

da-KTTJ _ N b s — Nbkg 1 

where our mass resolution allows us to have bins of width AM^ = 0.01 GeV 2 . 

The background content is found by fitting spectra of selected data sample with a 
superposition of Monte Carlo distributions describing the signal and background sources. 
The only free parameters of such fits are the relative weights of signal and backgrounds 
in data, computed as a function of M nn . 

The " bare" cross section cr^^ which is obtained after applying the corrections writ- 
ten in eq. 96, inclusive of FSR, is used to determine aj^: 

(98) c4r ds * (s) ^' 

•'Slam 

where the lower and upper bounds of the spectrum depend on specific measurement 
configurations. 

The measured avw is related to the pion form factor via the relation cr 7r7r (s) = 
7^-/3^ \F n (s)\ , where s is the center-of-mass energy squared, is the pion mass and 
(in — y/l — 4771^ / s is the pion velocity in the center-of-mass frame. Since |F„-(s)| 2 is 
measured by many experiments, it is customarily used for spectra comparisons. 

Two different selection regions are used: the first one is named "Small angle 7 selec- 
tion" (SA) since photons are required to be within a cone of 8 7 < 15° around the beam 
line (narrow cones in fig. 50, left), the second one is called "Large angle 7 selection" 
(LA) since there should be at least one photon at a polar angle of 50° < Q 1 < 130° (large 
central cones in fig. 50, left). In both cases the two charged pion tracks should have 
50° <9 n < 130°. 

10'2.2. Small angle 7 selection and data analysis. At small photon polar angle 7 , 
ISR events are vastly more abundant than FSR events. Requiring angular separation 
between the pions and the photon further suppresses FSR [135]. Thus, the first KLOE 
measurement [9] was performed using this event selection criterion. Since ISR-photons 
are mostly collinear with the beam line, a high statistics for the ISR signal events remains. 
On the other hand, a highly energetic photon emitted at small angle forces the pions also 
to be at small angles (and thus outside the selection cuts), resulting in a kinematical 
suppression of events with < 0.35 GeV 2 . 

In this analysis the photon is not explicitly detected. Its direction is reconstructed 
from the tracks' momenta by closing kinematics: ~ Pmiss — V<t> ~ (Ptt+ +Ptt-)- The 
separation between pion and photon selection regions greatly reduces the contamination 
from the resonant process e + e~ — > <j> — ► 7r + 7r~7r°, in which the ir° mimics the missing mo- 
mentum of the photon(s) and from the final state radiation process e+e~ — ► 7t + 7t~7fsr- 

Discrimination of 7r + 7r~7 from e + e~ — > e + e _ 7 events is done via particle identification 
based on the time of flight, on the shape and the energy of the clusters associated to 
the tracks. In particular, electrons deposit most of their energy in the first planes of 



(97) 




the calorimeter while minimum ionizing muons and pions release uniformly the same 
energy in each plane. Events with at least one of the two tracks not identified as electron 
are selected. This criterion results in a rejection power of 97% for e + e~j events, while 
retaining a selection efficiency of ~ 100% for tt + it~j events (see fig. 51, left). 

Contaminations from the processes e + e~ — > /i + /i~7 and 4> ~~ > 7r + 7r~7r° are rejected 
by cuts on the track mass variable, mtrk, defined by the four- momentum conservation, 
assuming a final state consisting of two particles with the same mass and one photon, 

and on the missing mass, m m i SS = \J E 2 X 
e + e~ — > 7r + 7r~ X. 



= y E x — |-Px| 2 , defined assuming the process to be 
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Fig. 51. - Left: e+e— 7 suppression using particle ID. Right: Signal and background distributions 
in the M tr k — M^-plane. 



KLOE first published in 2005 results obtained with 141 pb^ 1 collected in 2001 [9], 
which we will refer to as KLOE05 in the following. The statistical errors ranged from 
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~2% at the lower limit in s T to ^0.5% at the p peak. The experimental systematic 
uncertainties were mostly flat in s T and amounted to 0.9%. The luminosity was esti- 
mated using large-angle Bhabha scattering and contributed an additional, dominantly 
theoretical uncertainty of 0.6%. The cross section <J v + n - was obtained via eq. 96. The 
contribution to the systematic uncertainty on <7 T + T - from FSR was 0.3%, and the theory 
accuracy of H was 0.5%. Finally, for the evaluation of a{^ ad , it was necessary to remove 
the effects of vacuum polarization in the photon propagator for the process e + e~ — > tt + tt~ 
itself by correcting for the running of a em . This contributes 0.2% to the systematic un- 
certainty. Thus from KLOE05 analysis we obtained for the contribution to a™ in the 
energy interval 0.35<s<0.95 GeV 2 : 

(99) 0)7(0.35 < s < 0.95 GeV 2 ) = (388.7 ± 0.8 stat ± 3.5 syst ± 3.5 th ) x 10~ 10 . 

Inclusion of the KLOE05 data together with the CMD-2 data in the evaluation of 
oJj ad increased the discrepancy between the calculated and measured values of a M to 
(25.2 ± 9.2) x 10~ 10 (^2.7a) [136]. Furthermore, the precise KLOE spectrum exhibited 
a pronounced deviation from those derived from r lepton's spectral functions that led to 
the later being excluded in a M compilations [136]. 

KLOE in 2008 completed the analysis of 241 pb" 1 data taken in 2002, which we will 
refer to as KLOE08 in the following. These data benefit from cleaner and more stable 
running conditions of DA<I>NE resulting in less machine background. In particular the 
following changes are applied with respect to the data taken in 2001: i) an additional 
trigger level was implemented at the end of 2001 to eliminate a 30% loss due to pions 
penetrating up to the outer calorimeter plane and being misidentified as cosmic rays 
events; ii) the offline background filter, which contributed the largest experimental sys- 
tematic uncertainty (0.6%) to the published work [9], has been improved and includes 
now a downscale algorithm providing an unbiased control sample. This greatly facilitates 
the evaluation of the filter efficiency which increased from 95% to 98.5%, with negligible 
systematic uncertainty; iii) in addition, the knowledge of the detector response and of 
the KLOE simulation program has been improved. Thus the relative systematic errors 
on the extraction of a™ in the mass range [0.35,0.95] GeV 2 decreased from 1.3% for the 
2001 to 0.9% for 2002 data. Finally, we want to stress that the dominant term contribut- 
ing to the systematic error in KLOE08 analysis is given by the 0.5% uncertainty quoted 
by the authors of Phokhara [134]. The KLOE08 results are shown in fig. 52. On the left 
the number of tt + it~j events per 0.01 GeV 2 are shown. The ir + ir~-f differential cross 
section is shown in the right. Finally, while performing the 2008 analysis we revisited 
the published 2001 data analysis, named KLOE05 updated, and we found a bias in the 
evaluation of the trigger correction that affects mostly the low M ra region. Correcting 
for this effect and normalizing to the new Bhabha cross section we updated the published 
spectrum to compare with KLOE08. Fig. 53 , left shows the pion form factor evaluated 
from these two sets of data. Comparison among a* n Q0.35,0.95]GeV 2 ) values in units of 
10 10 evaluated with the KLOE small angle 7 selection analyses: 

KLOE05: 388.7 ± 0.8 stat ± 3.5 sys ± 3.5 th ; 

KLOE05 Updated: 384.4 ± 0.8 stat ± 3.5 sys ± 3.5 th ; 

KLOE08: 387.2 ± 0.5 stat ± 2.4 sys ± 2.3 th . 

We see that in terms of the net shift of the updated value is about one standard 
deviation below the KLOE05 one and it is also in excellent agreement with the KLOE08 
analysis. 
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Fig. 52. - Left: KLOE08 tt+tt^ events per 0.01 GeV 2 . Right: KLOE08 7T + 7r _ 7 differential 
cross section 



Finally, we make a comparison of the most recent aV* evaluations released by the 
CMD-2 [137] and SND [138] experiments with the KLOE08 [139], in the mass range 
e [630, 958] MeV. 
CMD-2: 361.5 ± 1.7 stat ± 2.9 sys - 
SND: 361.0 ± 2.0 stat ± 4.7 sys ; 
KLOE08: 356.7 ± 0.4 stat ± 3.2 sys ; 

The CMD-2 and SND dispersion integrals are performed with the trapezoid rule, the 
KLOE value is done extrapolating our bins to match the [630,958] MeV mass range, and 
summing directly the bin contents of the dcr 7r7r7 differential spectrum, weighed for the 
kernel function. 

The CMD-2 and SND values agree with the KLOE08 result within one standard 
deviation. Thus, inclusion of the latest KLOE results in the world average definitely 
established without a doubt the discrepancy between the measured and calculated values 
of a leading to hopes that there is a glimpse of which direction we should search for 
new physics. 

10'2.3. Large angle 7 selection and data analysis. The contribution to a{j ad from 
e + e~ — > 7r + 7r~ in the interval in s between threshold and 0.35 GeV 2 is approximately 
1000 x 10~ n . To explore the low-mass part of the spectrum, KLOE selects events 
requiring a photon in the calorimeter with energy greater than 50 MeV and 50° < B 1 < 
130°. With the photon explicitly detected, about 40% of the background from <f> — ► 
7r + 7r~7r° events is rejected by kinematic closure. A further variable, the 3-dimcnsional 
angle between the missing momentum and the photon, which is peaked at degrees 
for signal events, removes most of the remaining <f> — > 7r + 7r~7r° background. 

There are however two sources of irreducible background : FSR and <p — > /07 — > 
7r + 7r~7. FSR here is a very complicated problem: for 6 1 > 40° ISR and FSR events 
contribute nearly equally to the tt + tt~j spectrum in the tails of the p, and the accuracy 
of the generator used to obtain FSR corrections is critical. 

As discussed by Binner et al. in Ref. [135], the ISR-FSR interference results in a mea- 
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- Left: KLOE05 updated form factor in comparison with KLOE08. Right: comparison 
KLOE 2008, SND 2006, and CMD 2007 form factors. 



surable charge asymmetry, as well as a distortion in the mass spectrum. The measured 
charge asymmetry can be used as a useful gauge of the accuracy of the simulation. KLOE 
has performed comparisons of this type in the analysis of 4> — > Tt + n~j as discussed in 
the previous section on scalar mesons [139]. 
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Fig. 54. - left: mr cross section from small angle 7 selection, right: comparison between small 
and large 7 angle selections. 



We have made a measurement of cr W7r with photon emitted at large angle using the 
data collected in 2002. Since it is an independent event selection scheme, by comparing 
this result with the small angle 7 analysis, we can also test the present models of FSR 
contributions. Fig. 54 shows both spectra; note that the 7T7r threshold is reachable 
with the large angle selection. To make a quantitative comparison, at present we are 
limited by the unknown interference between the FSR process and the resonant decays 
4> — > ir + ir~~f as mentioned before. Therefore we limit the comparison to the range 



84 



F. BOSSI, E. DE LUCIA, J. LEE-FRANZINI, S. MISCETTI, M. PALUTAN ETC. 



where resonance contributions are negligible, that is, a M to [0.5,0.85] GeV 2 where from 
the "small angle selection" we obtain 255.4 ± 0.4 stat ± 2.5 sys and for the "large angle 



selection", 252.5 ± 0.6, 



± 5.1, 



in excellent agreement. 



The main source of systematic uncertainty in the large 7 angle selection is the fo (980) 
background subtraction. 

10'2.4. Pion form factor from R. We can also obtain the pion form factor directly 
from measuring the ratio, R, bin-by-bin of the observed pion and muon (radiative) cross 
sections. In the absence of photons from final state radiation, one finds [140] 



(100) 



4(1 



7isr(7)) 



1% 



d<r(n[i + tisr(t)) ' 



where is the muon mass, and /3 M , f3„ are the muon and pion velocities in the center- 
of-mass frame. 

The unavoidable presence of photons from final state emission of both the pions and 
the muons in the data leads to a deviation from the "ideal" eq. 100 which requires 
treatment of final state radiation (FSR) with care. The advantage of this method is 
that one does not need an absolute luminosity measurement, thus providing a totally 
independent measurement of the pion form factor from those obtained using small or 
large angle 7 selections. 

The experimental challenge of measuring the pion form factor using the ratio is the 
overwhelming cross section over tttt cross section at low momentum, see fig. 55, left. 
Furthermore, there is an overlap region between the 7r7r7 events and /X/X7 events which 
makes precision event counting/separation a delicate operation, fig. 55, right, shows the 
separation between 7r + 7r~7 and /i + /i _ 7 events, achieved selecting different mtrk regions 
around the missing mass peaks. The residual contamination of ■k + tt~tt° events are seen 
at high m tr k values. 
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Fig. 55. - Left: mr versus /i^i cross section as a function of energy. Right: . Discrimination of 
the ^t + /i~7 (filled left histogram) from the 7t + 7t~7 (filled right histogram) events. 



We are developing additional tools for ir/j, separation and the results from extracting 
|Ftt(s)| 2 spectra from using the ratio is extremely promising [141]. We are confident of 
obtaining aJJ 7 from yet another independent method by the end of 2008. 
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11. — Conclusions and outlook 

The contributions by KLOE to particle physics during the years 2002-2008 are of 
great importance. The study of the decays of the kaons have demonstrated unitarity of 
the CKM matrix with accuracy better than one per mil and universe ility of the couplings 
between leptons and quarks. New mass and lifetime values have been obtained. Studies 
of the T) meson have reached new levels of precision. KLOE has made the most accurate 
measurement of the hadronic cross section pertaining to the calculation of the muon 
anomaly, and has made the most detailed studies on the nature of light scalar mesons. 

A summary of all the masses, lifetimes, widths and branching ratios measured by 
KLOE to date is shown in table XIV. Results of tests on conservation of discrete sym- 
metries, as well as on the unitarity of the CKM matrix are shown in table XV. 

It is interesting to note that over half of the results noted above are based on analysis 
of one fifth of the total KLOE data sample. Several of the measurements have reached 
the state where additional statistics will not improve accuracy due to theoretical or inter- 
mediate state uncertainties as described in the specific sections. There are measurements 
that could profit with respect to statistical and systematic errors continuing the studies of 
the entire body of data. Furthermore, several new measurements are only possible using 
the whole data set whose analysis is expected to be completed in the coming couple of 
years. These are discussed in the following three subsections. In the fourth subsection we 
describe a future project KLOE-2 that entails minor upgrade of KLOE in order to com- 
plete the original physics program proposed in 1992 but had been stymied by DA$NE's 
limited luminosity. Thus it is clear that we have our hands full for the foreseable future! 

111. Kaon sector. - 

1. Studies on rare kaon decays can benefit of the use of the entire KLOE statistics. 
The accuracy in the determination of the branching ratio of Ks — > Tt + e~v(^) events 
and of the related charge asymmetry will scale with statistics with respect to those 
shown in sect. 5'4.2, thus reaching ~0.5% (fractional) and 5xl0~ 3 respectively. 

2. We have observed the Ks — ► tt + ^D{"f) decay channel and will aim at reaching a 
precision better than 2% on its BR. 

3. We plan to improve our limit on Ks — * 37r° decays (see sect. 5'4.4). Preliminary 
studies have shown that, by slightly modifying the analysis cuts, one can improve 
on background rejection, while leaving the signal efficiency almost unaffected. This 
means that we can hope to set the limit on the corresponding BR at the level of 
~2xl0~ 8 . 

4. Moreover, we can search for Ks — > 7r + 7r~7r° events, which proceed mainly via a 
CP-conserving, AI = 3/2 transition. A preliminary, but rather detailed, analysis 
shows that KLOE should be able to measure the corresponding branching ratio 
with a ^60% precision, comparable to the precision with which this BR has been 
measured by three experiments. 

5. The CP violating Kl — > 2n° events can be measured with a ^0.5% accuracy, thus 
setting at the few per mil level the precision on Re(e'/e) which KLOE can reach. 

6. The study of the main Kj, decay channels, will allow us to improve our accuracy 
in the measurement of the Kl lifetime by a factor of <~2. We can also improve the 
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Masses (MeV) 








K° 


497.583 ± 0.021 




V 


547.874 ± 0.032 


Lifetimes (ns) 








K L 


50.84 ± 0.23 




K ± 


12.347 ± 0.030 


Widths (keV) 








4> — ► e + e~ 


1.32 ± 0.06 


Branching Ratios 






4> decays 


4> -> 7/7 


(6.20 ± 0.27) xl0~ 5 




4> -> 7r°7r°7 


(1.07 ± 0.07) xl0~ 4 




(f> — ► 777T°7 


(7.01 ± 0.24) xl0~ 5 




<t> -f W7T° 


(4.4 ± 0.6) xl0~ 5 




(j> -> K°X0 7 


< 1.8 xl0~ 8 at 90% C.L. 


Kl decays 


Jfz, — ► 7re;/ 


0.4008 ± 0.0015 






0.2699 ± 0.0014 




tfz, 3tt° 


0.1996 ± 0.0020 




Jfz, — ► tt + tt~tt" 


0.1261 ± 0.0011 




Kl — ► 7r + 7r~ 


(1.963 ± 0.21) xlO" 3 




Kz, -> 77 


(5.569 ±0.077) xl0~ 4 


Ks decays 


Ks — ► 7T + 7T~ 


0.60196 ± 0.00051 




K s -► 7r°7r° 


0.30687 ± 0.00051 






(7 4- n dq"i v 1 n -4 




Ks — * 77 


(2.26 ± 0.13)xl0~ 6 




Ks 3tt° 


<1.2xl0" 7 at 90% C.L. 




A"s — > e + e _ (7) 


<9.3xl0" 9 at 90% C.L. 


K ± decays 


K+ -f m + ^(t) 


0.6366 ± 0.0017 




K+ -» 7r+7r°(7) 


0.2067 ± 0.0012 




K+ -f 7r°e+z/(7) 


0.04972 ± 0.00053 




K+ -> 7r°/i + ^(7) 


0.03237 ± 0.00039 




K+ - TT+TtV 


0.01763 ± 0.00034 


r\ decays 


77 — ► -K + ir~e + e~ 


(2.4 ± 0.4) xl0~ 4 




?7 — > 7T°77 


(8.4 ± 3.0) xlO -8 


ui decays 


Ct) — » 7T + 7r~7r° 


0.9024 ± 0.0018 




CJ — > 7T°7 


0.0809 ± 0.0014 



Table XIV. - Summary of KLOE results on masses, lifetimes, widths and branching ratios 
measurements. 



precision in the measurement of the scalar form factor for semileptonic decays by 
the same amount. 

7. The form factor slopes for K ± can be measured as well, using about 5.9xl0 6 Kf 3 
and about 2.5 x 10 6 K^ 3 decays from the entire data sample. It is important to note 
that in the scalar form factor measurement from K ± there is no 7r//j ambiguity. 

8. The accuracy on the lifetime can be improved by a factor ~ 1.5. As discussed in 
sect. 6, these are key ingredients for the determination of the CKM matrix element 



PRECISION KAON AND HADRON PHYSICS WITH KLOE 



87 



C symmetry 


A LR (r) -> TT + n-TT°) 

A Q (r) -» 7T + 7r~7T ) 

A a ( 71 > 77^" 7T 7T^ ^ 


< 1.6 xlO" 5 
(Omtl lr) xl0~ 2 
(-0.05 ± 0.11) xlO" 2 
f0 08+S ^l xlO" 2 


CP symmetry 


|e| (/rom BR(K L -> Tr+Tr",) 

R,e(e) (from Bel I- Steinberg ev relation) 

t] — y n + iv~ 


(2.216 ± 0.013) )xl0~ 3 
(159.6 ± 1.3) xlO -5 
< 1.3 xlO" 5 


CPT symmetry 

with CP symmetry 


t(K+)/t(K-) 
M(o + ) - M(o~) 

r(p + )-r(p-) 

lm(5) (from Bell- Steinberg er relation) 
Re S + Re x- = (As - A L )/A 
Ree - Rey = (As + Az,)/4 


(1.004 ± 0.004) 
(1.5 ± 1.1) MeV 
(1.8 ± 2.1) MeV 
(0.4 ± 2.1) x 10 5 
(-0.5±2.5)x 10~ 3 
(1.2±2.5)x 10~ 3 


Lepton universality 


r Me (from /+(0) K s for A^ and 


1 1.000 ± 0.008 

(2.55 ±0.05 ±0.05) x 10 


CKM matrix unitarity 


/+(0)K. 

|K.| a + |V„df -1 (Kdfrom [87]) 


0.2157±0.0006 
0.2766± 0.0009 
-0.0004 ± 0.0007(~ 0.6a) 



Table XV. - Summary of KLOE results on symmetry tests. 



9. Semi-rare charged kaon decay channels can also be studied. In particular, we are 
planning to measure the branching ratio of the decay into three charged pions. With 
the present data sample a few per mil statistical accuracy can be reached. Using 
the same data set the measurement of the branching ratio of the — > 7r°7r°e ± i' e 
decay can reach a 10% statistical accuracy. This will allow us to complete the 
measurement of the main branching ratios. 

11'2. Light meson spectroscopy. - Many analyses on light meson spectroscopy could 
be refined in the future. Although the level of accuracy for the measurements on the 
scalar sector has already reached its limit, a combined analysis of different final states, 
such as (j) — ► ir + ir~j/(f> — > 7r°7r°7, will give more insight into the structure of light scalar 
mesons. 

On the other hand, the pseudoscalar mesons could benefit from the larger statistics 
of the whole KLOE data sample. 

1. In particular, we can complete the study of the main n decay channels by measuring 
the BR's of r/ — > ir + ir~j, 77 — ► e + e~j with a statistical precision of few per mil. 

2. Moreover, the knowledge on numerous rare n decay channels can be improved. 
Apart from lowering the limits on C (P,CP) violating decays 77 — > 3^/ (r/ — > 
7r + 7r~, 7r°7r°), we could reach a statistical precision of 3.5% on the BR(?7 — > 7r + 7r~ 
e+e~) and measure for the first time the asymmetry between the irir and ee decay 
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planes. The knowledge of the systematica for the latter decay could reach a few % 
accuracy by relying on large data/MC control samples. An observation at the level 
of 100 events is expected for the double Dalitz decay rj — > e+e~e+e~. 

3. Finally, the BR(?7 — > 7r°77) will reach a 1.5% statistical accuracy. Improvements of 
the analysis based on a new clustering scheme and on a photon veto at small angle 
could reduce the related systematics and background contamination, allowing also 
the measurement of the 77 invariant mass spectra. 

4. For the rj' , we have in hand about 450,000 of these mesons, so precise measurements 
could now be carried out. We plan to do this in three ways: 

i. looking at the 7/ dynamics by studying the Dalitz plot for the rj — > ?77r + 7r~, 
■q' — > 7r + 7T~7r° processes; 

ii. measuring ratio of branching fractions for some of the main decay channels, 
such as rj' — > ir + ir~rj, rj' — > ir°ir°rj, rj' — > W7, thus improving their BR accuracy; 

iii. searching for never heretofore observed decay channels, whose upper limits are 
at the % level (rj' — > 7r + 7r~7r°, rj' — ► 7r + 7r~7r + 7r~) 

11'3. Pzon /orm factor at about 1 GeV . - In early 2006, KLOE also collected ^200 
pb _1 of data at y/s = 1000 MeV, at which cr(e ± — > — > 7r + 7r~7r°) drops to 5% of its peak 
value. These data by themselves should determine the contribution to a™ normalizing 
to absolute luminosity, with a statistical error of ^35 x 10 -11 . This will give us another 
check on our measurements using the previously mentioned three methods: small angle 
selection, large angle selection and mr-f to jiri-f ratio. 

11'4. KLOE-2. - In late 2007 an important accelerator experiment has begun at 
DA$NE, with the goal of increasing the machine's luminosity by a factor of at least 
three, from the implementation of an innovative interaction scheme [142]. 

The exciting perspective of collecting much more data within a few years, inspired 
many KLOE members, together with several new collaborators, to form a new enterprise 
with the goal of continuing the KLOE physics program with an upgraded detector [143]. 

This project, dubbed with a certain lack of fantasy KLOE-2, foresees the insertion 
of new subdetectors in the inner part of KLOE, between the beryllium beam pipe and 
the internal wall of the drift chamber. An inner tracking chamber, based on the GEM 
technology is being developed at LNF. Crystal calorimeters will help catching photons 
emitted at very low polar angles. These detectors will increase acceptance for many 
physics channels, especially for Ks, rj and rj' decays. The implementation of a detector 
to tag 77 interactions is also under study. 

KLOE-2 is an ambitious project, which requires a few years of preparation and a 
longer period of data taking: the goal is to acquire at least 20-30 fb _1 by year 2012. 

Appendix A. The KLOE Collaboration: inclusive to 2008 
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